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Vorwort des Herausgebers 


Die Fahrzeugtechnik ist kontinuierlich Veränderungen unterworfen. 
Klimawandel, die Verknappung einiger für Fahrzeugbau und -betrieb 
benötigter Rohstoffe, globaler Wettbewerb, gesellschaftlicher Wandel und 
das rapide Wachstum großer Städte erfordern neue Mobilitätslösungen, die 
vielfach eine Neudefinition des Fahrzeugs erforderlich machen. Die 
Forderungen nach Steigerung der Energieeffizienz, Emissionsreduktion, 
erhöhter Fahr- und Arbeitssicherheit, Benutzerfreundlichkeit und 
angemessenen Kosten sowie die Möglichkeiten der Digitalisierung und 
Vernetzung finden ihre Antworten nicht aus der singulären Verbesserung 
einzelner technischer Elemente, sondern benötigen Systemverständnis und 
eine domänenübergreifende Optimierung der Lösungen. 


Hierzu will die Karlsruher Schriftenreihe für Fahrzeugsystemtechnik 
einen Beitrag leisten. Für die Fahrzeuggattungen Pkw, Nfz, Mobile 
Arbeitsmaschinen und Bahnfahrzeuge werden Forschungsarbeiten 
vorgestellt, die Fahrzeugsystemtechnik auf vier Ebenen beleuchten: das 
Fahrzeug als komplexes, digitalisiertes mechatronisches System, die 
Mensch-Fahrzeug-Interaktion, das Fahrzeug in Verkehr und Infrastruktur 
sowie das Fahrzeug in Gesellschaft und Umwelt. 


Die wahrnehmbare Akustik eines Fahrzeugs ist gleich nach dem 
optischen Eindruck für die Wertanmutung und Akzeptanz eines Fahrzeugs 
von großer Bedeutung. Hierzu tragen insbesondere der Fahrzeugantrieb und 
das Reifen-Fahrwerk-System bei. Dort werden Elastomerbauteile aufgrund 
ihres Schwingungsdämmungseffekts in großer Anzahl eingesetzt. 
Massenträgheiten auf dem Übertragungsweg von Schwingungsquellen 
zum Gehör der Fahrzeugnutzer führen zu einem Tiefpassfilterverhalten. 
Dies hat in der Vergangenheit dazu geführt, dass das hochfrequente 
Übertragungsverhalten von Elastomerlagern oberhalb von 300 Hz meist 
vernachlässigt wurde. Durch die fortschreitende Geräuschreduzierung von 
Verbrennungskraftmaschinen und insbesondere durch die vermehrte 
Verwendung elektrischer Antriebsmaschinen werden hochfrequente 
Geräuschanteille zunehmend deutlicher wahrnehmbar, so dass eine 
Auslegung von Elastomerlagern bis in den Kilohertzbereich 


Vorwort des Herausgebers 


wünschenswert ist. Hierzu fehlen aber bislang die Methoden und das 
Wissen über die relevanten Einfluss- und Gestaltungsparameter. Hier 
setzt Herr Watzl mit seiner Arbeit an, in der er am Beispiel von 
Aggregatelagern Vorgehensweisen für Versuch, Modellbildung, 
Simulation und die gezielte Auslegung der Designparameter für 
Gummi-Metall-Buchsen und Hydrolager aufzeigt und in deren 
Anwendung eine erhebliche Verbesserung des Dämmungsverhaltens 
erreicht. 


Karlsruhe, im März 2022 Frank Gauterin 


Kurzfassung 


Im Zuge der Reduktion von Emissionen eines Kraftfahrzeuges werden so- 
wohl innermotorische, wirkungsgradsteigernde Maßnahmen ergriffen als 
auch die Triebstränge zunehmend elektrifiziert. Diese Maßnahmen führen 
zu einem verstärkten Schwingungseintrag bei hohen Frequenzen und bedin- 
gen eine Verbesserung des Isolationsverhaltens des Körperschallpfades vom 
Antrieb bis zum Fahrerohr. Hierzu ist es notwendig die Körperschalltransfe- 
reigenschaften der Lagerelemente auf die im Gesamtfahrzeugverbund vor- 
herrschenden Randbedingungen anzupassen. Dies betrifft insbesondere die 
Körperschallempfindlichkeiten der anschließenden Karosseriestruktur und 
das Anregungsspektrum des Antriebsaggregats. 

In der vorliegenden Arbeit wird daher das Körperschallübertragungsver- 
halten von hydraulischen und konventionellen Aggregatelagerelementen im 
akustisch relevanten Frequenzbereich bis circa 2000 Hz experimentell und 
numerisch untersucht. Hierzu wird zunächst eine Prüfstandsumgebung zur 
Ermittlung der dynamischen Körperschalltransfereigenschaften im Kilo- 
hertzbereich vorgestellt und anschließend die wichtigsten Einflussparame- 
ter auf das Bauteilverhalten experimentell untersucht. In einer anschließen- 
den Modellbildung wird das Systemverhalten durch ein Finite-Elemente- 
Modell abgebildet, um die Eigenmoden des betrachteten Systems genauer 
zu analysieren. Durch dieses Vorgehen gelingt die Identifikation mehrerer 
charakteristischer Resonanzen, sowohl von Buchsenlagern, als auch von 
hydraulischen Topflagern, welche den einzelnen Komponenten der Bauteile 
zugeordnet werden können. 

Darüber hinaus kann durch die erfolgte Modellbildung ein gezieltes Design 
des Systems durchgeführt werden, um diese Resonanzen auf die jeweili- 
gen Fahrzeugrandbedingungen anzupassen. Im Anschluss an die Modell- 
bildung wurde zusätzlich eine Design-Space-Exploration für mögliche Ela- 
stomerkomponenten von hydraulischen Topflagern durchgeführt. Hierdurch 
können die möglichen Optimierungspotentiale sowie die Parameter, welche 
zur Hebung dieser notwendig sind, aufgedeckt werden. Durch die in der 
vorliegenden Arbeit herausgearbeitete Modellierung der Bauteile kann be- 
reits in einer frühen Fahrzeugentwicklungsphase eine genaue Aussage über 
das Körperschalltransferverhalten und somit über das zu erwartenden Ge- 
räusch am Fahrerohr getroffen werden. Hierdurch ist es in einer virtuellen 


ii 
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Fahrzeugentwicklung möglich Optimierungspotentiale bezüglich des An- 
triebsgeräuschkomforts abzuleiten und die akustische Grundauslegung von 
Lagerelementen ohne Gesamtfahrzeug- und Komponentenversuche durch- 
zuführen. 


Abstract 


Due to emission reduction of passenger cars, efficiency of combustion engi- 
nes is increased and powertrains are more and more electrified. In order to 
do so, measures are taken, that increase high frequency engine caused vibra- 
tions. Therefore, isolation of the noise transfer path from the engine to the 
pessenger’s ear needs to be increased. This makes it necessary to tune engi- 
ne mounts depending on the properties of the overall vehicle. This includes 
in particular the sensitivity to structure borne noise of the vehicle body and 
excitation spectrum of the engine. 

In the present thesis the structure-borne sound transmission properties of 
hydraulic and conventional engine mounts is analyzed experimentally and 
numerically in a acoustically relevant frequency range up to about 2000 Hz. 
In order to do so, first of all a test bench environment for analyzing the 
dynamic sound transmission properties is presented. In subsequent steps the 
main influence parameters on the transmission properties are determined ex- 
perimentally. Analysis of the system behavior and individual modes using a 
finite element model is conducted in a following model development step. 
Using this procedure allows the identification of several characteristical re- 
sonances of conventional bushings and hydraulic engine mounts, which can 
be allocated to specific components of the mounts. 

Forthermore, using the developed model leads to a specific design of engine 
mount systems fitting the resonances to the boundary condtions of the ve- 
hicle. Subsequent to the model development a design space exploration of 
the elastomeric components of hydraulic engine mounts is executed. Here- 
by, the possibilities for optimization and paremeters, that are necessary for 
enhancements, can be determined. Using the model of engine mounts de- 
veloped in the present thesis leads to prediction of noise transfer behavior 
and the estimated engine noise levels at the pessenger’s ears in early ve- 
hicle development process stages. This allows improvements on the engine 
noise level and a design of engine mounts in a virtual product development 
without or less vehicle or component testing. 
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1 Einleitung und Zielsetzung 


Um eine Einführung in die vorliegende Arbeit zu geben, soll zunächst in 
Abschnitt 1.1 ein motivierender Überblick über die Bedeutung von Aggre- 
gatelagerungen im Gesamtfahrzeug gegeben werden. Anschließend sollen 
in Abschnitt 1.2 die Herausforderungen bei der Auslegung von Aggregate- 
lagerelementen dargelegt werden, um anschließend die für die vorliegende 
Arbeit relevanten Forschungsfragen zu definieren und daraus in Abschnitt 
1.3 eine Zielsetzung abzuleiten. Im darauffolgenden Abschnitt 1.4 soll ab- 
schließend die Strukturierung der Arbeit erläutert werden. 


1.1 Bedeutung der Aggregatelagerung im 
Gesamtfahrzeug 


Neben der Formgebung bzw. dem Design eines Fahrzeuges stellt der vi- 
broakustische Gesamteindruck eines Fahrzeuges die wichtigste unmittel- 
bar erlebbare Eigenschaft eines Automobils dar [ZEF* 13]. Zusätzlich trägt 
das Geräusch, insbesondere das Motorgeräusch, zu einer Emotionalisie- 
rung des Fahrzeuges bei und stellt somit ein wichtiges markenprägendes 
Kaufkriterium für den Kunden dar. Neben der Luftschallabstrahlung eines 
Antriebstranges findet in einem Fahrzeug auch eine Körperschallübertra- 
gung vom Antrieb in den Fahrgastraum statt. Der Hauptkörperschallpfad 
besteht dabei aus den Lagern, welche den Triebstrang in der Karosserie an- 
binden. Zur Isolation von unerwünschtem Körperschall von Verbrennungs- 
geräuschen, Kettentrieben und Zahneingriffen, sind diese Lagerelemente 
als nachgiebige Strukturen aus Elastomermaterialien ausgebildet. Durch die 
Nachgiebigkeit dieser Aggregatelager entsteht ein weiteres schwingungs- 
fähiges Subsystem im Fahrzeug, welches durch die Fahrbahnanregung zu 
tieffrequenten Schwingungen angeregt wird. Neben der Isolation des Kör- 
perschalls haben die Aggregatelager daher zusätzlich die Aufgabe diese 
tieffrequenten Schwingungen zu dämpfen. Da der Triebstrang, im speziellen 
bei hochmotorisierten Fahrzeugen, einen nicht unerheblichen Teil der Fahr- 
zeugmasse ausmacht, hat dessen Lagerung im Fahrzeug einen großen Anteil 
am Komforteindruck eines Fahrzeugs und auch an den fahrdynamischen Ei- 
genschaften. Somit liefern Aggregatelagerelemente einen wichtigen Beitrag 
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sowohl zum tieffrequenten Fahrkomfort als auch zum vibroakustischen Ge- 
samteindruck eines Fahrzeuges [Gen10, ZEF* 13, PBBF10]. 

Die Anforderungen an die Körperschallisolationseigenschaften eines Ag- 
gregatelagers steigen durch zunehmende Elektrifizierung bzw. Hybridi- 
sierung der Antriebsstränge noch weiter an, da durch den Entfall der La- 
dungswechselgeräusche eines Verbrennungsmotors eine Maskierung uner- 
wünschter, triebstranginduzierter Geräusche ausbleibt. 


1.2 Problemstellung und Definition der 
Forschungsfragen 


Die Auslegung von Aggregatelagerelementen erfolgt mit Hilfe moderner 
MKS-Rechentools heute größtenteils virtuell und mit dem Ziel der Opti- 
mierung tieffrequenter, fahrkomfortrelevanter Kriterien. Die vibroakusti- 
schen Eigenschaften können erst in der Hardwarephase im Gesamtfahrzeug 
abgeprüft und abgestimmt werden. Eine zielgerichtete Entwicklung bezüg- 
lich der Isolationseigenschaften mit größeren Änderungen an den Bauteilen 
ist aus Zeit- und Kostengründen dann meist nicht mehr möglich. Oftmals 
kommen zur Lösung spezifischer Probleme dann Tilger oder Sperrmas- 
sen zum Einsatz, welche Kosten- und Gewichtsnachteile für das Gesamt- 
fahrzeug mit sich bringen. Eine Auslegung der Bauteile in einer isolier- 
ten Prüfstandsumgebung oder einer virtuellen Umgebung kann hierbei das 
notwendige Frontloading in der Entwicklung bereitstellen, um auch die vi- 
broakustischen Anforderungen an Aggregatelagerelemente in den moder- 
nen Fahrzeugentwicklungsprozess einfließen zu lassen. Hierzu ist es jedoch 
notwendig, die wesentlichen Einflussparameter auf die Körperschallüber- 
tragungseigenschaften von Aggregatelagern zu kennen und diese gezielt 
auslegen zu können. Akustisch relevante Phänomene im Fahrzeug, welche 
durch Körperschall über die Aggregatelager übertragen werden, befinden 
sich üblicherweise in einem Frequenzbereich zwischen 200 Hz und 2kHz 
[PBBF10]. Die Geräuschanteile mit Frequenzen unterhalb dieses Bereichs 
können durch das Hörorgan nur schlecht wahrgenommen werden, diejeni- 
gen höherer Frequenzen werden durch sich mit steigender Frequenz verbes- 
sernde Karosserieeigenschaften deutlich geringer an das Insassenohr über- 
tragen. 


Ausgehend von der beschriebenen Problemstellung ergeben sich daher im 
Wesentlichen drei Forschungsfragen: 
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e Wie lassen sich die Körperschalltransfereigenschaften von Aggrega- 
telagern beeinflussen, um diese an die jeweiligen Anforderungen des 
Gesamtfahrzeuges anpassen zu können? 


e Welche Wechselwirkungen zwischen den Komponenten beeinflussen 
das Gesamtsystem? 


e Welche Grenzen für die Transfereigenschaften ergeben sich aus den 
Randbedingungen, die sich aus dem Einsatz der Bauteile in einem 
Kraftfahrzeug ableiten. 


1.3 Zielsetzung der Arbeit 


Aus diesen drei im vorherigen Abschnitt dargelegten Forschungsfragen 
lassen sich konkrete Zielsetzungen für die vorliegende Arbeit ableiten. 
Zum einen soll ein Verfahren erarbeitet werden, mit dessen Hilfe die Kör- 
perschalltransfereigenschaften der betrachteten Bauteile experimentell er- 
mittelt werden können. Im Anschluss daran sollen über Parameterstudien 
an verschiedenen Bauformen von Aggregatelagern die wesentlichen Ein- 
flussparameter auf das Körperschallübertragungsverhalten herausgearbeitet 
werden. Aufbauend auf diesen Parameterstudien soll eine Modellbildung er- 
folgen, welche das Körperschallübertragungsverhalten der Bauteile abbildet 
und die erarbeiteten Parameter berücksichtigt. Mit diesem Modell sollen im 
Anschluss mit Hilfe einer multikriteriellen Optimierung die Grenzen einer 
möglichen Auslegung der Bauteile aufgezeigt werden. Auf diesem Wege 
kann ein Auslegungsspielraum dargelegt werden, in welchem eine Optimie- 
rung der Körperschallübertragung durch die Aggregatelager erfolgen kann. 


1.4 Strukturierung der Arbeit 


Die Struktur der vorliegenden Arbeit orientiert sich an der im vorherigen 
Abschnitt dargelegten Zielsetzung. Zunächst soll in Kapitel 2 der Stand der 
Forschung und Technik zur Modellierung der Körperschalltransferpfade, 
zum dynamischen Verhalten von Elastomeren sowie zum Aufbau und der 
Herstellung von Aggregatelagern erörtert werden. Im darauf folgenden Ka- 
pitel 3 soll anschließend eine Methodik aufgezeigt werden, mit deren Hilfe 
die Körperschalltransfereigenschaften der Lagerelemente experimentell er- 
mittelt werden können. Darauf aufbauend werden im selben Kapitel die 
Parameterstudien an den Bauteilen erläutert und die wesentlichen Einflus- 
sparameter herausgearbeitet. 


1 Einleitung und Zielsetzung 


In Kapitel 4 soll das Körperschalltransferverhalten der Bauteile dann durch 
geeignete Modellierung nachgebildet werden. Hierzu wird zunächst eine 
kurze Einführung in die verwendete Modellierungsstrategie gegeben. Hier- 
zu wird an dieser Stelle nochmals auf die theoretischen Grundlagen einge- 
gangen. Da diese nur für das Kapitel der Modellbildung eine Rolle spielt, ist 
dieser Teil nicht im Stand der Forschung und Technik integriert. Anschlie- 
Bend wird die Modellierung und eine entsprechende Validierung dieser Mo- 
dellierung durchgeführt. 

Kapitel 5 führt, aufbauend auf der vorangegangenen Modellierung, eine 
Design-Raum-Exploration durch, um die Grenzen der Auslegung für Aggre- 
gatelagerbauteile aufzeigen zu können. Auch hier wird, wie im vorherigen 
Kapitel, nochmals eine kurze Einführung in die Grundlagen des verwende- 
ten Optimierungsalgorithmus gegeben. 

Abschließend soll in Kapitel 6 eine Schlussfolgerung und Zusammenfas- 
sung der Ergebnisse dieser Arbeit gegeben werden. Darüber hinaus soll 
ebenso gezeigt werden, inwieweit die in dieser Arbeit gewonnenen Erkennt- 
nisse den Stand der Technik erweitert haben. 
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Im folgenden Kapitel sollen zunächst die für die Bearbeitung der gestellten 
Forschungsfragen notwendigen Grundlagen dargelegt werden. Hierzu soll 
in einem ersten Schritt eine analytische Betrachtung des Körperschalltrans- 
ferpfades mit Hilfe von Vierpolen vorgestellt werden (Abschnitt 2.1). An- 
schließend soll in Abschnitt 2.2 auf die spezifischen Werkstoffeigenschaften 
der Elastomere eingegangen werden, welche für die betrachteten Bauteile 
von Relevanz sind. Daraufhin wird in Abschnitt 2.3 ein Überblick über den 
Aufbau verschiedener Aggregatelagerelemente und die Herstellung dieser 
gegeben. Zuletzt wird auf Basis der dargelegten Grundlagen die in Abschnitt 
1.2 formulierten Forschungsfragen konkretisiert. 


2.1 Modellierung von Transferpfaden und 
Transfereigenschaften 


Zur Betrachtung der Körperschallübertragung in Fahrzeugen soll in den fol- 
genden Abschnitten eine theoretische Einführung in die Beschreibung des 
dynamischen Verhaltens von Bauteilen und deren Schnittstellen gegeben 
werden. Aus der Betrachtung dieser lassen sich verschiedene Dämmungs- 
größen ableiten, um die Wirksamkeit von Körperschallisolatoren beurteilen 
zu können. Mit dem Hintergrund der Optimierung des Komforts und der 
Innenraumakustik in Fahrzeugen wird detailliert auf die Beschreibung ein- 
zelner Bauteile mittels mechanischer Vierpole eingegangen. 


2.1.1 Mechanische Kenngrößen der Transfereigenschaften von 
Bauteilen 


Die Beschreibung der Körperschalltransfereigenschaften eines Bauteils kann 
mit Hilfe verschiedener mechanischer Kenngrößen erfolgen. 

Wird ein Körper mit einer punktförmigen, harmonischen Kraft der Form 
F(t) = Fsin(@ t) beaufschlagt, so reagiert dieser mit einer resultierenden 
Deformation x(t) = ê sin(@t + ô). Hierbei sind F(t) das Zeitsignal der 
Kraft, F die Amplitude des harmonischen Kraftverlaufs, x(t) das Zeitsignal 
des Schwingweges und ĉ die Amplitude des harmonischen Schwingweges. 
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Durch eine komplexe Erweiterung der Sinusterme kann die Kraft bzw. De- 
formation des Körpers als komplexe Größe ausgedrückt werden zu: 


F(t) =F d” bzw. x(t) =£ e° de® (2.1) 
mit: dem Winkelversatz 6 zwischen dem Kraft- und Wegsignal. 


Um das Übertragungsverhalten zu charakterisieren, wird die Reaktion des 
Bauteils auf die jeweilige Anregung bezogen, die diese Reaktion zur Fol- 
ge hat. Dieser Zusammenhang wird im Zeitbereich durch Differentialglei- 
chung(-en) beschrieben und lässt sich im Frequenzbereich durch die Über- 
tragungsfunktion T(®) gemäß Gleichung 2.2 beschreiben [ZEF* 13]. 


A(@) 
H(@) = ——~ 2.2 
Hlo) = Fa (2.2) 
mit: A(@) und E(@) die Laplace- bzw. Fourier-Transformierte der 
Zeitsignale A(t) am Ausgang und E(t) am Eingang. 


Die gebräuchlichste Übertragungsfunktion von Bauteilen ist die mechani- 
sche Steifigkeit. Diese ist definiert als der Widerstand eines linear elasti- 
schen Körpers gegen eine mechanische Verformung bei einer einwirkenden 
Kraft (vgl. [KünO1]) und lässt sich somit aus den komplexen Größen der 
Kraft und der Deformation errechnen zu: 


c(@) = — (2.3) 


mit: c(@): komplexe dynamische Steifigkeit, F(@): Fourier-Transformierte 
des Kraftverlaufs, x(@): Fourier-Transformierte des Schwingweges, 


Betrachtet man die Deformation des Körpers an einer Stelle j bei einer 
Krafteinleitung an der Stelle i, so ergibt sich die Transfersteifigkeit c’/des 
Körpers: 


dÄ (o) = = (2.4) 


mit: c'/(@): komplexe Transfersteifigkeit vom Punkt i zum Punkt j, F ‘(@): 
Fourier-Transformierte des Kraftverlaufs an der Stelle i, x/(@): 
Fourier-Transformierte des Schwingweges an der Stelle j. 


Sind der Ort der Krafteinleitung und der Ort der Deformation identisch, so 
spricht man von der Ein- bzw. Ausgangssteifigkeit eines Körpers. 
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Wird statt des Schwingweges x an der Stelle j die Schwingschnelle v be- 
trachtet, so kann als Beschreibungsgröße für die Widerstandsfähigkeit des 
Körpers gegen eine äußere Erregung die mechanische Impedanz herangezo- 
gen werden (vgl. [KM85]). Diese errechnet sich demnach gemäß: 

F‘(®) 


ZU(o) = wo) (2.5) 


mit: Z‘/(@): komplexe Transferimpedanz vom Punkt i zum Punkt j, F! (œ): 
Fourier-Transformierte des Kraftverlaufs an der Stelle i, v/(@): 
Fourier-Transformierte der Schwingschnelle an der Stelle j. 


Neben der Steifigkeit und der Impedanz gibt es noch weitere Darstellungs- 
formen der Übertragungseigenschaften von Bauteilen. Die Nomenklatur 
dieser Übertragungsfunktionen entstammt zum Teil der Elektrotechnik und 
ist aus der Betrachtung von elektromechanischen Analogien entstanden. 
Eine Übersicht über die wichtigsten dieser Übertragungsfunktionen ist in 
Tabelle 2.1 dargestellt. 

Die Übertragungsfunktionen von idealisierten Bauteilen lassen sich gemäß 
den in Tabelle 2.1 dargestellten Beziehungen einfach bestimmen. Für ei- 
ne ideale Punktmasse ist die Kraft nach dem zweiten Newtonschen Axiom 
proportional zur Beschleunigung a, für einen idealen viskosen Dämpfer 
proportional zur Schwingschnelle v und für eine ideale, masselose Feder 
proportional zum Schwingweg x. Demnach ergeben sich die Reaktionskräf- 
te dieser idealen Komponenten mit den Konstanten m, k und c zu den in 
Tabelle 2.2 gegebenen Zusammenhängen. 

Über den Zusammenhang im Frequenzbereich zwischen Beschleunigung, 
Schwingschnelle und Schwingweg 


a(@) = jav(@) = (j@)’x(®) (2.6) 


mit: a(@): Fourier-Transformierte der Schwingbeschleunigung, a(@): 
Fourier-Transformierte der Schwingschnelle, x(@): Fourier-Transformierte 
des Schwingweges 


ergeben sich die Ubertragungsfunktionen dieser idealisierten Bauteile im 
Frequenzbereich gemäß der in Tabelle 2.3 aufgeführten Zusammenhänge. 
Die Übertragungsfunktionen dieser idealisierten Bauteile sind in Abbildung 
2.1 dargestellt. 

Die elektromechanischen Analogien erlauben es ebenfalls, die Größen, wel- 
che für die schwingungstechnische Betrachtung von Körperschallübertra- 
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kinematische 
Größe Übertragungs- : 
Berechnungs- funktion Berne 
grundlage 
i F(o) dynamische Steifigkeit, 
Schwingweg xo) =C(o) Hemmung , effektive Steife 
x 
x(@) _ dynamische Nachgiebigkeit, 
Fa) ~N()  Weichheit 
F(o) mech. Impedanz, mech. 
Schwingge- vo) — Z(®) Scheinwiderstand, Mitgang 
schwindigkeit, 
Schnelle Re{Z(a)} Resistanz 
Seu) Im{Z(@)} Reaktanz 
xo) _y (0) mech. Admittanz, Beweg- 
Flo) >= lichkeit, Mobilität 
Schwingbe- F(o) _y (o) dynamische Masse, Schein- 
schleunigung a(o) = masse, effektive Masse 
k=v=a a = A(@) Akzeleranz, Inertanz 


Tabelle 2.1: Nomenklatur der komplexen Ubertragungsfunktionen von mechani- 
schen Systemen [ZEF* 13, KM85, SEM08] 


gungspfaden notwendig sind, mit Hilfe von Beziehungen der Elektrotech- 
nik und Systemtheorie zu betrachten. Hierdurch können die Methoden und 
Werkzeuge, die in der elektrischen Schaltungsanalyse entwickelt wurden, 
für die Analyse von mechanischen Schwingungssystemen angewandt wer- 
den (vgl. [ZEF*13]). Im folgenden Abschnitt soll daher explizit auf die 
Beschreibung von Bauteilen mittels Vierpolen eingegangen werden. 


2.1.2 mechanische Vierpole 


In der Systemtheorie können mechanische oder elektromechanische Syste- 
me durch Zweitore bzw. Vierpole beschrieben werden. Die Ein- bzw. Aus- 
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Komponente Reaktionskraft 
ideale Punktmasse F(@) =ma(®) 


idealer viskoser Dämpfer F(@) =kv(q@) 
ideale Feder F(®)=cs(®) 


Tabelle 2.2: Reaktionskräfte idealer Bauteile 


idealer viskoser 


ideale Punktmasse Dämpfer (Dim- 


masselose Feder 


(Masse m) pferkonst.&) (Federkonst. c) 
Bir. fe EN oe E _— c 
dyn. Masse M am m u Ga Gey 
(j@)F _ B —_ E- č 
Impedanz Z m = jom ak = ja 
dyn. Steifig- io)? F i jo)E _.; 
mic aloe Sek Zee 


Tabelle 2.3: Übertragungsfunktionen idealisierter Bauteile [Sch79, Se105] 


gangsgröße dieser Vierpole sind die Kräfte vor und nach dem jeweiligen 
Element und je nach Betrachtungsform der Schwingweg bzw. die Schwing- 
schnelle am Ein- und Ausgang. Die Anwendung der Vierpol-Theorie auf 
mechanische Schwingsysteme erfolgte 1957 erstmalig durch Molloy [M0157]. 
Snowdon [Sno79], Sell [SEMOS, Sel05], Tesar [TF88] und Waller [WK75] 
geben ebenso eine Herleitung und verschiedene Anwendungsfälle auch für 
die Körperschallübertragung. 


Grundlagen der Theorie mechanischer Vierpole 


Vierpolparameter können in Kettenform dargestellt werden (vgl. [SEMO08], 
[Mol57]), bei der die Ausgangsgrößen in Abhängigkeit der Eingangsgrößen 
angegeben sind. Alternativ dazu kann auch die Steifigkeitsform oder die Im- 
pedanzform (vgl. [SEM08, Ger80]) benutzt werden. 

In dieser Darstellung stellen die Vierpolparameter Steifigkeiten dar und las- 
sen sich so einfach interpretieren. In den nachfolgenden Abschnitten und 
Kapiteln wird daher explizit auf die Darstellung in Steifigkeitsform einge- 
gangen. Abbildung 2.2 zeigt einen solchen mechanischen Vierpol in Steifig- 
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a) b) 
eek 11 o rer 
© NI ~. D Team 
N . >. 
N ~ =~ 
N 
N 
N 
N 
1 10 1001000 1 


10 1001000 1 10 100 1000 
Frequenz in Hz Frequenz in Hz Frequenz in Hz 
—— Punktmasse = = = - viskoser Dämpfer =-=-=-- Feder 


a):dynamische Masse; b):Impedanz; c):dynamische Steifigkeit 


Abbildung 2.1: Graphische Darstellung mechanischer Transferfunktionen ideali- 
sierter Bauteile [Sel05] 


keitsform. 


E 
iS 


} 


mit: F;: Kraft an der Eingangsseite des Vierpols; F: Kraft an der 
Ausgangsseite des Vierpols; x]: Schwingweg an der Eingangsseite des 


Vierpols; x2: Schwingweg an der Ausgangsseite des Vierpols; C: Vierpol 
Matrix des Vierpols in Steifigkeitsform 


Abbildung 2.2: Vierpol eines mechanischen Bauteils in Steifigkeitsform [Sel05] 


Die zugehörige Matrix-Darstellung C des Vierpols in Steifigkeitsform kann 
gegeben werden durch: 


(2.7) 


mit: c;;:Transfersteifigkeit des Vierpols; cı: Eingangssteifigkeit des 
Vierpols; c>,: Ausgangssteifigkeit des Vierpols 
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Für die Berechnung der Übertragungseigenschaften von komplexen Syste- 
men, die aus mehreren hintereinandergeschalteten Vierpolen bestehen, ist 
hingegen die Kettenform geeigneter. Eine Darstellung des mechanischen 
Vierpols in Kettenform ist in Abbildung 2.3 gegeben. 

Die mathematische Beschreibung des Vierpols in Kettenform ist gegeben 


Fy Le; 
[A] 

x 2 

ro os 


mit: F: Kraft an der Eingangsseite des Vierpols, welche von außen auf den 
Vierpol wirkt; F,: Kraft an der Ausgangsseite des Vierpols, welche vom 
Vierpol auf die dort angebundene Struktur wirkt; x): Schwingweg an der 
Eingangsseite des Vierpols; x,: Schwingweg an der Ausgangsseite des 
Vierpols; A: Vierpol Matrix des Vierpols in Kettenform 


Abbildung 2.3: Vierpol eines mechanischen Bauteils in Kettenform [Sel05] 


durch: 


F F F 
fy =A 2| _ a&i] 412 2 (2.8) 


xX] X2 49; an| |X% 


mit: A der Vierpol-Matrix eines Vierpols in Kettenform 


Bei der Darstellung in Kettenform ist besonders zu beachten, dass die Kraft 
F diejenige Kraft ist, welche vom Vierpol auf die dort angebundene Struk- 
tur wirkt (vgl. [Sel05, Fre80]). Die Kraft Fı hingegen ist diejenige Kraft, 
welche von der am Eingang angebundenen Struktur auf den Vierpol wirkt. 
Dadurch wird erreicht, dass bei identischer Richtung der Kräfte gemäß Ab- 
bildung 2.3 die Ausgangsgrößen des vorherigen Vierpols den Eingangs- 
größen des darauffolgenden Vierpols entsprechen. Hierdurch wird die Be- 
rechnung dieser Schaltung vereinfacht. Neben der Steifigkeitsform und der 
Kettenform existieren noch weitere Darstellungsformen, die mit Hilfe von 
Transformationsmatrizen ineinander umgerechnet werden können [Rub67]. 
Die Umrechnung von der Steifigkeits- in die Kettenform ist in Tabelle 2.4 
dargestellt. 


Kopplung mechanischer Vierpole 


Ein mechanisches System kann als ein Netzwerk aus mehreren Teilsy- 
stemen, welche miteinander gekoppelt sind, betrachtet werden. Um diese 
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Kettenform Steifigkeitsform 
Fy aj; an| |E2 Fy Cy En] |x 
Gl = e= 
x Ay} an| | Xo Fy Co Lo9} | Xo 
2 F c det(C) 
A= 11 212 ©] er 
= 1 29 
a9, an Di Or 
det(A) 
an = 
C= a9] a] Cil €12 
= I a2 
ay ay C2) £22 


Tabelle 2.4: Umrechnung der Vierpolmatrizen von Steifigkeits- in Kettenform 


Kopplung im Verhalten des Gesamtsystems zu beriicksichtigen, kann jedes 
seiner Teilsysteme durch einen Vierpol beschrieben werden. Im Anschluss 
daran kann das System mathematisch zu einem Gesamt-Vierpol zusammen- 
gefasst werden [Sel05]. Die einzelnen Vierpole müssen, soweit sie mit der 
selben Schwingschnelle bzw. dem selben Schwingweg beaufschlagt wer- 
den, parallel geschaltet werden. Dies ist nur dann sinnvoll, wenn die ein- 
zelnen Elemente räumlich nicht zu weit voneinander entfernt angeordnet 
sind. Bei der Lagerung einer großen Maschine oder Anlage gilt dies nur, 
solange einfache Hubschwingungen vorliegen. Die Auswirkung von Kipp- 
und Drehbewegungen, bei denen die Schwingwege der einzelnen Lagerele- 
mente unterschiedlich sind, lassen sich mit Hilfe der Vierpolmethode nicht 
berechnen [Sei99, Sel05]. 


Werden mehrere Teilsysteme als starr miteinander verbunden betrachtet und 
wird davon ausgegangen, dass diese alle mit dem selben Schwingweg be- 
aufschlagt werden und die Ein- bzw. Ausgangskräfte der Teilsysteme sich 
zu den jeweiligen Gesamtkräften addieren, spricht man von einer Parallel- 
schaltung der Vierpole. Zur Bestimmung der Gesamtsteifigkeitsmatrix von 
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n parallelen Teilvierpolen werden die einzelnen Steifigkeitsmatrizen gemäß 
Gleichung 2.9 aufsummiert [Fel62]: 


N 
c= JCI (2.9) 
i=0 


mit: C^: Summensteifigkeits vierpol-Matrix der Parallelschaltung, (e : 
Steifigkeitsvierpol-Matrix des i-ten Vierpols der Parallelschaltung 


Bei einer seriellen Schaltung von mehreren Vierpolen, bei der die Ausgangs- 
größen des vorherigen Vierpols jeweils die Eingangsgrößen des nachfolgen- 
den sind, kann eine mathematische Beschreibung des Gesamtvierpols ma- 
thematisch am einfachsten über die Kettenform dargestellt werden [Fel62]: 


N 
ass =] JA‘ (2.10) 
i=0 


mit: AS“: Summenkettenvierpol-Matrix der Reihenschaltung, € : 
Kettenvierpol-Matrix des i-ten Vierpols der Reihenschaltung 


Vierpolparameter idealer mechanischer Elemente 


Nachdem in Abschnitt 2.1.1 bereits die Transfereigenschaften von den idea- 
lisierten Bauteilen Punktmasse, viskoser Dämpfer und Feder erörtert wur- 
den, sollen im folgenden Abschnitt auch die Vierpolparameter dieser Bau- 
teile hergeleitet werden. 


Punktmasse: Wird die Punktmasse mit der Vorzeichenkonvention des 
Vierpols in Kettenform aus Abbildung 2.3 freigeschnitten und vorausge- 
setzt, dass die Verschiebungen vor und nach der Punktmasse identisch sind, 
so folgt nach dem zweiten Newtonschen Axiom: 


mx; = F + Fo (2.11) 


xX, =X. (2.12) 
Ein Umstellen von Gleichung 2.11 ergibt das Gleichungspaar: 


Fı=R+mi = Po +m(jo)*x. (2.13) 
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xj Hx (2.14) 


Somit ergibt sich die Vierpolmatrix in Kettenform zu: 


F F 1 (jo)’m| |E 

= = AMasse 22 = (j ) £2 (2.15) 
aj X2 0 1 % 

Über die in Tabelle 2.4 gegebenen Transformationsvorschriften lassen sich 

die Vierpolmatrizen in die Steifigkeitsform überführen. Da der Eintrag a5, 

Null ist, ergibt sich die Steifigkeitsform zu: 


er = = (2. 1 6) 


Was Gleichung 2.16 mathematisch zum Ausdruck bringt, wird bei nähe- 
rer Betrachtung auch anschaulich logisch: Versucht man eine Punktmasse 
gegen eine unendlich steife Wand zu drücken, so wird sich diese nicht ver- 
formen. Demnach werden die Steifigkeiten unendlich groß. 


viskoser Dämpfer: Für den idealen, viskosen Dämpfer ergibt sich mit 
der Dämpfungskonstanten k die Reaktionskraft des Dämpfers mit der Vor- 
zeichenkonvention aus Abbildung 2.3: 

F 


Fi = F =k — %2) > x1 = + (2.17) 
jok 


Zusammen mit dem Kräftegleichgewicht (Fi=F)) ergeben sich die Vierpol- 
parameter für den idealen viskosen Dämpfer mit der Dämpferkonstanten k 
in Kettenform zu: 


Ei zum Fo = 1 0 2 (2.18) 


1 
x % Tok % 


| 


In Steifigkeitsform ergeben sich die Vierpolparameter zu: 


a jak —j@k 
Ce = A a (2. 19) 
T jok —jøæk 
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Masselose Feder: Analog zum Vorgehen für den viskosen Dämpfer er- 
geben sich aus dem Kräftegleichgewicht und der Reaktionskraft der idealen 
Feder mit der Federkonstanten c zwei Gleichungen: 


Fi =F, (2.20) 


F, 
Fy =c(x —x2) > x1 = = +x (2.21) 
c 


und damit die Vierpolparameter in Kettenform: 


F F 1 0| |F 
al ~ |x : Xp 
Diese lassen sich wiederum in die Steifigkeitsform umrechnen: 
: cc 
Ce = (2.23) 
Ze cc 


Für alle idealisierten Bauteile entsprechen die Eingangs- und Ausgangsstei- 
figkeiten gleich den Transfersteifigkeiten. Für die ideale Feder und den vis- 
kosen Dämpfer gilt dies nur, solange die Masse der Bauteile keinen Einfluss 
auf deren Verhalten zeigt. 


Beschreibung des Körperschallpfades im Fahrzeug mit Hilfe 
von Vierpolen 


Im Fahrzeug wird ein Körperschallpfad vom Aggregat über das nachgiebige 
Aggregatelager bis in die angrenzende Karosseriestruktur gebildet, von wo 
aus eine Abstrahlung als Luftschall stattfindet, welcher von den Passagieren 
wahrgenommen werden kann. Maßnahmen zur Minderung dieses Körper- 
schalleintrages zielen entweder darauf ab, die Emission von Körperschall 
durch die Quelle zu reduzieren oder die Schwingungsenergie zur Quelle zu- 
rückzureflektieren, um bei einem mehrmaligen Durchlaufen von schallab- 
sorbierenden Strukturen entsprechend reduziert zu werden. Nach Kurze und 
Sell [Sel05, KM85] kommen dabei vor allem schwere Massen oder weiche 
Federn zum Einsatz. Um die Dämmungswirkung beschreiben zu können, 
werden bevorzugt zwei verschiedene Arten der Dämmungsgrößen verwen- 
det. Die Durchgangsdämmung betrachtet das Verhältnis von Kräften oder 
Schwingwegen, -schnellen oder -beschleunigungen ins Verhältnis zueinan- 
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der. Die Einfügedämmung vergleicht die Schwingungsgrößen auf der Emp- 
fängerseite der Struktur beim Einsatz eines Isolationselementes mit dem 
Zustand ohne dieses Isolationselement. Bei dieser Beschreibung der Däm- 
mung von Körperschall kann auch die Auswirkung des Isolationselemen- 
tes auf das Verhalten des Gesamtsystems beurteilt werden. Beide Betrach- 
tungen erlauben jedoch keine quantitativen Rückschlüsse auf die Auswir- 
kungen des Isolationselementes auf den körperschallinduzierten Luftschall 
am Fahrerohr ohne dass ein Referenzzustand bekannt wäre. Da Fahrzeug- 
entwicklungsprozesse zunehmend virtualisiert werden, besteht das Ziel, das 
Antriebsgeräusch bereits auf Basis von Simulation oder Komponentenver- 
suchen beschreiben zu können 

Die Modellierung des Körperschallübertragungspfades durch Vierpole er- 
möglicht es, die Wechselwirkungen zwischen Triebstrang, Karosseriestruk- 
tur und den Aggregatelagern besser zu verstehen und zu quantifizieren. An- 
sätze hierzu wurden bereits von Fülbier [Fiil01], Sell [Sel01, Sel05] und 
Krämer & Helber [KH93] unternommen. Im Folgenden soll diese Methode 
dargestellt werden, um herauszustellen, welche Vierpolparameter der Ag- 
gregatelager von besonderem Interesse sind. Dieser Körperschallpfad vom 
Aggregat bis zum Insassenohr ist in Abbildung 2.4 dargestellt. Für jedes Ag- 


Abbildung 2.4: Vierpole für Lager und Karosserie eines Körperschallpfades 
[Gen10] 


gregatelager und jede Raumrichtung existiert je ein solcher Übertragungs- 
pfad, welcher einen Teilschalldruck erzeugt. Alle Teilschalldrücke wieder- 
um summieren sich zum Gesamtschalldruck im Innenraum auf. Die Trans- 
ferpfadanalyse mit Hilfe von Vierpolen setzt voraus, dass die Übertragung 
in den drei Raumrichtungen sich nicht gegenseitig beeinflussen und die Sy- 
steme als linear angenommen werden können [Fre80]. Für die Modellie- 
rung des Aggregatelagers wird zunächst die Steifigkeitsform verwendet, da 
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die Vierpolparameter so einfach zu bestimmen und interpretieren sind. Dies 
führt zu folgender Vierpolgleichung: 


= laser Xx] ch, ch, X (2.24) 
Fol x} |5 ch] ao 


Hierbei stellen ĉii die Eingangssteifigkeit des Lagers, gb die Ausgangsstei- 
figkeit des Lagers und ci bzw. ci die Vorwärts- und Rückwärts-Transfer- 
steifigkeiten des Lagers dar. Der in Abbildung 2.4 dargestellte Vierpol für 
die Karosseriestruktur zeigt den Zusammenhang zwischen mechanischen 
Eingangsgrößen in die Karosserie (F,,x,) und den akustischen Ausgangs- 
größen (P39). Hierbei stehen F) und x7 fiir die Kraft und den Weg, mit 
der die Karosseriestruktur durch das Lagerelement angeregt wird. p, ist der 
Schalldruck am Empfängerort, q, der Schallfluss durch die Oberfläche eines 
abgeschlossenen Volumens, welches den Empfängerort umgibt. Für jeden 
Körperschallübertragungspfad entsteht so ein Karosserie-Vierpol, welcher 
die Übertragungseigenschaften der Struktur von der jeweiligen Einleitungs- 
stelle zum Empfängerort abbildet. Krämer und Helber [KH93] sowie Sell 
[Gen10, Sel05] geben diesen Vierpol in Hybridform an, da so die für die 
Beschreibung notwendigen Größen einfach zu bestimmen sind. Bei hybri- 
den Vierpolen sind die Größen auf der Empfänger- und der Senderseite des 
Vierpols nicht getrennt sondern je eine links bzw. rechts des Vierpols ange- 
ordnet. Dementsprechend ergibt sich die Vierpolgleichung zu: 


K K 
P3 = Karosserie d; = hii hiz d3 (2.25) 
w F, hd, hh] |F> 


mit: HKarosserie, der Vierpolmatrix der Karosserie in Hybridform mit den 
Einträgen hy „ha: 151, h22, p,: dem Teilschalldruck an einem 
Empfängerort im Fahrzeuginneren, qz: dem Schallfluss durch die 
Oberfläche eines geschlossenen Volumens, welches den jeweiligen 
Empfängerort beinhaltet, x,: dem Schwingweg am Punkt der 
Körperschalleinleitung in die Karosserie, F,: der dynamischen Kraft am 
Punkt der Körperschalleinleitung in die Karosserie 


Der Schallfluss ist dabei die definiert durch [ZEF* 13]: 


q= [vda = [yaa (2.26) 
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mit: v: dem Schallschnellefeld; dA: vektorielles Oberflächenelement 


Aus Gleichung 2.26 ist ersichtlich, dass der Schallfluss die über eine Ober- 
fläche integrierte Schallschnelle ist. Hierbei ist hervorzuheben, dass nur der- 
jenige Anteil der Schallschnelle berücksichtigt wird, welcher senkrecht zur 
Oberfläche steht. Dies wird durch die vektorwertige Darstellung des Ober- 
flächenelements und die Verrechnung als Skalarprodukt mit dem Schallfluss 
berücksichtigt. Für die Ermittlung des Schallflusses q, wird das Volumen, 
welches den Empfängerort umgibt, klein gegenüber der betrachteten Wel- 
lenlänge gewählt. Ist dieses Volumen zusätzlich noch frei von einer Schall- 
quelle, so wird 4, sehr klein werden. Vereinfachend wird daher im Folgen- 
den q, = 0 gesetzt. 


Die Vierpolparameter h%, und h%, des Karosserievierpols sind besonders re- 
levant für die Beschreibung des Transferpfades: h%, wird oft als vibroakusti- 
sche Übertragungsfunktion bzw. Noise-Transfer-Function (NTF) bezeichnet 
[Fre92]. Sie gibt abhängig von der karosserieerregenden Kraft an, wie viel 
Schalldruck beim Empfänger entsteht. Unter der der Vereinfachung, dass 
q; = 0 ist, ergibt sich h zu: 


ik, = 3 (2.27) 


Die vibroakustische Ubertragungsfunktion kann auf einfachem Wege zum 
Beispiel mit der Impulshammermethode ermittelt werden [ZEF* 13]. Die 
zweite wichtige Größe iy ist die Eingangsnachgiebigkeit der Karosserie. 
Sie stellt den Kehrwert der lokalen dynamischen Steifigkeit der Karosserie 
(LDS) dar und setzt die Deformation an der Einleitungsstelle ins Verhältnis 
zur anregenden Kraft. Mit der bereits getroffenen Vereinfachung q, = 0, 
kann die Eingangsnachgiebigkeit beschrieben werden als: 


X2 


hh = 
F, 


(2.28) 
Oft wird auch der Kehrwert der Eingangsnachgiebigkeit angegeben, in 
diesem Zusammenhang spricht Bargende [Bar03] von der Karosserie-Ein- 
gangssteifigkeit. Für die Beschreibung der Reihenschaltung von Lager- und 
Karosserie-Vierpol müssen beide in die Kettenform überführt werden. Für 
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den Lagervierpol kann dies durch die in Tabelle 2.4 gegeben Transformati- 
onsvorschrift erfolgen. Es ergibt sich demnach: 


ln „lager Fa — 1 ch det(C™****) Fy 
x > Xo ch 1 h x 


(2.29) 
1 let, (chek chch) | Eo 


= E 
c L 
“|i ch X2 


Für die Umformung des Karosserievierpols reicht es aus, die beiden skala- 
ren Gleichungen der Hybridform umzustellen (vgl. [Gen10]): 


p, = hiig, +k > (2.30) 
X = hid, HURE? (2.31) 
Aus 2.30 lässt sich umformen: 
1 h% 
E, = (2.32) 
h> h? 
Durch Einsetzen von 2.32 in 2.31 ergibt sich: 
pa h% h6 hŠ _ Wen. (2.33) 
fee hh =. 
Somit ergibt sich der Karosserie-Vierpol in Kettenform zu: 
K 
Pa = A Karosserie P3 = l. 1 hi P3 
= = q, hh h, det (HKarosserie) q, 
= > u (2.34) 
1 | 1 hy P3 


K 
hi2 |h% (eu) (a 


Der Gesamtvierpol des Übertragungspfades vom Aggregat über das Aggre- 
gatelager und die Karosseriestruktur ergibt sich gemäß Gleichung 2.10 zu: 


G =A A (2.35) 
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1 cn + h5ydet(C'*8*") -ch hf a delle "erden H Karosserie) 


K . 
ch, his 1+ chh —h§ _ ch det(HK%osserie) 
(2.36) 
Dies führt zur Beschreibung des Übertragungspfades durch: 


F 


A| _¢|23! mitga |En 812 


x d3 &ı Ên 


(2.37) 


Um den Teilschalldruck des Übertragungspfades bestimmen zu können, 
muss der Schwingweg xı vor dem Lager gemessen werden. Dies kann über 
Piezo-Beschleunigungssensoren und eine Integration meist einfach reali- 
siert werden. Geht man weiterhin davon aus, dass q; = 0 ist, so kann der 
Teilschalldruck errechnet werden durch: 

1 chihi 


P3 = 7% = bye (2.38) 
S ga lg 1+ ehhh 


Es reicht fiir die Berechnung des Teilschalldrucks demnach aus, vier der 
acht Vierpolparameter zu bestimmen [Sel05]. Diese sind: 


e ay Transfersteifigkeit des Aggregatelagers 
° ch: Ausgangssteifigkeit des Aggregatelagers 


° h% : Eingangsnachgiebigkeit der Karosserie 


hb: vibroakustische Übertragungsfunktion der Karosserie (Noise- 
Transfer-Function (NTF)) 


Bei der Betrachtung von Gleichung 2.38 fällt auf, dass mit steigender Trans- 
fersteifigkeit des Lagers und größer werdender vibroakuakustischer Über- 
tragungsfunktion der Karosserie, der Teilschalldruck ansteigt. Des weiteren 
steigt der Teilschalldruck stark an, wenn sich das Produkt chh dem Wert 
—1 nähert. Betrachtet man die komplexen Größen ch, und h4, in Polarko- 
ordinaten, so ergibt sich eine Polstelle in der Übertragungsfunktion, wenn 
gilt: 


! 


Ic3| (cos( pz, ) + isin(Ppx, ) |n32| (cos(P,x,) + fsin(P,x,) =—1 (2.39) 
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Hierbei sind Pek, und Prk, die Argumente der komplexen Größen ch und 


h%,. Durch Umformen von Gleichung 2.39 ergibt sich: 


Ich] A (=sinp sinQyK, + COSE COSP)K, 


(2.40) 
+ j(sinz, cospe + 054 SiN, « ) =-| 
Mit den trigonometrischen Zusammenhängen 
: ; 1 
sin(x)sin(y) = 5 (cos(x — y) — cos(x + y) (2.41) 
1 
cos(x)cos(y) = > (cos(x — y) + cos(x + y) (2.42) 
1 
sin(x)cos(y) = 5 (sin(x— y) +sin(x+y) (2.43) 
ergibt sich: 
BR ! 
Ich laf] (coslo, + mg) + Jsin(@z,+94,))+-1 24 


Damit Gleichung 2.44 eine wahre Aussage ergibt, miissen drei Bedingungen 
erfüllt sein: 

|| |h2] = 1 (2.45) 

sin(P.4 + P, ) =0> L 


£92 


+o =na,neN (2.46) 
COS(PrL, + pK) =-l > Pa +P, = (2n—1)a,nEN (2.47) 


Die erste Bedingung aus Gleichung 2.45 ist erfüllt, wenn der Betrag der Ein- 
gangssteifigkeit der Karosserie dem Betrag der Ausgangssteifigkeit des La- 
gers entspricht. Die zweite Bedingung aus Gleichung 2.46 ist erfüllt, wenn 
die beiden Größen ein ganzzahliges Vielfaches von 180° als Phasendifferenz 
aufweisen. Die letzte der drei Bedingungen ist erfüllt, wenn die Phasendif- 
ferenz ein ungeradzahliges Vielfaches von 180° beträgt. 

Somit kann zusammenfassend geschlossen werden, dass Gleichung 2.44 er- 
füllt ist, wenn der Kehrwert des Betrags der Eingangsnachgiebigkeit der 
Karosserie dem Betrag der Ausgangssteifigkeit des Lagers entspricht und 
zusätzlich die Phasendifferenz zwischen diesen Größen ein ungeradzahliges 
Vielfaches von 180° beträgt. Ist dies der Fall, könnte man von einer Impe- 
danzanpassung des Lagerelements auf die Karosserie sprechen. Üblicher- 
weise wird dieser Begriff in der Hochfrequenztechnik verwendet, um die 
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Anpassung der Last an davor geschalteten Elemente zu beschreiben. Ziel 
ist es hierbei, eine möglichst verlustfreie Übertragung an den Empfänger- 
ort sicherzustellen [LLLO4]. In dem hier dargelegten Fall wäre es darüber 
hinausgehend eine geschickte Zusammenlegung von Resonanzen des Lage- 
relements und der Karosseriestruktur, welche auf der Empfängerseite (dem 
Schalldruck p3) zu einer Überhöhung ins Unendliche führt. Dieses Verhal- 
ten wirkt jedoch genau dem gewünschten Effekt nach Körperschallisolation 
durch die Motorlager entgegen. Zudem sind die entsprechenden Phasenla- 
gen nicht beliebig einstellbar bzw. insbesondere diejenige der LDS durch 
Fertigungstoleranzen starken Streuungen unterworfen. 

Abbildung 2.5 stellt beispielhaft den Betrag der lokalen dynamischen Stei- 
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Abbildung 2.5: Eingangssteifigkeit und vibroakustische Transfereigenschaften 
einer Karosseriestruktur 


figkeit und der vibroakustischen Übertragungsfunktion an einem Anbin- 
dungspunkt eines Aggregatelagers dar. Hierin ist zu erkennen, dass mit zu- 
nehmender Frequenz die lokale dynamische Steifigkeit zunimmt und die 
vibroakustische Übertragungsfunktion abnimmt. Dies ist mit lokalen Reso- 
nanzen der Karosseriestruktur begründbar, welche dafür sorgen, dass sich 
mit zunehmender Frequenz einzelne Strukturteile abkoppeln und somit we- 
niger stark an der Übertragung von Körperschall beteiligt sind. Für den 
Transferpfad vom Aggregat bis zum Insassenohr bedeutet dies, dass der 
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übertragene Schall mit zunehmender Frequenz abnimmt. Im Gesamtfahr- 
zeug beschränken sich die körperschalldominierten akustischen Phänomene 
daher zumeist auf den Frequenzbereich bis etwa 2 kHz [PBBF10]. Bei hö- 
heren Frequenzen sind in der Regel Luftschallabstrahlungen der Antriebs- 
einheit der dominantere Geräuschpfad. In der vorliegenden Arbeit soll daher 
der Frequenzbereich bis 2 kHz als der akustisch relevante betrachtet werden. 


2.2 Charakteristische Eigenschaften von Elastomeren 


Zur Analyse der Transfereigenschaften von Elastomerlagern spielt das Ver- 
ständnis der Werkstoffmechanik von Elastomeren eine herausragende Rolle. 
In den nachfolgenden Abschnitten soll daher ein Überblick über die grund- 
legenden Eigenschaften dieser Werkstoffe gegeben werden. Dies geschieht 
besonders im Hinblick auf diejenigen Eigenschaften, welche die Elastomere 
von metallischen Konstruktionswerkstoffen unterscheiden. 


2.2.1 dynamisches Verhalten von Elastomeren 


Das dynamische Verhalten von Elastomeren ist gänzlich unterschiedlich 
zum Verhalten von metallischen Werkstoffen. Die Schub- und Elastizitäts- 
moduli sind deutlich geringer, die Bruch- bzw. Reißdehnungen sind wesent- 
lich höher und können mehrere hundert Prozent betragen. Dieses Verhalten 
geht vor allem auf den molekularen Aufbau der Elastomere zurück. Sie 
bestehen aus miteinander weitmaschig vernetzten Polymerketten, welche 
durch mehr oder weniger chemische Verbindungsstellen miteinander ver- 
bunden sind. 

Die makromolekulare Struktur der Polymerketten und die hohe Zahl von be- 
weglichen, drehbaren Kohlenstoff-Kohlenstoff Bindungen sorgt dafür, dass 
sich eine große inter- und intramolekulare Beweglichkeit einstellt. Durch 
diese hohe Beweglichkeit sind die Polymerketten nicht geradlinig mitein- 
ander verknüpft sondern mit sich selbst und anderen Polymerketten ver- 
knäult. Daher liegt ohne äußere Belastung der Werkstoff in einem Zustand 
hoher Entropie vor [Rötl3]. Eine geringe Krafteinleitung führt dabei zu 
einer mechanischen Deformation und einer Ausrichtung der Ketten in die 
Kraftrichtung. Durch dieses Verhalten sind große, reversible Dehnungen 
des Werkstoffs möglich [SSK11]. Da sich durch diese Ausrichtung der Ab- 
stand benachbarter Atome innerhalb der Molekülketten nicht ändert, bleibt 
die potentielle Energie des Systems auch unverändert. Die Ausrichtung und 
das Entknäulen der Ketten führt jedoch dazu, dass die Entropie im System 
bei einer äußeren Krafteinwirkung abnimmt [SH13]. Das System wird nach 
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dieser Deformation ohne das Einwirken von äußeren Kräften nun wieder 
bestrebt sein, einen Zustand maximaler Entropie anzunehmen, was dazu 
führt, dass sich die Ketten wieder in der zuvor bestehenden Knäulstruktur 
anordnen. Es kommt daher zu einer Elastizität aufgrund der Entropie im 
Werkstoff. In diesem Zusammenhang spricht Treloar auch von der soge- 
nannten Entropie- bzw. Gummiielastizität [Tre46, Tre46, Tre75]. 


Neben der Entropieelastizität ist für das Verhalten von Elastomerwerkstof- 
fen vor allem die hohe Materialdämpfung charakteristisch. Untersuchungen 
hierzu gehen zurück auf Gehman [GWS41] und Treloar [Tre42]. Wenn ein 
Bauteil aus einem Elastomerwerkstoff mit einer sinusförmigen Kraft beauf- 
schlagt wird, so bedingt die viskoelastische Eigenschaft dieser Werkstoffe, 
dass die Deformation dem Kraftverlauf nacheilt. Dies sorgt für eine Phasen- 
differenz zwischen der Spannung und der Dehnung in Elastomeren. Dieser 
Zusammenhang ist in Abbildung 2.6 dargestellt. Die dargestellte Phasendif- 
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(© Dehnung 


Abbildung 2.6: Zeitlicher Verlauf der Spannung und Dehnung bei einer sinusförmi- 
gen Belastung von Elastomeren 


ferenz zwischen Spannung und Dehnung resultiert aus der Energiedissipa- 
tion, welche sich bei jedem Schwingspiel im Werkstoff einstellt. Die Ener- 
giemenge, die pro Schwingspiel in Wärme umgewandelt wird, entspricht 
dabei der Hysteresefläche der über der Dehnung aufgetragenen Spannung 
[GL13]. 

Diese viskosen Dämpfungseigenschaften des Werkstoffs sorgen für das 
viskoelastische Materialverhalten. Daher ist das dynamische Modul von 
Elastomeren eine komplexe Größe E*, deren Realteil E’ denjenigen Span- 
nungsanteil darstellt, welcher mit der Dehnung in Phase ist, wohingegen der 
Imaginärteil E” denjenigen Spannungsanteil widerspiegelt, welcher einen 
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Phasenversatz von 90° zur Dehnung aufweist [Geh81]. Der E-Modul von 
viskoelastischen Materialien kann demzufolge dargestellt werden durch: 


E*=E'+ je" (2.48) 


E’ wird hierbei als Speichermodul, Æ” als Verlustmodul bezeichnet. Der 
Verlustmodul ist dabei derjenige, welcher für die Energiedissipation im 
Werkstoff verantwortlich ist und ist demnach ein Maß für die Dämpfungsei- 
genschaften des Werkstoffs. Für die Beschreibung dieser Eigenschaften hat 
sich eine weitere Darstellungsform etabliert. Die Größe 6 wird als Verlust- 
winkel bezeichnet und ist durch Gleichung 2.49 gegeben [Geh57, GL13]. 


tan(d) = — (2.49) 


Der Verlustwinkel bezieht demnach den Verlustmodul auf den Speichermo- 
dul. Der Zusammenhang zwischen E*,E’, E” und 6 ist in Abbildung 2.7 in 
der komplexen Ebene dargestellt. 
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Abbildung 2.7: Darstellung des komplexen Elastizitätsmoduls von Elastomeren 
[Geh57] 


Da reine Elastomere für die meisten technischen Anwendungen aufgrund 
der geringen Reißfestigkeit und Abriebfestigkeit nicht geeignet sind, wer- 
den chemische Zusatzstoffe und feste Füllstoffe verwendet. Die chemischen 
Zusatzstoffe haben dabei die Aufgabe, die Polymerisation gezielt zu beein- 
flussen und wirken dabei als Vernetzer oder Weichmacher. Die Füllstoffe 
haben entweder aktive Funktion zur Verbesserung der mechanischen Eigen- 
schaften oder inaktive Funktion zur verbilligenden Volumenvergrößerung 
[Gra03, Wra96]. 

Als Füllstoffe werden meist Ruß- oder Silikapartikel verwendet, welche sich 
beim Compoundieren der Rohmasse und beim Vulkanisieren des Elasto- 
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mers zu Aggregaten anordnen. Die Form und Größe dieser Aggregate hat 
wiederum einen enormen Einfluss auf das Werkstoffverhalten [SH13]. Eine 
schematische Darstellung der Struktur gefüllter Elastomere ist in Abbildung 
2.8 gegeben. 


Füllstoff-Füllstoff- 
Interaktion: 


Füllstoff-Elastomer-Interaktion: 


Aggregat Füllstoffcluster 


Primärparti- 
kel 
phys. Bindung 


Abbildung 2.8: Schematische Darstellung der Struktur gefüllter Elastomere nach 
Payne [Gra03, Pay62] 


Da zwischen den Molekülen und den Füllstoffen eine komplexe Interaktion 
stattfindet, sollen die wichtigsten Wechselwirkungen zwischen diesen nach- 
folgend näher erläutert werden. 


e Hydrodynamische Verstärkung: 

Durch die Verstärkung des Elastomers mit einem Füllstoff entsteht 
ein Volumenanteil des Werkstoffs aus nicht verformbarem Fiillstoff. 
Dies hat zur Folge, dass bei gleicher makroskopischer Verformung 
die mikroskopische Verformung der Polymerketten im Vergleich zum 
ungefüllten Zustand zunehmen muss [Gra03]. Das Resultat ist eine 
Zunahme des Speichermoduls und somit eine Zunahme der Steifig- 
keit des Werkstoffs. Dieser Effekt wurde erstmals von Guth und Gold 
[GG38] fiir Suspension von Partikeln in Fliissigkeiten untersucht und 
später von Smallwood auf verstärkte Polymernetzwerke übertragen 
[Sma44]. 
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e Füllstoff-Kautschuk-Interaktion 
Durch Nebenvalenz-Effekte zwischen Füllstoffen und dem Polymer 
und chemischen Bindungen der Molekülketten an die Füllstoffpar- 
tikel nimmt die freie Kettenlänge ab und die Vernetzungsdichte des 
Polymernetzwerks steigt an. Hierdurch wird ähnlich zur hydrodyna- 
mischen Verstärkung bei gleicher makroskopischer Deformation des 
Werkstoffs die mikroskopische Verstreckung der Molekülketten er- 
höht [Gra03, SH13]. Es kommt demzufolge ebenso zu einer Verstei- 
fung des Werkstoffs durch die Füllstoff-Kautschuk-Wechselwirkung. 


Füllstoff-Füllstoff-Interaktion 

Zusätzlich zur Interaktion mit den Polymerketten, interagieren die 
Füllstoffpartikel auch untereinander. Beim Erreichen einer kritischen 
Füllstoffkonzentration kommt es nicht nur zur Ausbildung von Füll- 
stoffaggregaten sondern zur Bildung von Fiillstoffclustern. Bei wei- 
terer Erhöhung des Füllstoffanteils kommt es zur Ausbildung eines 
Netzwerkes aus Füllstoffen. Dieses Füllstoffnetzwerk liegt parallel 
zum Netzwerk der Polymerketten vor und hat eine zusätzliche ver- 
steifende Wirkung [Pay62]. 


Die Füllstoff-Füllstoff-Interaktion hat eine starke Abhängigkeit des E-Moduls 
des Werkstoffs von der Belastungsamplitude zur Folge. Dieser Effekt ist be- 
kannt als Payne-Effekt und soll im Folgenden erläutert werden. 


Payne-Effekt 


Die Untersuchungen von Payne [Pay62, Pay63, Pay64] zur Amplitudenab- 
hängigkeit der dynamischen Moduli von mit Ruß gefüllten Elastomerwerk- 
stoffen zeigen eine starke Abnahme des Speichermoduls mit wachsender 
Dehnungsamplitude. Diese Nichtlinearität des Werkstoffverhaltens wird in 
der Literatur als Payne-Effekt bezeichnet. Der Wert des Speichermoduls nä- 
hert sich dabei einem Grenzwert G/,. Für sehr kleine Dehnungsamplituden 
nähert sich der Speichermodul ebenfalls einer waagerechten Asymptote Go. 
Der Verlustmodul G” hat im Bereich der größten Änderung von G’ sein 
Maximum G/,.. Die Verläufe des Speicher- und Verlustmoduls über der 
(Schub-)Dehnungsamplitude y sind in Abbildung 2.9 dargestellt. Payne 
fiihrt diese Abnahme des Speichermoduls auf das Zusammenbrechen des 
Füllstoffnetzwerkes zurück. Die verstärkende Wirkung der Füllstoff-Füllstoff- 
Interaktion nimmt dadurch stetig ab. Der Grenzwert G/, ist dadurch lediglich 
durch die Anteile des Polymernetzwerks, der hydrodynamischen Verstär- 
kung und der Füllstoff-Kautschuk-Wechselwirkungen bestimmt. 
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Abbildung 2.9: Darstellung der Abhängigkeit des dynamischen Schubmoduls von 
der Dehnungsamplitude nach Payne [Pay62, Pay63, Pay64] 


Die verstärkende Wirkung der Füllstoff-Füllstoff-Interaktion nimmt mit zu- 
nehmendem Volumenanteil der Füllstoffe zu. Dies hat zur Folge, dass die 
Differenz zwischen G) und G/, mit zunehmendem Füllstoffanteil größer 
wird und demnach auch die Abnahme des Speichermoduls durch den Payne- 
Effekt mit zunehmendem Füllstoffanteil ansteigt. 

Payne [Pay60] hat diesen Effekt für Naturkautschuk-Vulkanisate mit un- 
terschiedlichen Volumenanteilen von Ruß als Füllstoff untersucht. Hierzu 
wurden zylindrische Proben mit einer Frequenz von 0,5 Hz bei verschiede- 
nen Dehnungsamplituden zyklisch belastet. Die gemessenen Werte für die 
E-Moduli sind in Abbildung 2.10 über der Dehnungsamplitude dargestellt. 
Der Payne-Effekt regeneriert sich innerhalb weniger Sekunden nach der 
dynamischen Belastung, was dazu führt, dass diesem nicht durch eine Vor- 
konditionierung des Werkstoffs begegnet werden kann [RKS*11, HLO9]. 


Mullins-Effekt 


Die Interaktionen zwischen Füllstoff und Polymerketten bedingen noch eine 
weitere Nichtlinearität des Werkstoffverhaltens von gefüllten Elastomeren. 
Während einer zyklischen Belastung des Werkstoffs kommt es bei den er- 
sten Belastungszyklen zu einer zunehmenden Entfestigung des Prüfkörpers. 
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Abbildung 2.10: Darstellung des Payne-Effekts für Naturkautschuk-Vulkanisate 
mit unterschiedlichen Volumenanteilen von MAF-Ruß [Pay62, 
Pay63, Pay64] 


Diese Abnahme des dynamischen Moduls ist bekannt als Spannungserwei- 
chung oder Mullins-Effekt [MT65, Joh92]. Der Mullins-Effekt wird dadurch 
erklärt, dass einzelne chemische Verbindungen in den Polymerketten durch 
die erste maximale Dehnung aufbrechen und sekundäre Bindungen zwi- 
schen dem Polymernetzwerk und den Füllstoffen aufbrechen [Sto99]. Die- 
ses Aufbrechen von chemischen Bindungen hat eine bleibende Deformation 
des Werkstoffs zur Folge, welche auch nach einer vollständigen Entlastung 
noch erhalten bleibt. 

In Abbildung 2.11 ist das Spannungs-Dehnungs-Diagramm für eine Be- 
lastung eines gefüllten Elastomers durch drei aufeinanderfolgende Bela- 
stungszyklen dargestellt. Es kann hierin erkannt werden, dass nach der er- 
sten Belastung eine residuale Dehnung erhalten bleibt, welche sich nach 
dem zweiten Zyklus noch erhöht. Ebenso zeigt sich die Entfestigung des 
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Werkstoffs, da beim zweiten und dritten Belastungszyklus die maximale 
Spannung geringer ausfällt als bei dem jeweiligen vorhergehenden. Nach 


> 


Spannung o 


ne E 


Abbildung 2.11: Darstellung des Mullins-Effekts im Spannungs-Dehnungs- 
Diagramm nach [RS12] 


einer vom Polymernetzwerk und dem Füllgrad des Werkstoffs abhängigen 
Anzahl von Belastungszyklen kommt es zu einem stationären Zustand, von 
dem aus keine weitere Entfestigung mehr stattfindet. In diesem Zustand 
wird der Werkstoff als vorkonditioniert betrachtet [Med78]. Eine weitere 
Erweichung tritt erst dann ein, wenn die Dehnungsamplitude über den be- 
reits zuvor erreichten Maximalwert hinausgeht. Diani et al. [DFG09] ha- 
ben Untersuchungen zur Reversibilität des Mullins-Effekts angestellt und 
eine vollständige Erholung des Moduls nur nach längerer Temperaturein- 
wirkung festgestellt, wohingegen die residuale Dehnung nach längerer Zeit 
auch bei Raumtemperatur wieder vollständig verschwindet. Es kann festge- 
halten werden, dass der Mullins-Effekt durch die maximale Deformation des 
Werkstoffs bestimmt wird, welcher dieser während seiner Belastungshisto- 
rie erfahren hat [Lio06]. Für den Bauteilversuch an Elastomerbauteilen für 
Fahrzeuge muss aufgrund des Mullins-Effektes daher beachtet werden, dass 
eine entsprechende Konditionierung erfolgen muss. Im konkreten Fall von 
Aggregatelagerbauteilen bedeutet dies, dass die Bauteile zunächst mit der 
maximalen im Betrieb auftretenden Dehnung beaufschlagt werden sollten, 
um einen für den Fahrbetrieb relevanten Konditionierungszustand einzustel- 
len. Dies wird in der Regel erreicht, indem das Bauteil quasistatisch bis in 
die Endanschläge deformiert wird [Tre15]. 
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Einfluss der Temperatur auf das Werkstoffverhalten 


Die Besonderheiten in der Mikrostruktur von Elastomeren haben als Folge, 
dass die Eigenschaften von Elastomeren stark von der Temperatur des Werk- 
stoffs abhängen. Charakteristisch für das Verhalten von Elastomeren ist eine 
starke Zunahme des Elastizitätsmoduls bei abnehmender Temperatur. Diese 
Zunahme geschieht in einem kleinen Temperaturbereich und hat eine starke 
Abnahme der Reißdehnung zur Folge, was zu einem glasartigen Verhalten 
des Werkstoffs führt. Aus diesem Grund wird die Temperatur, bei der die- 
se Zunahme des Moduls erfolgt, als Glasübergangstemperatur bezeichnet. 
Ursache für dieses Verhalten ist eine Abnahme der Beweglichkeit einzelner 
Kettensegmente [SD86]. Die Bewegungsfreiheitsgrade innerhalb der Poly- 
merketten aufgrund von Rotationen innerhalb der Bindungen zwischen den 
Monomereinheiten werden dabei eingeschränkt [BS92]. Die Reduktion der 
Kettenbeweglichkeit führt dazu, dass die Entropieelastizität des Werkstoffs 
nicht mehr gegeben ist und das Verhalten nahezu ausschließlich durch die 
Energieelastizitat bestimmt wird [SD86]. Die Abhängigkeit des komple- 
xen Elastizitätsmoduls von der Temperatur ist schematisch in Abbildung 
2.12a dargestellt. Abbildung 2.12a zeigt, dass im Hauptanwendungsbereich 


(b) 


Temperatur in log K Frequenz in log Hz 


Abbildung 2.12: Schematische Darstellung der Abhängigkeit des komplexen Mo- 
duls von der Temperatur (a) und der Frequenz (b) nach [Gro13] 


von Elastomeren, oberhalb der Glasübergangstemperatur, die dynamischen 
Moduli stetig abnehmen. Daraus lässt sich folgern, dass mit zunehmender 
Temperatur die Werkstoffe immer weiter erweichen. 


Die Untersuchungen von Williams et al. [WLF55] zeigen, dass ein weiterer 
Effekt das dynamische Verhalten von Elastomeren beeinflusst. Eine Erhö- 
hung der Belastungsgeschwindigkeit bzw. der Belastungsfrequenz sorgt für 
eine zunehmende Versteifung des Werkstoffs. Auch hierbei kommt es bei 
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hohen Frequenzen zu einem starken Anstieg der dynamischen Moduli. Die- 
ses Verhalten ist in Abbildung 2.12b dargestellt. Es kann festgehalten wer- 
den, dass das Verhalten von Elastomeren durch die Bauteiltemperatur und 
die Belastungsfrequenz in gleichem Maße beeinflusst wird. Eine Erhöhung 
der Belastungsfrequenz entspricht dabei einer Verringerung der Temperatur 
des Werkstoffs und umgekehrt. Daher wird dieses Verhalten auch als Zeit- 
Temperatur-Verschiebung bezeichnet. Die nach Williams, Landel und Ferry 
[WLF55] benannte WLF-Gleichung stellt diesen Zusammenhang mathema- 


tisch dar. f 
I 2.50 
f ) OG+(T-n) ha 


Hierbei steht ar für das Verhältnis der Belastungsfrequenz fo bei Referenz- 
temperatur Tọ zu einer äquivalenten Belastungsfrequenz f bei einer beliebi- 
gen Temperatur T. Die Werte C1 und C2 sind materialabhängige Konstanten, 
wobei Cı dimensionslos ist und die Konstante C> die Einheit Kelvin besitzt. 
Die WLF-Gleichung ermöglicht es, das Verhalten von Elastomeren bei einer 
Temperaturänderung zu beschreiben oder das Verhalten bei Frequenzen zu 
bestimmen, die mit den zur Verfügung stehenden Messgeräten nicht erreicht 
werden können. 


Ci (T _ To) 


log(ar) = log( 


2.2.2 Modellierung des Werkstoffverhaltens von Elastomeren 
im akustisch relevanten Frequenzbereich 


Die Beschreibung der viskoelastischen Materialeigenschaften von Elasto- 
meren in Simulation mit finiten Elementen kann nur unter bestimmten Vor- 
aussetzungen erfolgen. Diese grundlegenden Annahmen und die Material- 
modellierung mittels der für diese Arbeit verwendeten FE-Software DS Si- 
mulia ABAQUS soll in diesem Abschnitt erläutert werden. 

Die in Abschnitt 2.2 erläuterten Eigenschaften, die charakteristisch für das 
Verhalten von in Aggregatelagerbauteilen verwendeten Elastomeren sind, 
können im akustisch relevanten Frequenzbereich am einfachsten über ein 
lineares, viskoelastisches Materialmodell abgebildet werden. Dieses Mate- 
rialmodell setzt voraus, dass die Deformation der einzelnen Elemente klein 
ist und dass die Schubanteile (deviatorische Anteile) im multiaxialen Span- 
nungszustand unabhängig von den volumetrischen Spannungsanteilen sind 
[HKS14, Baal0]. 

Geht man analog zur Darstellung in Abschnitt 2.2.1 von einer sinusförmi- 
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gen Schubdeformation y(t) mit der Amplitude % des Werkstoffs mit der in 
Gleichung 2.51 dargestellten Form aus, 


y(t) = ne!” (2.51) 


so ergibt sich die resultierende Schubspannung zu: 
a(t) = (G(@) + jG" (@) woe (2.52) 


Hierin sind T(t) die resultierende Schubspannung, G’ und G” die Schub- 
speicher- und Verlustmoduli. 


Um das viskoelastische Materialverhalten darzustellen, werden rheologi- 
sche, zeitabhängige Modelle eingesetzt. Zur Abbildung dieses Verhaltens 
sind bereits eine große Anzahl von Modellen entwickelt worden, welche ei- 
ne beliebig komplexe Kombination aus Feder, Dämpfer und Reibelementen 
verwenden. In ABAQUS wird die Viskoelastizität sowohl im Zeit- als auch 
im Frequenzbereich über eine Parallelschaltung von n Maxwell-Elementen 
und einem Federelement modelliert. Dieses Modell wird auch als generali- 
siertes Maxwell-Modell bezeichnet [HKS14, GL13]. Wird das generalisier- 
te Maxwell-Modell einem Dehnungssprung unterzogen, so relaxiert das ge- 
samte Modell mit der Zeit gegen einen konstanten Spannungswert [Wril3]. 
Jedes einzelne der parallelgeschalteten Maxwell-Elemente hat dabei eine ei- 
gene Relaxationszeit und eine Relaxationsstärke. Durch eine hohe Anzahl 
von Maxwell-Elementen kann das Materialverhalten approximiert werden. 
Jedes Maxwell-Element repräsentiert dabei eine Stützstelle im Relaxations- 
spektrum des Materials [KT11]. Nach [HKS14] kann das Relaxationsver- 
halten des generalisierten Maxwell-Modells beschrieben werden durch: 
N 
G(t) = Go + È Ge (2.53) 
zu 


i= 


Hierbei ist G» das Langzeit-Schubmodul, G; das Relaxationsmodul des 
i-ten Maxwell-Elements, 7; der Relaxationskoeffizient des i-ten Maxwell- 
Elements und n die Anzahl der Maxwell-Elemente im generalisierten Maxwell- 
Modell. 

In ABAQUS wird der normierte Relaxationsmodul eingefiihrt, der definiert 


ist als: 


Gt). N 
g(t) = © mit Go = Ge +26 (2.54) 
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mit Go dem Anfangsschubmodul zum Zeitpunkt 7 = 0. Somit ergibt sich für 
das normierte Gesamtrelaxationsmodul [Ber05]: 


=1- 3 Gi(1-e”/®) (2.55) 


Aus dynamisch-mechanischen Analysen ist eine Ermittlung der Prony- 
Terme nicht direkt möglich, eine Ermittlung von dynamischem Speicher- 
Schubmodul und Verlust-Schubmodul ist in der Regel einfacher. Über die 
Fourier-Transformation von Gleichung 2.53 kann ein Zusammenhang zwi- 
schen dem frequenzabhängigen komplexen Schubmodul und den Prony- 
Parametern gewonnen werden [HKS 14, BMD07]: 


N Gio? 
G'(@ =Go{l—) GJ+G 2.56 
ol È + rear I+ Pa? (2.56) 
N Gito 
G"(@) = Go Y ——— 257 
(o) Li eel (2.57) 


Die Parameter der Prony-Reihe können durch M experimentell ermittelte 
Wertepaare für Speicher- und Verlustmodul (G4, ;p und Glip) über eine Mi- 
nimierung der Fehlerfunktion X bestimmt werden [HKS14, FCR11]: 


1 
G 


ro [(G' (0i) — Gexp(@i))7 + (G" (0) — Gexp(@i))?] 2.58) 


imes 


gN 


Dieses Vorgehen kann analog für die Ermittlung der Prony-Parameter des 
komplexen Kompressions- bzw. Elastizitätsmoduls erfolgen. Unter der An- 
nahme eines isotropen Materialverhaltens kann eine Ermittlung des Kom- 
pressions- und/oder des Elastizitätsmoduls auch auf Basis der komplexen 
Querkontraktionszahl erfolgen [KT11]. 

Der Zusammenhang zwischen dem Schubmodul und den weiteren Materi- 
alparametern ergibt sich gemäß [BMP02] zu: 


1 _ 3(1—2v) 3KE 


G= 7040)" aaa R-E 


(2.59) 


mit: G: Schubmodul; E: Elastizitätsmodul; K: Kompressionsmodul und v: 
Querkontraktionszahl 
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Elastomere verhalten sich bei kleinen Dehnungsamplituden und niedriger 
Belastungsfrequenz bei dynamisch-zyklischer Beanspruchung nahezu in- 
kompressibel [MS04]. Dieses Verhalten bedingt, dass die Querkontrakti- 
onszahl den Wert 0,5 annimmt und das Kompressionsmodul sehr groß wird. 
Da ein Wert für v von 0,5 für numerische Instabilität in Simulationsmodel- 
len führt, wird üblicherweise ein Wert von 0,495 — 0,498 angenommen 
[DM99, Bro06]. 

Betrachtet man die dynamischen Moduli als komplexe Größen, so ist auch 
die Querkontraktionszahl komplex und hat einen geringen Imaginäranteil, 
welcher aber nach Pritz [Pri00] klein ist. 

Die Untersuchungen von Pandini et al. [PP08], Caracciolo et al.[CG96] 
und Pritz [Pri00, Pri98] zeigen, dass die Querkontraktionszahl über der Fre- 
quenz und der Temperatur nicht konstant ist. Pandini zeigt in seinen Be- 
obachtungen, dass der Zusammenhang, welcher tiber die WLF-Gleichung 
gegeben ist (vgl. Abschnitt 2.2.1) in gleichem Maße auch für das Verhal- 
ten der Querkontraktionszahl gegeben ist. Dies bedeutet, dass auch für die 
Querkontraktionszahl eine Äquivalenz zwischen einer Erhöhung der Bela- 
stungsfrequenz und einer Absenkung der Temperatur des Werkstoffs exi- 
stiert. Abbildung 2.13 zeigt den Verlauf des Betrages der Querkontrakti- 
onszahl eines Elastomer-Werkstoffes über der Belastungsfrequenz. Es ist 
zu erkennen, dass der Betrag der komplexen Querkontraktionszahl des Ela- 
stomers über der Frequenz abnimmt und bis zu einer Frequenz von circa 
100Hz ungefähr dem Wert 0,5 entspricht. Dieses Verhalten ist, wie auch 
der Verlauf des Speicher- und Verlustmoduls, abhängig von der jeweiligen 
Kautschukmischung. 

Abbildung 2.14 zeigt die komplexe Querkontraktionszahl eines Styrol- 
Butadien-Elastomers für verschiedene Frequenzen bei Raumtemperatur in 
der komplexen Ebene. 

Es kann erkannt werden, dass der Betrag der Querkontraktionszahl bei nie- 
derfrequenter Anregung sehr nahe am Wert 0,5 liegt und der Imaginärteil 
nahe Null ist. Mit steigender Frequenz nimmt der Betrag der Querkontrak- 
tionszahl ab und der Imaginärteil nimmt negative Werte an, welche aber 
für den Frequenzbereich bis 100 kHz zwischen —0,01 und 0 liegen. Für den 
Betrag der Querkontraktionszahl ist der Imaginärteil daher vernachlässigbar 
gering. 

Um das Verhalten von Elastomerwerkstoffen für die Analyse von Systemen 
im akustisch relevanten Frequenzbereich in einer Simulationsumgebung ab- 
bilden zu können, muss zwingend das Verhalten der Querkontraktionszahl 
berücksichtigt werden. Dies kann erreicht werden, wenn sowohl das Schub- 
modul als auch das Kompressionsmodul nach dem in Gleichungen 2.56, 
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Abbildung 2.13: Verlauf des Betrags der Querkontraktionszahl über der Bela- 
stungsfrequenz des Werkstoffs für einen Acryl-Elsatomer bei 
Raumtemperatur nach [Pri00] 


Im(v) 


-0,05 


Abbildung 2.14: Darstellung der Querkontraktionszahl eines Elastomers auf Styrol- 
Butadien-Basis in der komplexen Ebene fiir verschiedene Fre- 
quenzen bei Raumtemperatur nach [Pri00] 
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2.57 und 2.58 beschriebenen Verfahren an das reale Werkstoffverhalten an- 
gepasst werden. 


2.3 Systembeschreibung von Aggregatelagerbauteilen 


Zum Verständnis des Aufbaus und der Funktionsweise von Aggregatelager- 
bauteilen soll in den folgenden Abschnitten ein Überblick über die für die 
Lagerung des Triebstrangs im Fahrzeug üblicherweise verwendeten Bautei- 
le gegeben werden. Zunächst soll jedoch eine Einführung in die Aggrega- 
telagerung und die Aufgaben, welche diese im Gesamtfahrzeug überneh- 
men, eingegangen werden. Im Anschluss an die Systembeschreibung der 
verschiedenen Bauteile wird auf die Herstellung von Elastomerbauteilen zur 
Aggregatelagerung eingegangen, da die Prozessführung bei der Herstellung 
der Elastomerkomponenten einen großen Einfluss auf das Verhalten des Ge- 
samtbauteils vor allem im akustischen Frequenzbereich hat. 


2.3.1 Aufgaben der Aggregatelagerung 


Als Aggregatelager werden diejenigen Bauteile bezeichnet, welche den 
Triebstrang mit der Karosseriestruktur des Fahrzeuges verbinden [ZEF* 13]. 
Der Triebstrang umfasst dabei den Motor-Getriebe-Verbund sowie je nach 
Antriebstrangkonzept Differentiale und Gelenkwellen. Diese Bauteile müs- 
sen mit dem Fahrzeug verbunden und ihre Gewichtskräfte müssen abge- 
stützt werden. 

Der Motor-Getriebe-Verbund kann entweder direkt an der Karosserie oder 
über einen Hilfsrahmen mit der Fahrzeugstruktur verbunden werden. Dif- 
ferentiale werden häufig auf einem mit einer weiteren Isolationsebene zum 
Fahrzeug gelagerten Achsträger, dem sogenannten Fahrschemel, gelagert 
[SHB13]. Abbildung 2.15 zeigt beispielhaft die schematische Darstellung 
der Aggregatelagerung eines Antriebstranges für ein Fahrzeug mit Standar- 
dantrieb. 

Für dieses Triebstrangkonzept besteht die Aggregatelagerung aus zwei Mo- 
torlagern (1), einem Getriebelager (2), einem Gelenkwellenmittellager (3) 
und den vier Hinterachsdifferential- bzw. Hinterachsgetriebelagern (4). In 
Fahrzeugen mit starker Motorisierung und einem entsprechend hohen Mo- 
tormoment, welches sich hauptsächlich über die beiden Motorlager abstützt, 
kommt in manchen Fällen noch eine Momentenstütze zum Einsatz, welche 
die Drehbewegung des Aggregats begrenzt [HEG13]. Die Hauptaufgaben 
der Aggregatelagerung bestehen dabei aus [ZEF*13, Hem02]: 
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OF Motorlager ®: Getriebelager ©: Gelenkwellenmittellager 
©) : Hinterachsgetriebe-/Differentiallager 


Abbildung 2.15: Schematische Darstellung der Aggregatelagerung eines Fahrzeugs 
mit Standardantrieb [ZEF* 13] 


e mechanische Befestigung des Aggregats: 


— Abstützen von Gewichtskräften; 
— Ausgleich von Toleranzen im Rohbau; 


— Begrenzung der Auslenkung des Triebstrangs aus der Ruhelage 


e Abstützen von dynamischen Belastungen: 


— Abstützen des Motor-/Getriebemoments; 


— Abstützen von Trägheitskräften in Folge von Reifen-Fahrbahn- 
Erregung; 


Bereitstellung von Dämpfung für die Moden des Aggregates; 


Abstützen von Trägheitskräften in Folge von Fahrzeuglängs- 
und Querdynamik 


Abstützen von Crashlasten 


Isolation von triebstrangerregten Schwingungen: 
— Isolation der oszillierenden Massenkräfte; 
— Isolation von Gaskräften des Verbrennungsprozesses; 
— Isolation von hochfrequenten Schwingungen aus Nebenantrie- 


ben und Verzahnungen/Kettentrieben 
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Die Anforderungen, die sich aufgrund der genannten Aufgaben der Aggre- 
gatelagerung an die Bauteile ergeben, sind teilweise im Widerspruch zuein- 
ander. So sind für das Abstützen von hohen Kräften steife Lagerelemente 
von Vorteil, was jedoch mit einer Verschlechterung des Isolationsverhaltens 
einhergeht. Um die Vielzahl an Anforderungen zu erfüllen, haben sich für 
die Aggregatelagerung einige Bauformen von Lagerelementen durchsetzen 
können. Diese sollen im Folgenden näher erläutert werden. 


2.3.2 Buchsenlager 


Die häufigste und einfachste Bauform von Aggregatelagern ist die Buch- 
se. Buchsen zeichnen sich vor allem durch eine sehr kleine Bauform und 
die geringen Bauteilkosten aus. In den nachfolgenden Abschnitten wird de- 
tailliert auf den Aufbau und die bekannten Eigenschaften dieser Bauteile 
eingegangen. 


Aufbau und Materialien 


Im Allgemeinen besteht eine Buchse aus einem Außenring, einem axial 
überstehenden Kern und einem Elastomeranteil, der diese beiden Bauteile 
miteinander verbindet. In Abbildung 2.16 ist ein Schnitt durch ein solches 
Buchsenlager dargestellt. Die Anbindung an die umgebenden Strukturen 


(1):Außenring; (2): Kern; (3): Durchgangsbohrung; (4): Elastomeranteil 


Abbildung 2.16: Darstellung eines Schnitts durch ein Buchsenlager 
erfolgt für den Außenring üblicherweise durch eine Übermaßpassung, für 
den Kern mit Hilfe einer Schraubenverbindung durch die enthaltene Durch- 


gangsbohrung [Tre15]. 
Der Außenring und der Kern werden in der Regel aus Aluminium gefertigt. 
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Es gibt jedoch auch Sonderformen, bei der aus Gewichtsgründen und aus 
Gründen der thermischen Isolation Außenring und/oder Kern aus glasfa- 
serverstärkten Kunststoffen hergestellt werden. Dies ist vor allem dann der 
Fall, wenn die umgebenden Strukturen im Fahrbetrieb hohe Temperaturen 
aufweisen, welche die Eigenschaften des Elastomers beeinflussen oder es zu 
einer thermo-oxidativen Schädigung des Materials kommt. In diesen Fällen 
werden auch synthetische Elastomere auf Basis von Silikonkautschuk oder 
Ethylen-Propylen-Dien-Kautschuk (EPDM) verwendet. Diese weisen im 
Vergleich zu Naturkautschuk eine höhere Dauereinsatztemperatur auf. Für 
die meisten Anwendungen für Aggregatelagerbauteile und Fahrwerkskom- 
ponenten kommen jedoch aufgrund der guten mechanischen Eigenschaften 
und der geringen dynamischen Verhärtung des Werkstoffes Naturkautschuk- 
Elastomere (NR) zum Einsatz [Tre15, Hem02]. 

Die Übertragungseigenschaften von Buchsenlagern werden nahezu aus- 
schließlich durch die Eigenschaften dieses Elastomeranteils definiert. Die 
Einflussfaktoren auf dieses Verhalten sollen daher im nächsten Abschnitt 
erläutert werden. 


Eigenformen und Einflussfaktoren 


Die Hauptauslegungsmerkmale von Buchsenlagern sind die quasistatischen 
und dynamischen Transfersteifigkeiten in die drei translatorischen Raum- 
richtungen und je nach Anwendungsfall auch die rotatorischen Steifigkei- 
ten. Diese sind vor allem bei einer kardanischen Belastung des Bauteils, 
welche für die Nebenfederraten bei einem Einsatz in Fahrwerkskomponen- 
ten verantwortlich sind, von Relevanz [Blol0]. Werden Buchsen als Ag- 
gregatelager eingesetzt, sind die Verdrehwinkel der Bauteile aufgrund der 
großen Abstände zwischen den Lagerelementen klein und die kardanische 
Belastung wird meist vernachlässigt [SR11]. 

Die Abstimmung der Steifigkeitsverhältnisse der einzelnen Raumrichtungen 
erfolgt über die geometrische Formgebung der Gummikontur. In der Praxis 
haben sich für Aggregatelagerbuchsen sogenannte Viersteg- und Zweisteg- 
Buchsenlager durchgesetzt [PBA13, FDRS 16, GKVGNO07]. Abbildung 2.17 
zeigt diese beiden typischen Bauformen. 

Die Ausführung der Buchse über eine Zweisteg-Geometrie hat zur Folge, 
dass das Lagerelement in derjenigen Raumrichtung, in der die Stege ange- 
ordnet sind (x-Richtung in Abbildung 2.17), eine deutlich höhere Steifigkeit 
aufweist, als in die dazu orthogonale Richtung (z-Richtung in Abbildung 
2.17). Bei der Ausführung mit vier Stegen, ist die Steifigkeit der beiden ra- 
dialen Richtungen ähnlich. Die in Abbildung 2.17 b) abgebildete Buchse 
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Abbildung 2.17: Darstellung eines Zweisteg-Buchsenlagers (a) und eines vier-Steg- 
Buchsenlagers (b) 


hat aufgrund der Winkelstellung der Stege verschiedene Steifigkeiten in die 
x- und z-Richtung. Die Differenz zwischen diesen ist jedoch geringer aus- 
geprägt als bei der abgebildeten Zweisteg-Buchse. 

Die Abstimmung der Steifigkeitsverhältnisse wird bei Buchsenlager im We- 
sentlichen von der Anzahl der Gummistege, dem Winkel der Gummistege zu 
den einzelnen Raumrichtungen sowie den Querschnitten der Gummisehnen 
bestimmt. Göbel [G6b69] hat Grundlagen zu diesen Auslegungen erarbeitet 
und zeigt die Einflüsse verschiedener Designparameter der Gummisehnen 
auf das quasistatische Verhalten. 


(Quasi-)Statische und niederfrequente Eigenschaften 


Das quasistatische Verhalten der Bauteile zeichnet sich durch einen großen 
Bereich linearer Steifigkeit aus, welcher durch eine progressive Zunah- 
me der Steifigkeit bei großen Auslenkungen beschränkt wird. Abbildung 
2.18 zeigt die quasistatische Kraft-Weg-Hysterese und die daraus resul- 
tierende Transfersteifigkeit des in Abbildung 2.17 dargestellten Viersteg- 
Buchsenlagers. Der Bereich linearer Steifigkeit von ca. 1200 N/mm er- 
streckt sich bei diesem Bauteil auf ca. +3 mm. Für Auslenkungen außerhalb 
dieses Bereichs nimmt die Steifigkeit stark zu, da der Kern gegen die auf 
dem Außenring vulkanisierten Anschläge gepresst wird. Aufgrund der ge- 
ringen Kompressibilität des Elastomerwerkstoffs kommt es zu dem starken 
Anstieg der Steifigkeitswerte. Das Kraft-Weg-Diagramm weist eine Hyste- 
resefläche auf, welche durch die Werkstoffdämpfung bedingt ist. Die Ver- 
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lustarbeit, welche bei einer quasistatischen Belastung auftritt, ist jedoch so 
gering, dass diese in theoretischen Betrachtungen meist vernachlässigt wird 
[PH02]. 

Aufgrund der hydrodynamischen Verstärkung von gefüllten Elastomeren 
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Abbildung 2.18: (Quasi-)statische Kraft-Weg-Hysterese und Steifigkeit eines Buch- 
senlagers 


kommt es jedoch mit einer zunehmenden Belastungsgeschwindigkeit so- 
wohl zu einer Zunahme der dynamischen Steifigkeit als auch des Verlust- 
winkels. Die dynamischen Transfersteifigkeiten fiir niedrige Frequenzen 
und verschiedene Anregungsamplituden des bereits beschriebenen 4-Steg- 
Buchsenlagers können Abbildung 2.19 entnommen werden. Hierbei kann 
festgestellt werden, dass mit einer Abnahme der Anregungsamplitude eine 
Zunahme der Steifigkeit einhergeht. Dies ist durch den in Abschnitt 2.2.1 
erläuterten Payne-Effekt begründbar. Mit einer Abnahme der Auslenkungs- 
amplitude wirken sich die verstärkenden Effekte durch die Füllstoffe stärker 
aus, da es mit einer hohen Deformation zu einer zunehmenden, reversi- 
blen Auflösung des Füllstoff-Netzwerkes kommt. Für Elongationen, welche 
kleiner als 0,1 mm sind, kommt es zu keiner weiteren Zunahme der Trans- 
fersteifigkeit mehr, was ebenfalls mit den Untersuchungen von Payne in 
Einklang steht. 


Betrachtet man das Verhalten über der Frequenz, so kann festgestellt wer- 
den, dass die Transfersteifigkeit mit zunehmender Frequenz monoton an- 
steigt und ab etwa 10Hz einen konstanten Gradienten aufweist. Dieser 
Anstieg ist durch die viskoelastischen Eigenschaften des Elastomers be- 
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Abbildung 2.19: Transfersteifigkeiten eines Buchsenlagers bei Anregung mit ver- 
schiedenen Wegamplituden tiber der Frequenz 


dingt. Durch die Zunahme der Belastungsgeschwindigkeit mit der Frequenz 
kommt es durch die viskosen Anteile des Werkstoffverhaltens zu einer Ver- 
steifung des Materials. Zusammenfassend kann festgestellt werden, dass das 
Verhalten von Buchsenlagern im tieffrequenten Bereich vor allem von der 
Elastomermischung abhängt und sich mit zunehmender Frequenz und ab- 
nehmender Belastungsamplitude eine Versteifung des Werkstoffs einstellt. 
Die Versteifung ist dabei wiederum abhängig von der Elastomermischung, 
da diese vor allem durch die Füllstoff-Interaktionen bedingt wird (vgl. Ab- 
schnitt 2.2.1). 


2.3.3 Hydraulisch dämpfende Aggregatelager 


Das zunehmende Gewicht von Motor-Getriebe-Verbunden und die steigen- 
den Anforderungen bezüglich des Schwingungskomforts in Fahrzeugen ha- 
ben in den 80er Jahren des 20. Jahrhunderts dazu geführt, dass die bis dahin 
aus einem Elastomertragkörper bestehenden Aggregatelager durch eine hy- 
draulische Dämpfung erweitert wurden [Hem02]. Der Aufbau und die Ei- 
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genschaften dieser hydraulisch dämpfenden Aggregatelager sollen in den 
nächsten Abschnitten erläutert werden. 


Aufbau und Materialien 


Die hydraulisch dämpfenden Aggregatelager dienen dazu, dem Zielkon- 
flikt zwischen der Notwendigkeit hoher Dämpfung für tieffrequenten Fahr- 
komfort und den Isolationseigenschaften von niedrigdämpfendem Elasto- 
mer durch eine Trennung der Funktion zu begegnen [Ham85]. Die last- 
tragende Gummifeder (Tragfeder) wird dabei aus einem hochelastischen, 
niedrigdämpfenden Elastomer hergestellt. Die Dämpfung bei niedrigen Fre- 
quenzen und großen Amplituden wird von einem separaten hydraulischen 
System erzeugt. Abbildung 2.20 zeigt schematisch den Aufbau eines sol- 
chen hydraulischen Motorlagers. In Abbildung 2.21 ist ein Schnitt durch 
ein solches Lager dargestellt. 

Das Lager besteht aus der Tragfeder, unter welcher sich das hydraulische 
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Abbildung 2.20: Schematische Darstellung des Aufbaus eines hydraulisch dämp- 
fenden Motorlagers [WGP17] 


System befindet. Der aggregateseitige Anbindungspunkt befindet sich da- 
bei in der Regel in der Mitte der Tragfeder, die karosserieseitige Anbin- 
dung an einem Gehäuse, welches den äußeren Teil der Tragfeder und das 
Düsen-Membran-System fasst. Das Düsen-Membran-System schließt zu- 
sammen mit der Innenseite der Tragfeder eine hydraulische Arbeitskammer 
ein. Unterhalb des Düsen-Membran-Systems befindet sich die sogenann- 
te Ausgleichskammer, welche durch eine biegeweiche Rollmembran gegen 
die Umgebung abgedichtet ist. Das gesamte hydraulische System ist mit ei- 
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Abbildung 2.21: Darstellung eines hydraulisch dämpfenden Motorlagers 


nem Dämpfungsmedium (in der Regel Glykol) gefüllt [KT85, SA85]. 


Findet zwischen Aggregat und Karosserie eine Relativbewegung statt, so fe- 
dert die Tragfeder ein bzw. aus und das Volumen der Arbeitskammer ändert 
sich. Damit dies stattfinden kann, muss sich ein Volumenstrom zwischen 
Ausgleichs- und Arbeitskammer durch das Düsen-Membran-System ein- 
stellen. Die Flüssigkeit kann dabei durch zwei Gitter strömen, zwischen de- 
nen sich eine Entkoppelmembran befindet. Diese hat einen definierten Frei- 
weg zwischen den Gittern und kann für einen Volumenausgleich von der 
Arbeitskammer in die Ausgleichskammer sorgen. Ist der Relativweg groß 
genug, so wird die Ausgleichsmembran gegen eines der beiden Gitter ge- 
drückt und es kann kein Volumenausgleich mehr über die Membran stattfin- 
den. Es wird dann in der Arbeitskammer Druck aufgebaut, was ein Blähen 
der Tragfeder zur Folge hat. Infolge dessen strömt das Dämpfungsfluid dann 
durch den ebenfalls im Düsen-Membran-System enthaltenen Dämpfungs- 
kanal. Aus dem Flüssigkeitsvolumen im Dämpfungskanal und der Feder- 
steifigkeit der blahenden Tragfeder entsteht bei dynamischer Belastung ein 
Schwingungssystem. Durch das Flächenverhältnis der Tragfeder zum deut- 
lich geringeren Kanalquerschnitt entsteht ein hydraulisches Übersetzungs- 
verhältnis. Dies führt dazu, dass die kleine Flüssigkeitssäule im Resonanz- 
fall wie eine Masse von mehreren Kilogramm wirkt [SA85]. 

Dieses Resonanzsystem kann bei entsprechender Abstimmung dazu genutzt 
werden, Dämpfung für den Motor-Getriebe-Verbund bei dessen Starrkör- 
permoden bereitzustellen. Die schwingende Masse im Kanal wirkt dabei 
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Abbildung 2.22: Darstellung eines hydraulisch dämpfenden Motorlagers mit Leer- 
laufschaltung 


als Tilgermasse fiir das Aggregat. Es wird daher tiblicherweise auf einen 
Frequenzbereich zwischen ca. 6 und 16 Hertz ausgelegt, wo sich die Aggre- 
gatemoden in der Regel befinden [Tre15, Hem02]. 


Eigenformen und Einflussfaktoren 


Neben der zumeist in Personenkraftfahrzeugen eingesetzten, im vorheri- 
gen Abschnitt erläuterten, Bauform von hydraulischen Topflagern existiert 
noch eine große Anzahl von Sonderbauformen. Diese übernehmen zumeist 
zusätzliche Anforderungen oder sind auf spezielle Anwendungsgebiete zu- 
geschnitten [Trel5]. Zumeist zeichnen sich diese Sonderbauformen über 
spezielle Hydrauliksysteme aus, bei denen über das Zu- bzw. Abschalten 
von Hydraulikkanälen ein spezielles Verhalten erzielt werden kann. Exem- 
plarisch sollen diese Sonderbauformen an einem hydraulisch dämpfenden 
Lager mit Leerlaufschaltung gezeigt werden, welches im Abbildung 2.22 in 
einer Schnittdarstellung abgebildet ist. 

Der Aufbau ähnelt stark dem im vorherigen Abschnitt erläuterten hydrau- 
lisch dämpfenden Motorlager, weist aber zusätzlich einen elektrischen Ak- 
tuator zum Öffnen eines Hydraulik-Bypass-Kanals auf. Durch diesen By- 
passkanal kann durch das gleiche Funktionsprinzip, welches dem zuvor 
beschriebenen Dämpfungskanal zu Grunde liegt, ein Resonanzsystem ge- 
schaffen werden. Dieses Resonanzsystem kann mittels des Aktuators zuge- 
schaltet werden. Es kann dazu verwendet werden, bei einer abstimmbaren 
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Frequenz die Transfersteifigkeit des Bauteils abzusenken. Diese Absenkung 
der Steifigkeit wird üblicherweise auf die Frequenz der dominanten Motor- 
ordnung bei Leerlaufdrehzahl abgestimmt. Hierdurch wird das Leerlaufver- 
halten von Verbrennungsmotoren besser von den Fahrzeuginsassen entkop- 
pelt [ZEF* 13]. 

Die Auswirkungen dieser Sonderbauformen beschränken sich jedoch übli- 
cherweise auf den Frequenzbereich zwischen 15 und 50 Hz. Da das schalt- 
bare Resonanzsystem in dem im Rahmen dieser Arbeit betrachteten, aku- 
stisch relevanten Frequenzbereich überkritisch ist, zeigt sich kein gesonder- 
ter Einfluss dessen auf das Körperschallübertragungsverhalten. Die Unter- 
suchungen der vorliegenden Arbeit zum Körperschallübertragungsverhal- 
ten beschränken sich daher auf die im vorherigen Abschnitt beschriebenen 
Bauformen. Die erarbeiteten Ergebnisse sind jedoch auch auf die Sonder- 
bauformen übertragbar. 


(Quasi-)Statische und niederfrequente Eigenschaften 


Die Eigenschaften eines hydraulischen dämpfenden Aggregatelagers im 
Frequenzbereich bis etwa 50 Hz sind bereits vielfach untersucht worden 
und es existieren mechanische Ersatzmodelle sowohl zur Auslegung dieser 
Bauteile als auch zur Abbildung des Systemverhaltens in Mehrkörpersimu- 
lationen. Im Folgenden soll daher dieses Verhalten erläutert werden und 
ein Überblick über den Stand der Forschung zu diesen Bauteilen gegeben 
werden. Abbildung 2.23 zeigt Betrag und Phase der komplexen Transfer- 
steifigkeiten eines hydraulischen Motorlagers bei einer Sinusanregung mit 
verschiedenen Wegamplituden zwischen 1 und 30 Hz. Im Frequenzbereich 
bis fünf Hertz kann erkannt werden, dass mit abnehmender Anregungsam- 
plitude der Betrag der komplexen Steifigkeit zunimmt. Dies lässt sich eben- 
so wie beim Verhalten des Buchsenlagers mit dem Payne-Effekt begründen, 
da eine abnehmende Anregungsamplitude auch geringere Dehnungen im 
Elastomerwerkstoff hervorruft. 

Ab etwa fünf Hertz haben die hydraulischen Komponenten einen starken 
Einfluss auf das Bauteilverhalten. Bei circa 15 Hz stellt sich ein Maximum 
in der Phase der dynamischen Steifigkeit ein. Dies hat seine Ursache in 
dem Resonanzverhalten des Fluides im Hydraulikkanal. Diese Fluidmasse 
wird dann überkritisch betrieben und koppelt daher im höheren Frequenz- 
bereich ab und sorgt somit dafür, dass sich kein volumenausgleichender 
Fluidstrom durch den Kanal mehr einstellen kann. Infolgedessen kommt es 
zu einem starken Anstieg des Betrages der dynamischen Steifigkeit ober- 
halb dieser Resonanzfrequenz. Bei einer Anregung mit einem Schwingweg 
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Abbildung 2.23: Transfersteifigkeiten eines hydraulischen Motorlagers bei Anre- 
gung mit verschiedenen Wegamplituden tiber der Frequenz 


von 0,1 mm an der Eingangsseite des Bauteils ist diese Überhöhung deutlich 
geringer ausgeprägt. Dies hat seine Ursache darin, dass für diese sehr kleine 
Amplitude ein Fluidstrom über das Membransystem möglich ist. Dies führt 
dazu, dass das Hydrauliksystem für diese Amplitude nahezu komplett ent- 
koppelt ist. 

Je größer der Schwingweg an der Eingangsseite wird, desto geringer wird 
der Anteil, welcher über das Membransystem entkoppelt ist. Daher steigt 
die Steifigkeit oberhalb der Resonanz des Kanals mit steigender Amplitude 
zunächst an. Dieses Verhalten wird jedoch vom Payne-Effekt der Tragfeder 
überlagert. Da der Payne-Effekt eine Abnahme der Steifigkeit der Elasto- 
merkomponenten mit zunehmender Amplitude bewirkt, kommt es bei den 
dargestellten Übertragungssteifigkeiten zu einem Steifigkeitsmaxium bei 
einer Anregungsamplitude von 0,3 mm. Dieses Verhalten ist von Spurk & 
Andrä [SA85], Yu et al. [YNDO1], Zhang & Shangguan [ZS06] und Lee 
& Singh [LS08] ausgiebig untersucht. Pfeffer [PHO2] und Lee [LS08] ha- 
ben mechanische Ersatzmodelle erarbeitet, welche das Verhalten der Bau- 
teile abbilden. Tao [TZXZ15] zeigt einen Ansatz der Systemmodellierung 
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mit Hilfe fraktionaler Zeitableitungen. Shangguan & Lu [SLO4] verwenden 
einen Ansatz zur Beschreibung des Bauteilverhaltens, der auf einer gekop- 
pelten Fluid-Struktur-FE-Simulation basiert. Diese Modelle befassen sich 
jedoch ausschließlich mit dem Verhalten der Bauteile im tieffrequenten Be- 
reich bis etwa 50 Hz. Geisberger [Gei00], Colgate et al. [CCC*95] und Tro- 
ge et al. [TSH* 15] haben das Verhalten der Bauteile bei Frequenzen bis ca. 
200 Hz charakterisiert und Modelle zum Verhalten dieser entwickelt, wel- 
che aber nur für Sonderformen der hydraulischen Lager valide Ergebnisse 
liefern. 


2.3.4 Herstellung von Elastomer-Verbundbauteilen 


Da die Elastomerkomponenten sowohl bei Buchsenlagern als auch bei hy- 
draulisch dämpfenden Lagern einen großen Einfluss auf die Bauteileigen- 
schaften des Aggregatelagers haben, soll in diesem Abschnitt auf die Her- 
stellung dieser Bauteile eingegangen werden. Die Prozessführung bei der 
Herstellung von Elastomeren hat einen großen Einfluss auf das Verhalten 
des Werkstoffs. Da der Großteil aller Aggregatelager noch immer auf Ba- 
sis von Naturkautschuk gefertigt werden, werden die Vernetzungsreaktionen 
und die Vulkanisation auf Basis dieses Werkstoffs erläutert. Die Herstellung 
von Naturkautschuk besteht im Wesentlichen aus der Mischungsherstellung 
und -aufbereitung und der anschließenden Vernetzungsreaktion, der Vulka- 
nisation [Tre15, CMOS]. 


Mischungsherstellung und -aufbereitung 


Der Ausgangsstoff fiir nahezu alle Naturkautschuk-Elastomere ist der Feld- 
latex der Havea-Bäume. Dieser ist eine Dispersion von cis-1,4-Polyisopren 
in einer wässrigen Phase [RM05]. Dem Feldlatex werden in weiteren Verar- 
beitungsschritten durch Trocknung und Zugabe von Enzymen Wasser und 
Pflanzenproteine entzogen [Trel5]. Abbildung 2.24 zeigt die chemische 
Struktur des Isopren-Monomers und des cis-1,4-Polyisopren-Polymers, wel- 
ches in der späteren Vernetzungsreaktion zum Elastomernetzwerk vulkani- 
siert wird. Die Kettenlänge des Polyisoprens ist im Naturkautschuk so hoch, 
dass sich der Naturkautschuk nicht ohne einen vorherigen Prozessschritt 
weiterverarbeiten lässt. Ein von Hancock [Han57] entwickeltes Verfahren 
sorgt durch hohe mechanische Scherkräfte während eines Knetprozesses für 
einen Abbau der Kettenlänge. Während dieser Mastikation werden die lan- 
gen Polyisopren-Ketten zerrissen. 

Die Mastikation dient auch dazu, den Rohkautschuk mit den restlichen Mi- 
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Abbildung 2.24: Chemische Struktur von Isopren und cis-1,4-Polyisopren 


schungsbestandteilen zu vermengen und die weiteren Mischungsbestand- 
teile homogen im Kautschuk zu dispergieren [WMO5, Tre15, Han57]. Ei- 
ne vollständig Homogenität der Mischung ist jedoch nicht erreichbar. Da- 
her kommt es bei vulkanisierten Bauteilen zwangsläufig zu Inhomogenitä- 
ten, die durch die Mischungsherstellung beeinflusst wird. Da die Mischung 
bei der Herstellung von Elastomerbauteilen zur Aggregatelagerung im An- 
schluss in einem Injection-Moulding-Verfahren in eine Vulkanisationsform 
eingespritzt wird, werden diese Inhomogenitäten jedoch zufällig im Bau- 
teilvolumen verteilt vorliegen [Tre15]. 

Die Mischung umfasst neben dem abgebauten Naturkautschuk noch eine 
Vielzahl, je nach Mischung variierender, Bestandteile [RWMOS, Tre15]: 


e Füllstoffe: Ruß, Kieselsäure, Kaolin, Glimmer, Talkum, Kalk, Titan- 
dioxid 


e Weichmacher: Ester, niedermolekulare Polyethylene 
e Oxidations-, Ozon- und Alterungsschutzmittel 
e Verarbeitungshilfen: Klebrigmacher, Aktivatoren 


e Vulkanisationssystem 
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Die Füllstoffe haben zum einen die Aufgabe, die mechanischen Eigenschaf- 
ten des Vulkanisates gezielt einzustellen (aktive Füllstoffe), zum anderen 
dienen sie oft auch als verbilligende Zusätze, um den Volumenanteil des 
Naturkautschuk zu reduzieren (passive Füllstoffe). Die Weichmacher dienen 
dazu, die Glasübergangstemperatur zu senken und bei niedrigen Temperatu- 
ren die Kristallisationsfähigkeit des Vulkanisates zu reduzieren [RWMOS]. 
Aufgrund der Doppelbindungen in den Polyisoprenketten kommt es bei Na- 
turkauschukvulkanisaten zu einer Oxidation, welche durch geeignete Oxi- 
dationschutzmittel verhindert werden soll. In der Regel kommen hier Wach- 
se zum Einsatz, welche aus dem Vulkanisat an die Oberfläche diffundieren 
und so eine passivierende Schicht bilden [Hem02]. Die Verarbeitungshilfen 
dienen dazu, die Verarbeitungsviskosität der Rohmischungsmasse so ein- 
zustellen, dass diese bestmöglich in die Vulkanisationsformen eingebracht 
werden können. 

Der wichtigste Teil der Kautschukmischung ist das Vulkanisationssystem. 
Dieses umfasst neben elementarem Schwefel auch Vulkanisationsbeschleu- 
niger und/oder Vulkanisationsverzögerer sowie Schwefeldonoren. 

Auf den Vulkanisationsvorgang und die Aufgabe der Schwefeldonoren soll 
im nächsten Abschnitt genauer eingegangen werden. 


Vulkanisation 


Die Polymerisation des Naturkautschuks geht zurück auf das von Charles 
Goodyear entwickelte Verfahren der Vulkanisation von Naturkautschuk mit 
Schwefel [Goo53]. Nach wie vor wird auch heute der während der Masti- 
kation abgebaute Naturkautschuk mit Schwefel vermengt und unter hohem 
Druck und Temperatur ausvulkanisiert. Die während dieser Vulkanisation 
ablaufenden Reaktionen haben einen entscheidenden Einfluss auf die Eigen- 
schaften des späteren Vulkanisates und sollen daher im folgenden erläutert 
werden. 


Vernetzungsreaktionen der Schwefelvernetzung 

Die exakte Beschreibung der während der Vulkanisation von Naturkatu- 
schuk auftretenden Vernetzungsreaktionen sind vielfach Gegenstand der 
aktuellen Forschung und bis heute noch nicht vollumfänglich verstanden 
und beschrieben (vgl. [MEE05, Ree82]). Zur Veranschaulichung der Ver- 
netzungsreaktion soll daher im Folgenden exemplarisch die Substitutions- 
reaktion der Polyisoprenketten näher betrachtet und erläutert werden. Der 
Ablauf dieser Vernetzungsreaktion zwischen den Polyisopren-Ketten des 
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Naturkautschuks und dem Schwefel ist in Abbildung 2.25 dargestellt. 
Unter der erhöhten Temperatur der Vulkanisation bilden sich aus dem Ele- 


CH; CH; 
C— CH Sx: C—=CH + HS,: 
C A a / 
erh BANA any CH CH, nn 
Polyisopren 
Sx 
a 
C == CH 
Ann~ CH 
Sx 


iv CH CH, nn 
Sx 
| CH3 
CH—C 
/ 
YW CH2 CH; nn 


Abbildung 2.25: Chemische Reaktion bei der Vulkanisation von Naturkautschuk 
mit Elementarschwefel nach [CM05] 


mentarschwefel freie Schwefelatome, welche mit den Methylengruppen der 
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Polyisopren-Ketten reagieren. Hierdurch kommt es zur Abspaltung eines 
Wasserstoffatoms, was dazu führt, dass die verbleibende Methingruppe ein 
freies Elektron besitzt. Dieses kann wiederum mit weiteren verfügbaren 
Schwefelatomen oder Molekülen reagieren. Die reaktive S, Gruppe besitzt 
wiederum ein freies Elektron, welches in der Lage ist die Doppelbindung 
eines Isopren-Moleküls aufzuspalten und mit dem darin enthaltenen Koh- 
lenstoffatom eine Schwefelbrücke zu bilden [CMO5]. Durch diese Schwe- 
felbrücken kann eine Polyisopren-Kette mit einer anderen oder sich selbst 
verknüpft werden. In Abbildung 2.25 sind für die beiden beteiligten Polyi- 
sopren zwei verschiedene Reaktionsmechanismen erkennbar. Die Aufspal- 
tung der Doppelbindung und die Ausbildung der Schwefelbrücke an dieser 
(im dritten Reaktionsschritt hinzugekommenen) Isoprengruppe ist eine Po- 
lyaddition. Die Reaktion, welche die andere Isoprengruppe durchläuft, ist 
hingegen eine Substitutionsreaktion, da ein Wasserstoffatom abgespalten 
wird und durch eine Schwefelbindung ersetzt wird [SH13]. Für die Aus- 
bildung dieser Schwefelbrücken sind noch weitere mögliche Positionen in 
einer Isoprengruppe für die Substitutionen von Wasserstoffatomen denkbar. 
Die Positionen, an denen eine solche Abspaltung möglich ist, sind in Abbil- 
dung 2.26 dargestellt. 


C == CH 
ww CH CH2 ww 
| | 
(a) (b) 


Abbildung 2.26: Möglichkeiten der Substitution von Wasserstoffatomen zur Aus- 
bildung von Schwefelbindungen nach [SH13] 


Für jede Schwefelbrücke, welche sich im Netzwerk ausbildet, ist eine belie- 
bige Kombination von Polyadditionsreaktion und/oder Substitutionsreakti- 


on denkbar. Die Art der Schwefelbrücken, welche sich in einem Vulkanisat 
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ausbilden, hängt dabei von dem eingesetzten Vulkanisationssystem und der 
Vulkanisationstemperatur ab [Röt13, SH13]. 

Die Vulkanisationssysteme bilden oft Kautschuk-Beschleuniger-Zwischen- 
systeme, wobei die Beschleuniger zunächst über eine Schwefelbrücke mit 
einer Isoprengruppe reagieren und sich anschließend so abspalten, dass ein 
reaktives Schwefelatom an dieser verbleibt. Um die Konzentration von reak- 
tiven Schwefelionen während der Vulkanisation zu erhöhen, kommen auch 
Schwefeldonoren zum Einsatz, welche bei den erhöhten Temperaturen, un- 
ter denen die Vulkanisation stattfindet, Schwefel abspalten[Tre15]. 

Die Bereitstellung von Schwefel hat einen entscheidenden Einfluss auf 
die Stabilität der Schwefelbrücken. Besteht diese nur aus einem einzigen 
Schwefelatom, welches die beiden Ketten direkt miteinander verbindet, 
wird diese als monosulfidische Schwefelbrücke bezeichnet, welche eine 
hohe Temperaturstabilität aufweist. Enthält die Schwefelbrücke zwei oder 
mehrere Schwefelatome, wird in der Literatur von di- bzw. polysulfidischen 
Brücken gesprochen, welche sich unter Temperatur zersetzen können. Dies 
führt zu einer Schädigung der Elastomermatrix und einer Abnahme der Ver- 
netzungsdichte [Trel5, CM05, SH13]. 


Zeitlicher Verlauf der Vernetzungsreaktion 

Die chemischen Reaktionen bei der Vulkanisation des Elastomers eine be- 
stimmte Reaktionszeit benötigen, um die Vernetzungen der einzelnen Poly- 
merketten auszubilden, soll im Folgenden der zeitliche Verlauf dieser Reak- 
tionen vorgestellt werden. Der wichtigste Parameter, welcher die mechani- 
schen Eigenschaften des späteren Vulkanisates bedingt, ist die Vernetzungs- 
dichte des Polymers. Diese kann über eine Vulkameterprüfung während der 
Vulkanisation des Werkstoffs ermittelt werden, was somit Aufschluss über 
den zeitlichen Verlauf der Vernetzungsreaktion gibt [DINO4]. Das Verfah- 
ren zur Ermittlung der Vernetzungsdichte basiert auf der Proportionalität 
zwischen dem Schubmodul des Vulkanisates und der Vernetzungsdichte bei 
kleinen Deformationen [Röt13]. Vulkameter bestehen aus einer beheizbaren 
Reaktionskammer, einem Aktor zur zyklischen Deformation der Probe und 
einem Kraft- bzw. Drehmomentmessglied. Die Kautschukmischung wird in 
die Reaktionskammer gegeben und dort bei einer konstanten Temperatur 
vulkanisiert. Die Geometrie der Probekammer wird bewusst so gewählt, 
dass sie ein kleines Volumen und eine große Oberfläche aufweist, um ei- 
ne homogene, konstante Temperatur während der Vulkanisation zu gewähr- 
leisten. Während dieses Vulkanisationsverlaufes wird die vom Vulkanisat 
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übertragene Kraft gemessen. Über den zeitlichen Verlauf der gemessenen, 
harmonisch schwankenden Kraft wird eine Hüllkurve gelegt. Die sich er- 
gebende Hüllkurve wird als Vernetzungsisotherme bezeichnet [DINO4]. Ein 
Beispiel für eine solche Vernetzungsisotherme ist in Abbildung 2.27 dar- 
gestellt. Bei Betrachtung der Reaktionsisothermen fällt auf, dass das vom 
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! ti i 
g" ty *i 
Ma: Minimalmoment; M&: End-Moment; 
Mr: Übertragenes Moment zum Zeitpunkt tx; 
ti: Inkubationszeit; t; Umsatzzeit/Vulkanisationszeit; 


ty — ti Reaktionszeit 


Vulkanisationszeit t 


Abbildung 2.27: Vernetzungsisotherme einer Kautschukmischung ohne Reversion 
nach [Röt13] 


Vulkanisat übertragene Moment zunächst abnimmt und den Minimalwert 
M, annimmt. Diese Abnahme ist mit dem Erwärmen der eingefülltem Mi- 
schung in der beheizten Form zu begründen, wodurch die Viskosität zu- 
nächst abnimmt. Im Anschluss daran findet die eigentliche Vernetzungsre- 
aktion statt, was zu einem starken Anstieg der Vernetzungsdichte im Poly- 
mer und somit auch des Übertragenen Moments führt. Kommt es nicht zu ei- 
ner konkurrierenden Reaktion oder zu einem thermischen Abbau des Netz- 
werkes, so nähert sich die Vernetzungsisotherme der horizontalen Asym- 
ptote mit M = M» [R6t13]. Die Zeit bis zum Beginn der Reaktion, wird als 
Inkubationszeit t; bezeichnet. Sie wird über den Schnittpunkt der Tangente 
am Wendepunkt der Vernetzungsisothermen mit der horizontalen Geraden 
durch das Minimum der Vulkanisationsisothermen bestimmt [DINO4]. 
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Bei der Vernetzung mit Schwefel kommt es simultan zur Vernetzungsreak- 
tion zu einem thermischen Abbau der Schwefelbrücken. Dies hat zur Folge, 
dass bei einer längeren Vulkanisationsdauer die Vernetzungsdichte wieder 
abnimmt [Röt13, Tre15]. Abbildung 2.28 zeigt den Verlauf der Vernetzungs- 
isothermen für einen Vulkanisationsvorgang mit einem Reversionsverhal- 
ten. Es ist zu erkennen, dass im Vergleich zur Vernetzungsisothermen ohne 


M, max 


M; 


= Ubertragenes Moment M 


a 


Vulkanisationszeit t 


tx 
Ma: Minimalmoment; M,max: Maximalmoment; 
M;: Übertragenes Moment zum Zeitpunkt tx; 
ti: Inkubationszeit; ty Umsatzzeit/Vulkanisationszeit; 
ty —t; Reaktionszeit 


Abbildung 2.28: Vernetzungsisotherme einer Kautschukmischung mit Reversion 
nach [R6t13] 


Reversion keine horizontale Asymptote mehr existiert, sondern sich ein Ma- 
ximum in dem übertragenen Moment Mmax und somit in der Vernetzungs- 
dichte einstellt. Im Maximum ist die Bildungsgeschwindigkeit von Vernet- 
zungsstellen gleich der Zerfallsgeschwindigkeit. Am Beispiel der Schwe- 
felvulkanisation von Naturkautschuk zerfallen nach überschreiten des Ma- 
ximums mehr Schwefelbrücken als neue gebildet werden [R6t13, CM05]. 
Dies betrifft vor allem die polysulfidischen Brücken [Tan91, MLMD13]. 
Die Untersuchungen von Mansilla et al. ['MMMS15] zeigen, dass das Rever- 
sionsverhalten von der Vulkanisationstemperatur abhängt. Abbildung 2.29 
zeigt die Vernetzungsisotherme einer Naturkautschukmischung über der Re- 
aktionszeit bei zwei verschiedenen Vulkanisationstemperaturen (433 K und 
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443 K). In dieser Abbildung zeigt die Zeitachse, nicht wie zuvor die Vulka- 
nisationszeit, sondern die Reaktionszeit, da auch die Inkubationszeit von der 
Temperatur abhängig ist [Tre15]. Die Darstellung über der Reaktionszeit er- 
möglicht daher eine bessere Vergleichbarkeit des Verlaufs der Isothermen. 
Es ist zu erkennen, dass der Gradient des übertragenen Momentes bei der 
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Abbildung 2.29: Vernetzungsisotherme einer NR-Kautschukmischung bei zwei 
Verschiedenen Temperaturen nach [MMMS15] 


Vulkanisationstemperatur von 443 K höher ist und der Anstieg zeitlich frü- 
her einsetzt. Daraus kann geschlussfolgert werden, dass die Reaktionsge- 
schwindigkeit höher ist. Dies lässt sich mit der Arrhenius-Gleichung er- 
klären, die besagt, dass chemische Reaktionen bei erhöhten Temperaturen 
schneller ablaufen [Arr89]. Ebenso ist zu beobachten, dass das Maximum 
für die Vulkanisation mit 443K geringer ist und das Reversionsverhalten 
ausgeprägter ist, als bei Vulkanisation mit 443 K. Röthemeyer und Sommer 
[SH13, R6t13] begründen dies mit der Aktivierungsenergie der Abbaureak- 
tionen der Schwefelbrücken. Diese ist deutlich größer als die der Vernet- 
zungsreaktion. Dadurch wird der Reversionseffekt durch ein Absenken der 
Vulkanisationstemperatur geringer. 

Dies deckt sich mit den Untersuchungen von Milani et al. [MLMD13] & 
Leroy et al. [LSD13], worin an NR-Mischungen ab einer Vulkanisations- 
temperatur von etwa 403K ein Reversionsverhalten beobachtet wird, wel- 
ches mit der Temperatur stark zunimmt. 

Für den industriellen Herstellungsprozess hat diese Temperaturabhängigkeit 
der Reversion einen enormen Einfluss, da aus wirtschaftlichen Gründen die 


57 


2 Stand der Forschung und Technik 


Vulkanisationszeit gering gehalten werden soll. Ebenso ist die Temperatur 
vor allem bei Bauteilen mit großen Wandstärken über den Bauteilquerschnitt 
nicht homogen. Dies kann zu prozessbedingten Inhomogenitäten im Bauteil 
führen, welche sich mit den Inhomogenitäten, die bei der Mischungsherstel- 
lung entstehen, überlagern. Im Gegensatz zu diesen, liegen die prozessbe- 
dingten Inhomogenitäten jedoch nicht stochastisch im Bauteil verteilt vor, 
sondern sind bedingt durch die Temperaturverteilung während des Vulkani- 
sationsprozesses im Bauteil lokal vorhanden. 

Die erhöhte Temperatur während der Vulkanisation liegt im Inneren des 
Bauteils aufgrund der thermischen Isolation durch die Randschichten län- 
ger an als am Rand, was Reversionseffekte fördert. Es ist daher damit zu 
rechnen, dass im inneren eines Bauteils die Vernetzungsdichte geringer ist 
als in den Randschichten. Im Bauteil enthaltene thermische Trägheiten in 
Form von metallischen Einlegern oder Anbindungsteilen, wie Kernen, Au- 
Benringen, etc. sorgen dafür, dass auch an den Anbindungsstellen zu die- 
sen Bauteilen die Temperatureinwirkung länger andauert als an Flächen, die 
nach der Vulkanisation durch Luftkontakt auskühlen können. 


2.4 Konkretisierung der Forschungsfrage 


Die in den vorherigen Abschnitten dargelegten Grundlagen zum Stand der 
Forschung und Technik zeigen, dass die Körperschalltransfereigenschaften 
von Aggregatelagern aus Elastomerkomponenten durch mechanische Vier- 
pole beschrieben werden können. Insbesondere die (komplexe) Transfer- 
steifigkeit dieser Lagerelemente stellt einen maßgeblichen Einfluss auf das 
Antriebsgeräusch im Fahrzeuginnenraum dar. Die in Abschnitt 1.2 vorge- 
stellten Forschungsfragen lassen sich also konkretisieren zu: 


e Von welchen Parametern hängt die dynamische Transfersteifigkeit 
von konventionellen Buchsenlagern und hydraulischen Topflagern 
ab? 


e Welche Komponenten des Lagerelementes beeinflussen die dynami- 
sche Transfersteifigkeit wie stark und in welchem Frequenzbereich? 


e In welchen Grenzen und durch welche Stellgrößen lässt sich die dy- 
namische Transferfteigkeit beeinflussen? 


e Wie kann ein Lagerelement auf eine gewünschtes Verhalten ausge- 
legt werden, um es an die Fahrzeugrandbedingungen (Anregung des 
Antriebsaggregates, Empfindlichkeiten der Karosseriestruktur) anzu- 
passen? 
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Um eine gezielte Auslegung des Körperschallübertragungspfades vom Trieb- 
strang bis zum Fahrerohr durchführen zu können, ist eine experimentelle Er- 

mittlung der Vierpolparameter der daran beteiligten Systeme erforderlich. 

Die Übertragungseigenschaften der Karosseriestruktur können sowohl ex- 

perimentell, als auch numerisch bereits bis in den akustisch relevanten Fre- 

quenzbereich bestimmt werden [KH93, Fre92]. Die Ermittlung des Über- 

tragungsverhaltens von elastischen Lagerelementen in diesem Frequenz- 

bereich hingegen ist kaum Gegenstand wissenschaftlicher Untersuchungen 

und auch in der Praxis nicht verbreitet. Im folgenden Kapitel soll daher auf 
die Grundlagen der Bestimmung dieser Parameter eingegangen werden und 

wichtige Einflussgrößen auf das Übertragungsverhalten auf Basis experi- 

menteller Studien erarbeitet werden. 


3.1 Theoretische Grundlagen der experimentellen 
Ermittlung von Vierpolparametern 


Zur experimentellen Ermittlung der Vierpolparameter eines elastischen La- 
gerelementes sind drei wesentliche Grundelemente notwendig. Zuerst ist 
dies das Lagerelement selbst, eine Schwingungsquelle auf der Eingangsseite 
und eine Empfängerstruktur auf der Ausgangsseite [DIN11a]. Die Schwin- 
gungsquelle dient dazu, eine Deformation auf der Eingangsseite des zu cha- 
rakterisierenden Systems einzuleiten. Auf der Empfängerseite kann dadurch 
die Systemantwort auf diese Anregung bestimmt werden. Abbildung 3.1 
zeigt ein Blockdiagramm eines solchen Systemaufbaus. 

Bei einer Modellierung dieses Systems mit Hilfe von Vierpolen ergibt sich 
für den Vierpol des elastischen Elements in Steifigkeitsform (vgl. Abschnitt 
2.1.2) 


= (3.1) 
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Fi Fy 


Schwingungs- elastisches Empfängerstruk- 
quelle x] Element X2 tur 


mit: F1:Schnittkraft an der Eingangsseite; x;:Schwingweg an der 
Eingangsseite; F:Schnittkraft an der Ausgangsseite; x2:Schwingweg an 
der Ausgangsseite 


Abbildung 3.1: Blockdiagramm des Systems Schwingungsquelle, elastisches Ele- 
ment, Empfängerstruktur nach [DINO06] 


Wird die Empfängerstruktur des Prüfaufbaus derart ausgelegt, dass die 
Schwingwege auf der Ausgangsseite klein sind, sodass x2 << xı ist, so 
kann angenommen werden, dass x2 ~ 0 ist. Gleichung 3.1 kann unter dieser 
Annahme vereinfacht werden zu: 


Fy ĉi £12) |% C1141 
= = (3.2) 
Fy © Ca} | 0 Cy 1X] 


Es kann somit aus dem eingangsseitigen Schwingweg, der eingangsseitigen 
Schnittkraft sowie der ausgangsseitigen Schnittkraft die Vorwärts-Transfer- 
steifigkeit c); und die Eingangssteifigkeit c,, des elastischen Elements be- 
stimmt werden: 

Fy Fy 


Cy. = > Cy = . 
sT A £21 a 


(3.3) 


Die beiden anderen Matrixeinträge des Vierpols c], und c,, können über 
ein Vertauschen der Ein- und Ausgangsseite des elastischen Elements er- 
mittelt werden, wobei die Vorwärts- und Rückwärtstransfersteifigkeiten un- 
ter Annahme der Reziprozität identisch sein müssen. Somit kann über ein 
Vertauschen der Ein- und Ausgangsseite und einen Vergleich der ermittel- 
ten Transfersteifigkeiten auch die Qualität und Validität des Messergebnisse 
überprüft werden. 

Da die Größen F,, F, und x, in der Regel als zeitdiskrete Messdaten vorlie- 
gen, die Übertragungseigenschaften hingegen im Frequenzbereich beschrie- 
ben werden, ist eine Ermittlung der Vierpolparameter über die Fourier- 
Transformierten der Ein- und Ausgangsgröße notwendig. Um eine Nor- 
mierung der Größen zu schaffen, empfiehlt sich eine Berechnung der Über- 
tragungseigenschaften des elastischen Elements mit Hilfe der spektralen 
Leistungsdichten. Böhme [Böh13] nutzt zur Ermittlung der Übertragungs- 
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funktion T (f) eines linearen, zeitinvarianten Systems den Quotienten aus 


Auto- und Kreuzleistungsdichtespektrum aus Gleichung 3.4: 


Py 
Fy) = a 


Hierbei sind P,,(f) das Kreuzleistungsdichtespektrum und P,x(f) das Au- 
toleistungsdichtespektrum, welche nach dem Wiener-Chintschin-Theorem 
[Chi34] über die Kreuz- bzw. Autokorrelationsfunktion (r,, bzw. ryx) für 
den Fall kontinuierlicher Zeitsignale dargestellt werden kann durch: 


. (3.4) 


Pa(f) = i 2 ra (te Pde (3.5) 
Pah) =f tls Par G.6) 
mit: 
il Tp /2 
nl) = im = f „Ost + dat 3.7) 


O 1 mR. 
nalt) = lim — f Host + rar (3.8) 
7 Tr Tr —Tr /27 


Hierbei ist x(t) das Zeitsignal am Eingang des Systems, y(t) das Zeitsignal 
am Ausgang des Systems und F(t) die komplex Konjugierte von y(t). Die 
Autokorrelationsfunktion und das Autoleistungsspektrum sind aufgrund der 
Multiplikation von x(t + T) mit der komplex konjugierten x(t) reellwertig. 
Für den Fall von zeitdiskreten Messsignalen ergeben sich die Integrale ent- 
sprechend als Summen und für die Zeitverschiebung 7 eine Verschiebung 
um eine Anzahl von diskreten Zeitschritten. Für die mathematische Herlei- 
tung sei an dieser Stelle auf Böhme [B6h13] und Geering [Gee13] verwie- 
sen. 

Die Transfersteifigkeit der untersuchten elastischen Elemente ergibt sich 
demnach zu: 

_ Pax (f ) 


ca (f) = P.. (P) (3.9) 


Da eine messtechnische Erfassung kleiner Schwingwege aufgrund der Emp- 
findlichkeit der Messsensoren oft nicht möglich ist, kann alternativ auch ei- 
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ne Messung der Schwingbeschleunigungen erfolgen. In diesem Fall ergibt 
sich zunächst eine dynamische Transfermasse. 


P ay 
m; (f) = eh 5 


Im Frequenzbereich kann diese zu einer Transfersteifigkeit integriert werden 
(vgl. Tabelle 2.3): 


(3.10) 


= Pra (f) 


aw) 
c = <. (jo (3.11) 
21 (f) Paja (f) (j ) 
Analog ergibt sich für die Eingangssteifigkeit: 
Pra (f) 2 
c = —— ‘(jo (3.12) 
u(f) Bin (f) (j ) 


3.2 Versuchsaufbau und Methodik der experimentellen 
Analyse 


Die Laborverfahren zur Messung der vibro-akustischen Transfereigenschaf- 
ten elastischer Elemente sind in DIN EN ISO 10846 erläutert [DIN11a]. 
Hierbei werden prinzipiell zwischen dem direkten und dem indirekten Ver- 
fahren zur Bestimmung der translatorischen Übertragungseigenschaften un- 
terschieden. Für die Ausführungen zu dem Prinzip nach dem direkten Ver- 
fahren sei an dieser Stelle auf Anhang A.1 verwiesen. 


3.2.1 mechanischer und messtechnischer Prüfstandsaufbau 


Bei beiden Verfahren wird das elastische Lagerelement eingangsseitig mit 
einer Deformation beaufschlagt und die Reaktionskraft an der Ausgangs- 
seite auf diese Deformation ermittelt. Das Übertragungsverhalten des zu 
charakterisierenden Elements kann bei Kenntnis der Anregung und der Re- 
aktionskraft ermittelt werden, wenn die Bewegung an der Ausgangsseite 
hinreichend klein wird [DIN1 1a, DIN1 1b, DINO6]. Das direkte und das in- 
direkte Verfahren unterscheiden sich vor allem darin, diese ausgangsseitige 
Bewegung möglichst gering zu halten. 

Da beim Prüfaufbau nach dem direkten Verfahren die obere Prüffrequenz 
durch die erste Eigenfrequenz des verwendeten Prüfrahmens begrenzt wird, 
nutzt das indirekte Verfahren eine schwingungstechnisch entkoppelte Sperr- 
masse. Abbildung 3.2 zeigt den schematischen Aufbau eines Prüfstandes 
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nach dem indirekten Verfahren. Bei dem dargestellten Aufbau wird eine 


Sly |= 


aaa 


1: höhenverstellbare Traverse; 2: Federn zur dynamischen Entkopplung 
und Einstellung der Vorlast; 3: Sperrmasse; 4: Kraftmesssystem 
(ausgangsseitig) 5: Prüfrahmen; 6: Prüfkörper; 7: Kraftmesssystem 
(eingangsseitig); 8: Federn zur dynamischen Entkopplung und Einstellung 
der Vorlast; 9: Schwingerreger 


Abbildung 3.2: Schematischer Aufbau eines Prüfstandes nach dem indirekten Ver- 
fahren [DINO06] 


Sperrmasse (3) an Federn befestigt und das sich daraus ergebende Schwing- 
system durch eine große Masse und weiche Federn so abgestimmt, dass die 
Eigenfrequenz der translatorischen Starrkörpermode möglichst tief liegt. 
Dadurch, dass die Masse des Prüflings gegenüber der Sperrmasse sehr klein 
ist und die Steifigkeit des Prüflings endlich ist, ergeben sich sehr tiefe Starr- 
körpermoden des Systems. Für Frequenzen weit oberhalb dieser Eigenfre- 
quenzen ist die Sperrmasse dynamisch entkoppelt und es gilt x2 << x). 
Als Schwingerreger wird im Rahmen dieser Arbeit ein elektrodynamischer 
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Aktor verwendet, da dieser gegenüber servohydraulischen Systemen eine 
höhere Dynamik und obere Grenzfrequenz aufweist. Zur Aufbringung einer 
Vorlast werden zusätzliche Federn (8) verwendet, da eine statische Vorlast 
durch den elektrodynamischen Schwingerregern dessen maximale Leistung 
einschränkt. Diese Federn stützen die Vorlast, welche auf der Schwingspule 
des Schwingerregers lastet auf dem Gehäuse ab. Um eine möglichst hohe Ei- 
genfrequenz der Schwingspule des elektrodynamischen Schwingerregers zu 
erreichen, sollte dieser Aufrecht stehend betrieben werden. Die zusätzlichen 
Maßnahmen, um den Schwingerreger in der Horizontalen oder hängend zu 
betreiben, senken die Eigenfrequenz der Schwingspule in der Regel ab. Da- 
her wird bei diesem Verfahren der Schwingerreger unterhalb des Prüflings 
angeordnet. 

Die Starrkörpermoden der seismischen Masse in ihrer Lagerung begrenzen 
den Frequenzbereich nach unten, in dem das indirekte Messverfahren va- 
lide Ergebnisse liefern kann. Nach oben wird der Messbereich durch die 
Frequenz der ersten Kontinuumseigenschwingung der Prüfadapter und der 
seismischen Masse begrenzt. Eine Auslegung der Sperrmasse mit möglichst 
hoher Eigenfrequenz steht dabei im Zielkonflikt mit der Forderung nach 
einer möglichst großen Masse zur Absenkung der Frequenz der Starrkör- 
permode selbiger. 

Bei dem für diese Arbeit verwendeten Prüfaufbau beträgt die Sperrmas- 
se 60kg. Die Lagerung dieser verfügt über eine Summensteifigkeit von 
100 N/mm, was zu einer Starrkörpermode bei etwa 40 Hz führt. Die Vierpol- 
Parameter der zu untersuchenden Bauteile können daher ab etwa 100 Hz 
bestimmt werden. Das indirekte Verfahren eignet sich aufgrund der hohen 
Eigenfrequenz der als Widerlager verwendeten Masse im hochfrequenten 
Bereich besser zur Charakterisierung von Lagerelementen und wird daher 
für die Ermittlung der Vierpolparameter in der vorliegenden Arbeit gegen- 
über dem direkten Messverfahren bevorzugt. 


Zur Ermittlung der Vierpolparameter nach der im vorherigen Abschnitt be- 
schriebenen Signalanalyse bedarf es der messtechnischen Erfassung ver- 
schiedener Signale sowie einer Regelung der Schwingbeschleunigung des 
elektrodynamischen Schwingerregers. Der hierzu notwendige Messaufbau 
und die Regelstrategie sollen im Folgenden erläutert werden. 

Zur Ansteuerung des elektrodynamischen Schwingerregers mit dem Ziel ei- 
ne definierte Elongation des Prüfkörpers zu erzeugen besteht die Möglich- 
keit, den Schwingweg, die Schwingschnelle oder die Schwingbeschleuni- 
gung als Führungsgröße zu nutzen. In der Praxis erweist sich die messtech- 
nische Erfassung von Schwingbeschleunigungen in Form von piezoelektri- 
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schen Beschleunigungssensoren bei kleinen Schwingwegen als zuverlässi- 
ger als eine Erfassung von Weg- oder Schnellesignalen. Daher wird für den 
für die vorliegende Arbeit verwendeten Prüfaufbau die Schwingbeschleuni- 
gung Gr als Regelgröße verwendet. Prinzipiell wäre auch eine kraftgeregel- 
te Systemanregung denkbar, was aber zur Folge hätte, dass zur Ermittlung 
der Transfersteifigkeit die Schwingbeschleunigung auf der Ausgangsseite a2 
bestimmt werden muss. Hierfür gilt jedoch bei dem gewählten Prüfaufbau 
â ~ 0, was dieses Regelungskonzept für den gewählten Aufbau untaug- 
lich macht. Um mit der gewählten Regelgröße die Transfersteifigkeit des 
Prüflings ermitteln zu können ist zusätzlich die Bestimmung der Ausgangs- 
seitigen Kraft F erforderlich. 

Bei Verwendung der Schwingbeschleunigung ar auf der Eingangsseite als 
Führungsgröße für die Regelung ergibt sich der in Abbildung 3.3 dargestell- 
te Regelkreis. Die Regelgröße ar wird dabei ständig über die Riickkopplung 
mit der Sollgröße ar verglichen. Der Regler kann somit über die Variati- 
on der Reglerübertragungsfunktion auf die Regelabweichung e = ar — ar 
reagieren. Die so ermittelte Stellgröße as ist somit Eingangsgröße der Re- 
gelstrecke, wo sowohl innere Störeffekte d; (Resonanzen des Prüfaufbaus, 
Trägheitseffekte, Nichtlinearitäten des Prüflings, etc.) als auch äußere Stör- 
effekte d, (Signalrauschen, Schwingungseintragungen aus der Umgebung, 
Störfelder, etc.) auf das System wirken. Die daraus resultierende Regelgröße 
ar wird als Ausgang der Regelstrecke erneut mit der Führungsgröße vergli- 
chen. 


Störgrößen 


di | da 
Regelabweichung Stellgröße 


Führungsgröße € a Regelgröße 
= Regler 3,1 Regelstrecke 
aF aR 
Rückkopplung 


Abbildung 3.3: Blockschaltbild zur Regelung des elektrodynamischen Schwinger- 
regers 
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3.2.2 Validierung der Prüfstandsergebnisse 


Um eine Überprüfung der Validität der Prüfstandsergebnisse durchführen 
zu können, sollen im Folgenden Abgleiche mit einem weiteren Prüfstand 
sowie der Übertragungseigenschaften im Gesamtfahrzeug erfolgen. Durch 
dieses Vorgehen kann zum einen bestätigt werden, dass die quantitativen 
Aussagen, welche durch das im vorherigen Abschnitt erläuterte Prüfverfah- 
ren generiert werden, valide sind. Zum anderen kann hierdurch auch die 
Übertragbarkeit der am Prüfstand ermittelten Transfereigenschaften in das 
Gesamtfahrzeug überprüft werden. 


Abgleich mit Prüfständen nach dem direkten Messverfahren 


Die Ermittlung der Transfersteifigkeiten von Bauteilen nach dem in An- 
hang A.l vorgestellten direkten Verfahren ist seit vielen Jahren Stand der 
Technik. Entsprechend valide sind die Ergebnisse, welche diese Prüfstän- 
de liefern. Die Ermittlung von Transfersteifigkeiten nach dem indirekten 
Verfahren ist für die Charakterisierung von Aggregatelagern bisher in der 
wissenschaftlichen Literatur nicht verwendet worden. Ebenso stellt die ex- 
perimentelle Analyse von Aggregatelagerelementen im Kilohertzbereich in 
der industriellen Praxis eine Seltenheit dar. Dies hängt auch damit Zusam- 
men, dass die Konstruktion von Adaptionen und der Aufbau eines solchen 
Prüfstandes, insbesondere die Lagerung der Sperrmasse, mit sehr hohem 
Aufwand verbunden sind. Die Lagerung der Sperrmasse stellt deshalb ei- 
ne große Herausforderung dar, da vor allem die Rotationsfreiheitsgrade der 
Sperrmasse zu Taumelbewegungen führen können, welche sich im Messbe- 
reich des Prüfstandes befinden können. Um die Qualität des Prüfstandes hin- 
sichtlich der Eignung zur Bestimmung der Transfersteifigkeiten überprüfen 
zu können, soll zunächst ein Abgleich mit einem Prüfstand nach dem di- 
rekten Messverfahren erfolgen, welcher bis zu einer Frequenz von 400 Hz 
valide Ergebnisse liefert. Mit dem im vorherigen Abschnitt vorgestellten 
Prüfstand nach dem indirekten Verfahren ergibt sich demnach der Frequenz- 
bereich zwischen 100 und 400 Hz, um die ermittelten Transfersteifigkeiten 
zwischen diesen beiden Prüfständen abzugleichen. Der Abbildung 3.4, wel- 
che die nach dem direkten und indirekten Verfahren ermittelten komple- 
xen Transfersteifigkeiten in Betrag und Phase vergleichend gegenüberstellt, 
kann entnommen werden, dass die Übereinstimmung zwischen den beiden 
verwendeten Prüfständen im überlappenden Frequenzbereich sehr hoch ist. 
Dies gilt insbesondere für den Betrag der dynamischen Transfersteifigkeit. 
Der Verlauf der Phase zeigt beim indirekten Messverfahren bei etwa 105 und 
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210 Hz eine geringfügig ausgebildete Unstetigkeit und auch größere Abwei- 
chungen von den Messwerten nach dem direkten Verfahren. Die Ursache für 
die Phasenabweichung ist bei 105 Hz eine Querschwingung der Sperrmas- 
se, welche sich nicht vollständig unterdrücken lasst. Bei 210 Hz handelt es 
sich um eine Resonanz des elektrodynamischen Schwingerregers, welche 
mit Hilfe der im vorherigen Abschnitt dargestellten Regelstrecke nicht voll- 
ständig ausgeregelt werden kann. 

Ab etwa 350 Hz nimmt die Phase bei der nach dem direkten Verfahren er- 
mittelten Transfersteifigkeit deutlich zu, was bei dem Verlauf der nach dem 
indirekten Verfahren ermittelten Phase nicht zu erkennen ist. Hier zeigt sich 
die Resonanz des Prüfrahmens des Prüfstandes nach dem direkten Verfah- 
ren, welcher die Messergebnisse beeinflusst und den Messbereich des Prüf- 
standes nach oben hin begrenzt. Zusammenfassend kann festgehalten wer- 
den, dass die Übereinstimmungen zwischen den beiden Messverfahren sehr 
hoch sind und die jeweiligen Abweichungen sich dem Aufbau der Prüfstän- 
de zuordnen lassen. Es zeigt sich auch deutlich, dass für Frequenzen ober- 
halb von 350 Hz das indirekte Messverfahren bessere Ergebnisse liefert. Da 
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Abbildung 3.4: Vergleich der ermittelten Transfersteifigkeiten eines Getriebelagers 
nach dem direkten und indirekten Verfahren 
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das indirekte Prüfverfahren im tieffrequenten Bereich bis 400 Hz valide Er- 
gebnisse liefert, soll auch die Qualität der Messergebnisse im höherfrequen- 
ten Bereich validiert werden. Ein Abgleich mit anderen Prüfständen ist auf- 
grund der Tatsache, dass diese nur in sehr begrenzter Anzahl existieren, nur 
bedingt möglich. Es soll daher im Folgenden eine Korrelation zwischen dem 
Körperschallübertragungsverhalten von Bauteilen im Fahrbetrieb zu der er- 
mittelten dynamischen Transfersteifigkeit aufgezeigt werden. 


Korrelation der Übertragungseigenschaften zwischen 
Prüfstand und Fahrzeug 


Zur Beurteilung der Korrelation zwischen Fahrzeug und Prüfstand wird eine 
sinusförmige Schwingbeschleunigung mit einer über der Frequenz konstan- 
ten Amplitude von 30 m/s? verwendet, was in etwa der mittleren Schwing- 
beschleunigung bei einem Teillastbetrieb eines Verbrennungsmotors ent- 
spricht. Um eine Korrelation des Körperschallübertragungsverhaltens von 
Bauteilen zwischen Fahrbetrieb und Prüfstand abbilden zu können, ist es 
aus Gründen der Signalerfassung am einfachsten, die Übertragungsfunktion 
der Schwingbeschleunigung im Fahrzeug mit der am Prüfstand ermittelten 
Transfersteifigkeit zu vergleichen. Abbildung 3.5 zeigt die Übertragungs- 
funktion der 7;,;, der Schwingbeschleunigung ¥ı auf der Aggregateseite 
zur Schwingbeschleunigung x» auf der Karosserieseite an den Anbindungs- 
stellen eines Getriebelagers im Fahrzeug und dessen am Prüfstand ermittelte 
Transfersteifigkeit über der Frequenz. 

Es ist zu erkennen, dass die Übertragungsfunktion dort hohe Werte auf- 

weist, wo die Transfersteifigkeit ebenfalls hoch ist. Die Körperschallüber- 
tragung ist also weitgehend proportional zu der experimentell ermittelten 
Transfersteifigkeit. Lediglich bei ca. 250 Hz zeigt die Übertragungsfunktion 
der Schwingbeschleunigungen ein Maximum, welches nicht auf die Trans- 
fersteifigkeit zurückgeführt werden kann. Die Ursache für dieses Maximum 
liegt in einem Minimum der lokalen Eingangssteifigkeit an der Empfän- 
gerseite, was dort zu hohen Beschleunigungswerten führt und somit ein 
Maximum der Übertragungsfunktion bewirkt. 
Die guten Übereinstimmungen zwischen der am Prüfstand ermittelten Trans- 
fersteifigkeit und dem Übertragungsverhalten des Bauteils im Fahrzeug 
sprechen für valide Prüfstandsergebnisse und die Möglichkeit, Bauteilop- 
timierungen in der Prüfstandsumgebung durchführen zu können. Im Fol- 
genden sollen die beiden wichtigsten Bauformen von Aggregatelagern mit 
der beschriebenen Methodik experimentell untersucht werden und wichtige 
Einflussfaktoren auf das Bauteilverhalten analysiert werden. 
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Abbildung 3.5: Darstellung der Ubertragungsfunktion der Schwingbeschleuni- 
gungen im Fahrzeug und der Transfersteifigkeit des elastischen 
Elements am Prüfstand über der Frequenz 


3.3 Ermittlung des Übertragungsverhaltens von 
Buchsenlagern 


An denjenigen Aggregatelagerpositionen, welche durch einen hochfrequen- 
ten Körperschalleintrag durch das Antriebsaggregat besonders hohe akusti- 
sche Anforderungen erfüllen müssen, kommen in der Regel einfach Buch- 
senlager zum Einsatz. Diese Lagerpositionen sind meist die Anbindungs- 
punkte des Hauptgetriebes und des Achsgetriebes bzw. Differentials. Durch 
die Zahneingriffsgeräusche von schnell drehenden Getriebewellen kommt 
es zu einem Körperschalleintrag, welcher einen Frequenzbereich bis meh- 
rere Kilohertz umfasst. Nicht hydraulisch bedämpfte Buchsenlager eignen 
sich für diesen Einsatz, da hierbei meist geringe niederfrequente dynami- 
sche Lasten abgestützt werden und die niedrige dynamische Verhärtung der 
Lager aufgrund der nicht vorhandenen Hydraulik die akustischen Übertra- 
gungseigenschaften positiv beeinflusst. Im Folgendenen sollen die Einflüs- 
se auf die Körperschallübertragungseigenschaften dieser Bauteile in der im 
vorherigen Abschnitt erläuterten Prüfstandsumgebung untersucht werden. 
Hierzu wird der Einfluss des eingesetzten Elastomerwerkstoffes, der auf die 
Bauteile wirkenden statischen Vorlast, der Bauteiltemperatur und der Ein- 
fluss des Anregungssignals untersucht. 
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3.3.1 Einfluss des Elastomerwerkstoffes 


Ein entscheidender Anteil der Körperschallübertragungseigenschaften aller 
Elastomerlager ist auf den eingesetzten Elastomerwerkstoff zurückzufüh- 
ren, da dieser die einzige mechanische Verbindung von Erreger- und Emp- 
fängerseite des Bauteils darstellt. Das Schwingungsverhalten dieses Ela- 
stomeranteils bestimmt daher maßgeblich das Bauteilverhalten. Um die- 
sen Einfluss zu untersuchen, werden vier verschiedene Nauturkautschukmi- 
schungen mit Shore-A-Härten von 45 bis 60 Sh verwendet und in der selben 
Vulkanisationsform vulkanisiert. Anschließend werden diese vier Varian- 
ten in der zuvor vorgestellten Prüfstandsumgebung analysiert. Da die ex- 
perimtentelle Analyse von Buchsenlagern eine Herausforderung bezüglich 
der Vorrichtungen zur Montage an den Prüfstand stellt, kann aufgrund von 
Resonanzen in diesen Adaptern eine experimentelle Analyse nur bis 2 kHz 
erfolgen. Die so ermittelten Transfersteifigkeiten sind in Abbildung 3.6 dar- 
gestellt. 

Es ist zu erkennen, dass die Transfersteifigkeit im betrachteten Frequenzbe- 
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Abbildung 3.6: Vergleich der experimentell ermittelten Transfersteifigkeit von 
Buchsenlagern verschiedener Naturkautschukvulkanisate 
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reich durch einen Resonanzeffekt in den untersuchten Bauteilen beeinflusst 
wird. Diese ist zum einen durch ein Maximum im Betrag der Transferstei- 
figkeit, welches von einem starken Abfall des Betrags gefolgt wird, zum 
anderen durch ein Minimum der Phase der komplexen Transfersteifigkeit 
bestimmt. Für die Varianten mit 45 und 50 Sh ist dieses Verhalten im unter- 
suchten Frequenzbereich ersichtlich. Bei der untersuchten Variante mit einer 
Shore-A-Härte von 55 Sh ist ab etwa 1750 Hz ein Abfall in der Phase zu er- 
kennen und auch ein Anstieg des Betrages. Das Maximum des Betrages und 
das Minimum der Phase liegen daher vermutlich geringfügig außerhalb des 
untersuchten Frequenzbereiches. Für die Variante mit 60 Sh liegt der Reso- 
nanzeffekt deutlich außerhalb des analysierbaren Frequenzbereiches, da der 
Betrag nahezu linear ansteigt, was lediglich die dynamische Versteifung des 
Elastomers bewirkt und die Phase über der Frequenz nahezu konstant bleibt. 
Zusammenfassend kann geschlussfolgert werden, dass das Körperschall- 
übertragungsverhalten von Buchsenlagern sich durch einen Resonanzeffekt 
kennzeichnen, der eine starke Abhängigkeit vom verwendeten Werkstoff 
aufweist. Da das Bauteilverhalten dieser einfachen Bauteile nur von den 
Eigenmoden der beteiligten Komponenten beeinflusst werden kann, deutet 
dieses Verhalten darauf hin, dass es sich bei der betreffenden Bauteilmode 
um eine Kontinuumsschwingform des Elastomeranteils handelt. 


3.3.2 Einfluss der statischen Vorlast 


Beim Einsatz von Buchsenlagern zur Lagerung eines Aggregates bei ver- 
schiedenen Fahrzuständen, insbesondere bei der Abstützung verschieden 
hoher Antriebsmomente des Triebstrangs über die Lagerelemente, ist die 
quasistatische Auslenkung der Bauteile starken Veränderungen unterwor- 
fen. Im Folgenden soll daher das Verhalten der Bauteile bei unterschiedlich- 
en statischen Vorlasten untersucht werden. 

Abbildung 3.7 zeigt die dynamische Transfersteifigkeit eines Buchsenlagers 
über der Frequenz in einem Bereich von ON bis 3000N. Das untersuchte 
Bauteil weist eine statische Steifigkeit von 1000 N/mm, entsprechend ergibt 
sich bei 3000 N eine Elongation von 3mm, was einer maximalen Dehnung 
von 178 % im Bauteil entspricht. 

Bei Betrachtung des Betrages der dynamischen Transfersteifigkeit kann 
festgestellt werden, dass sich zwischen ON und 2500N keinerlei Änderun- 
gen in der Bauteilsteifigkeit einstellen. Die minimalen Abweichungen zwi- 
schen den einzelnen Transfersteifigkeiten, welche auf den Resonanzbereich 
des Bauteils bei etwa 1400 Hz beschränkt sind, lassen sich auf eine geringe 
Erwärmung des Bauteils während des Prüfablaufs zurückführen (vgl. Ab- 
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Abbildung 3.7: Dynamische Transfersteifigkeit eines Viersteg-Buchsenlagers bei 
verschiedenen Vorlasten 
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schnitt 3.3.3). Der Betrag der dynamischen Transfersteifigkeit bei einer Vor- 
last von 3000N ist circa 7 % höher als der Betrag bei geringeren Vorlasten, 
was darin begründet liegt, dass sich das Lagerelement bei dieser Vorlast 
bereits im progressiven Anteil der quasistatischen Federkennlinie befindet 
und somit auch über eine höhere Grundsteifigkeit verfügt. Eine Ermittlung 
der dynamischen Transfersteifigkeiten bei höheren Vorlasten ist daher auch 
nicht sinnvoll, da sich das Bauteil in diesem Fall bereits im Anschlag be- 
finden würde und deutlich höhere dynamische Steifigkeit aufweisen würde. 
Dies geht mit einer deutlichen Verschlechterung des Körperschallübertra- 
gungsverhaltens einher und muss daher bei der Auslegung der statischen 
Steifigkeit und des linearen Freiwegs bei einer Grundauslegung der Lager- 
elemente berücksichtigt werden. 

Die Phase der dynamischen Transfersteifigkeit zeigt ein ebenso unverän- 
dertes Verhalten über der Vorlast wie der Betrag. Die Ausprägung des Mi- 
nimums bei circa 1500 Hz, welches durch die Bauteilresonanz verursacht 
wird, zeigt jedoch eine leichte Veränderung. Es zeigt sich, dass mit zuneh- 
mender Vorlast die Ausprägung des Minimums leicht reduziert wird und es 
zu einer geringen, kontinuierlichen Erhöhung der Phasenwerte über der Vor- 
last kommt. 

Zusammenfassend kann festgehalten werden, dass die Variation der Vorlast 
und somit der statischen Elongation des Bauteils einen sehr geringen Ein- 
fluss auf das Körperschallübertragungsverhalten der Buchsenlager hat. Dies 
gilt, solange das Bauteil sich im linearen Bereich der statischen Federkennli- 
nie befindet. Da der progressive Bereich der statischen Federkennlinie nicht 
mehr als Auslegungsbereich für Komfortlastfälle gilt, wird eine Verschlech- 
terung der vibroakustischen Transfereigenschaften in diesem Bereich ak- 
zeptiert. Es kann demnach schlussgefolgert werden, dass für die Auslegung 
maßgeblichen Komfortlastfälle keine signifikante Änderung des Übertra- 
gungsverhaltens des untersuchten Buchsenlagers stattfindet. Da die Varian- 
tenvielfalt bezüglich der Geometrie gering ist, kann dieses Verhalten ver- 
allgemeinernd auf die gesamte Bauteilklasse von Komfort-Buchsenlagern 
übertragen werden. 


3.3.3 Einfluss der Bauteiltemperatur 


Da Aggregatelagerbauteile in unmittelbarer Nähe zu den Antriebsstrang- 
komponenten im Fahrzeug verbaut sind, sind diese während des Betriebs er- 
höhten Temperaturen ausgesetzt. Der Einfluss der Temperatur auf das Werk- 
stoffverhalten von Elastomeren bedingt, dass sich das Bauteilverhalten ent- 
sprechend verändert. Die Ergebnisse aus Abschnitt 3.3.1 legen nahe, dass 
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eine Temperaturvariation einen großen Einfluss auf das Körperschallüber- 
tragungsverhalten hat. 

Die Betriebstemperatur der Aggregatelager hängt von der Umgebungstem- 
peratur des Fahrzeuges, dem Betriebspunkt des Triebstrangs, der Wärme- 
leitfähigkeiten des Motortragarmes und der Umströmung der Elastomerele- 
mente mit Kühlluft, also der Fahrgeschwindigkeit des Fahrzeuges, ab. Um 
die Körperschallübertragung durch die Aggregatelager gezielt auslegen zu 
können ist es daher erforderlich, das Bauteilverhalten in einem breiten Tem- 
peraturbereich zu betrachten. 

Abbildung 3.8 zeigt die komplexe Transfersteifigkeit eines Vier-Steg-Buch- 
senlagers bei Temperaturen zwischen 281 K und 345 K in Betrag und Phase. 
Es kann festgestellt werden, dass im Frequenzbereich bis etwa 1000 Hz der 
Betrag der Transfersteifigkeit mit zunehmender Temperatur abnimmt. Dies 
deckt sich mit den den Erwartungen, die sich aus dem Verhalten von Elasto- 
merwerkstoffen nach Abschnitt 2.2.1 ergeben. Ebenso zeigt sich, dass sich 
die Eigenmoden des untersuchten Buchsenlagers mit zunehmender Tempe- 
ratur zu niedrigeren Frequenzen hin verschieben. Dies ist bedingt durch eine 
Abnahme des dynamischen Moduls bei einer nahezu konstanten Dichte des 
Elastomers. 

Die Bauteilprüfungen mit einer definierten Temperaturregelung sind auf- 
grund der geringen thermischen Trägheit des Elastomers und der Eigen- 
erwärmung durch die Strukturdämpfung mit großem Aufwand verbunden. 
Eine Abschätzung des Bauteilverhaltens bei veränderter Temperatur kann 
daher auch über die WLF-Gleichung (vgl. Abschnitt 2.2.1) erfolgen. 

Um die Parameter Cı und C2 der WLF-Gleichung zu ermitteln, wird der 
Betrag der dynamischen Steifigkeit im Frequenzbereich zwischen 100 Hz 
und 1000 Hz durch eine lineare Regression approximiert und dann zu höhe- 
ren Frequenzbereichen extrapoliert. Tabelle 3.1 zeigt die Parameter dieser 
linearen Regression: 

Als Referenztemperatur zur Ermittlung der WLF-Parameter werden 296 K 
gewählt, da bei dieser Temperatur eine Ermittlung der Transfersteifigkeiten 
von Bauteilen am einfachsten zu realisieren sind. Weiterhin wird als Refe- 
renz derjenige Frequenzpunkt gewählt, bei dem die jeweilig Transfersteifig- 
keit einen Betrag von 4000 N/mm annimmt. Dies ist für die Referenztem- 
peratur bei fo = 3987 Hz der Fall. Somit ergibt sich die WLF-Gleichung aus 
Gleichung (2.50) zu: 


S87 _ Ci(T—296K) (3.13) 


l = 
os( C + (T -296K) 
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Abbildung 3.8: Dynamische Transfersteifigkeit eines Viersteg-Buchsenlagers bei 
verschiedenen Temperaturen 
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Temperatur inK Steigung in (N/mm)/Hz Ordinatenabschnitt in N/mm 


281 1,10 2616 
287 0,75 2278 
296 0,59 1642 
312 0,34 1361 
322 0,28 1183 
326 0,26 1165 
345 0,19 1088 


Tabelle 3.1: Daten der linearen Regression der Transfersteifigkeiten 


Die beiden Parameter lassen sich über je zwei Temperaturstützstellen er- 
mitteln. Tabelle 3.2 zeigt die ermittelten WLF-Parameter für die einzelnen 
Stützstellen und die sich daraus ergebenden Mittelwerte. 


Temperatur inK Frequenz in Hz log(#) C CoinKk 


281 1349 1,08 -2,84 54,76 
287 2272 0,56 -2,85 54,87 
296 3987 0 - - 
312 7695 -0,66 2,83 53,6 
322 10012 -0,92 2,85 54,29 
326 10932 -1,01 2,84 53,98 
345 15473 -1,36 2,86 54,73 
Mittelwerte -2,84 54,37 


Tabelle 3.2: Experimentell ermittelte Parameter der WLF-Gleichung 
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3.3 Ermittlung des Übertragungsverhaltens von Buchsenlagern 


3.3.4 Einfluss des Anregungssignals 


In Abschnitt 3.2.2 konnte gezeigt werden, dass eine Sinusbelastung mit 
konstanter Schwingbeschleunigung an dem vorgestellten Prüfstand das Ver- 
halten eines Bauteils im Fahrzeug abbilden kann. Da der Payne-Effekt (vgl. 
Abschnitt 2.2.1) für Elastomerwerkstoffe jedoch eine starke Amplitudenab- 
hängigkeit vorhersagt, soll zunächst der Einfluss der Anregungsamplitude 
ermittelt werden. 


Zur Untersuchung der Abhängigkeit der Bauteileigenschaften von der Am- 
plitude einer Sinusanregung wird die Transfersteifigkeit von Buchsenla- 
gern unterschiedlicher Naturkautschukmischungen bei einer Frequenz von 
300Hz und verschiedenen Amplituden ermittelt. Abbildung 3.9 stellt die 
Ergebnisse dieser Untersuchungen für Buchsenlager aus Naturkautschuke- 
lastomeren zwischen 40 und 60 Shore-A bei Anregungen mit einem Sinus- 
signal mit einer Wegamplitude zwischen 0,1 um und 100 um dar. Die obere 
Grenze von 100 um entspricht derjenigen Amplitude, welche bei hoher Mo- 
torlast als antriebserregter Körperschall im betrachteten Frequenzbereich 
maximal zu erwarten ist. 

Bei Betrachtung von Abbildung 3.9 ist zu erkennen, dass für alle untersuch- 
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Abbildung 3.9: Dynamische Transfersteifigkeit verschiedener 4-Steg-Buchsenlager 
bei 300 Hz und verschiedenen Anregungsamplituden/- 
Dehnungasamplituden 
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3 Experimentelle Untersuchung von Aggregatelagern 


ten Bauteile nahezu keine Abhängigkeit von der Amplitude der Anregung 
besteht. Betrachtet man die auf der zweiten Abszisse aufgetragene maxima- 
le Dehnung, die bei dieser Anregung im Bauteil auftritt, so ist diese deutlich 
kleiner als ein Prozent. Der Payne-Effekt tritt, abhängig vom Füllstoffan- 
teil, üblicherweise bei Dehnungen von über einem Prozent auf, sodass die 
sehr geringen Dehnungen dazu führen, dass die Deformationen alle unter- 
halb der kritischen Deformationsgrenze sind, ab der das Füllstoff-Netzwerk 
aufreißt und die Füllstoff-Füllstoff-Interaktionen nicht mehr gegeben sind. 
Es kann daher geschlussfolgert werden, dass die Amplitude, bei der die 
Transfereigenschaften der Bauteile an einem Prüfstand ermittelt werden 
keinen Einfluss auf die Körperschalltransfereigenschaften haben, solange 
die im Elastomerbauteil auftretenden Dehnungen deutlich unter einem Pro- 
zent sind, was für Aggregatelagerelemente aufgrund ihrer großen Bauform 
und der geringen Schwingwege bei hohen Frequenzen gegeben ist. Kommt 
es jedoch zu Überlagerungseffekten mit fahrbahnerregten Schwingungen, 
so kann diese Grenze auch überschritten werden. Bei Fahrmanövern mit 
großen, fahrbahninduzierten Schwingungen im Frequenzbereich der Ag- 
gregatemoden (circa 5 bis 15 Hz) können Auslenkungen der Aggregatelager 
von mehreren Millimetern mit entsprechend hohen Dehnungen auftreten. 
Für diese Lastfälle spielen in der Praxis triebstrangerregte Geräuschphäno- 
mene eine untergeordneten Rolle, weshalb eine akustische Auslegung auf 
diese Lastfälle nicht als sinnvoll zu betrachten ist. Für die Modellierung des 
Bauteilverhaltens in Kapitel 4 kann daher die Amplitudenabhängigkeit ver- 
nachlässigt werden, solange die Dehnungen im Bauteil hinreichend klein 
bleiben. 


3.4 Ermittlung des Übertragungsverhaltens von 
hydraulischen Topflagern 


Zur Ermittlung des Übertragungsverhaltens hydraulischer Topflager sollen 
ebenso, wie bereits im vorherigen Abschnitt, für Buchsenlager wichtige 
Einflussfaktoren ermittelt werden, um diese in der im nachfolgenden Ka- 
pitel beschriebenen Modellbildung berücksichtigen zu können. Da die Ela- 
stomerkomponenten in hydraulischen Topflagern in ihren Mischungseigen- 
schaften denen von Buchsenlagern sehr ähnlich sind, bestehen die selben 
Einflüsse der statischen Vorlast und der Anregungsamplitude, wie sie be- 
reits im vorherigen Abschnitt für diese Bauform aufgezeigt sind. Auf die 
Einflüsse der statischen Vorlast und der Anregungsform wird für die hy- 
draulischen Lager daher nicht mehr explizit eingegangen und deren Einfluss 
vernachlässigt. 
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3.4 Ermittlung des Übertragungsverhaltens von hydraulischen Topflagern 


Für die Wechselwirkungen zwischen Hydraulik und den Elastomerkompo- 
nenten ist die Elastomermischung, aus welcher diese hergestellt sind, von 
Interesse. Da die Eigenschaften des Elastomers auch, wie bereits in Ab- 
schnitt 3.3.3 gezeigt, von der Temperatur abhängen, soll auch dieser Einfluss 
untersucht werden. Darüber hinaus ist ein weiterer wichtiger Parameter für 
das Gesamtsystemverhalten der Bauteile der Einfluss der Hydraulik. Im fol- 
genden sollen daher diese drei Haupteinflussfaktoren auf ihre Auswirkung 
bezüglich der Transfersteifigkeit von hydraulischen Topflagern untersucht 
werden. 


3.4.1 Einfluss des Elastomerwerkstoffes 


Da die Elastomer-Tragfeder einer der wichtigsten Komponenten eines hy- 
draulischen Lagers darstellt, soll über eine Variation der Härte des Bau- 
teils der Einfluss dieser ermittelt werden. Hierzu dienen Tragfedern gleicher 
Geometrie, jedoch unterschiedlicher Elastomermischungen, sodass sich im 
fertigen Vulkanisat Shore-Härten zwischen 45 Sh und 62 Sh einstellen. Zur 
gezielten Variation der Elastomertragfedern bei gleichbleibenden hydrauli- 
schen Eigenschaften finden die Untersuchungen mithilfe von modifizierten 
Lagergehäusen statt. Diese lassen einen Austausch einzelner Komponenten 
und ein anschließendes Wiederbefüllen zu. So können stochastische, pro- 
duktionsbedingte Streuungen in den Eigenschaften der einzelnen Bauteile 
bei der Variation der Elastomertragfedern vermieden werden. Zur Analyse 
werden die Bauteile mit Hilfe der in Abschnitt 3.2 erläuterten Prüfumge- 
bung untersucht. Die daraus resultierenden Transfersteifigkeiten der Bautei- 
le können Abbildung 3.10 entnommen werden. 

Die untersuchten Bauteile zeigen einen für hydraulische Topflager charak- 
teristischen Verlauf, welcher sich in ähnlicher Ausprägung bei vielen hy- 
draulischen Topflagern finden lässt. Die Transfersteifigkeit ist durch zwei 
bis drei charakteristische Resonanzen geprägt, welche stark durch die Ela- 
stomermischung der Tragfeder beeinflusst werden. Diese Resonanzen zei- 
gen sich durch Maxima im Betrag der Transfersteifigkeit, welche entspre- 
chende Maxima in der Phase bedingen. Hierbei zeigt sich, dass die Abso- 
lutwerte der Transfersteifigkeit erwartungsgemäß mit einer zunehmenden 
Werkstoffhärte der Tragfeder ansteigen. Ebenso zeigt sich, dass die Bauteil- 
resonanzen durch einen härteren Tragfederwerkstoff zu höheren Frequenzen 
verschoben werden. Dies führt dazu, dass die Bauteilresonanz, welche bei 
der untersuchten Variante mit 45 und 53 Sh bei ca. 2100 Hz liegt, bei der 
Variante mit 62 Sh außerhalb des Messbereichs liegt. Demzufolge ist diese 
dritte Bauteilresonanz bei dieser Variante nicht mehr erkennbar. 
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Abbildung 3.10: Transfersteifigkeiten von hydraulischen Topflagern mit Tragfedern 
unterschiedlicher Shore-Härten 
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3.4 Ermittlung des Übertragungsverhaltens von hydraulischen Topflagern 


Der starke Einfluss des Tragfederwerkstoffes auf das Bauteilverhalten legt 
nahe, dass die Tragfedern am Schwingverhalten des Gesamtsystems betei- 
ligt sind und dieses als Federelement beeinflussen. Ein Versteifung dieses 
Federelementes durch eine Werkstoffänderung führt demnach zu einer Fre- 
quenzerhöhung des Resonanzphänomens. Bei den verwendeten Elastome- 
ren handelt es sich, wie in Abschnitt 2.2.1 erläutert und in Abschnitt 3.3.3 
für das Bauteilverhalten von Buchsenlagern untersucht, um stark temperatu- 
rabhängige Werkstoffe. Da der Elastomerwerkstoff einen signifikanten Ein- 
fluss auf das Systemverhalten zeigt, wird dieser Einfluss in einer ähnlichen 
Größenordnung auch für eine Temperaturänderung der Bauteile erwartet. 
Dieser Einfluss soll im folgenden Abschnitt näher untersucht werden. 


3.4.2 Einfluss der Bauteiltemperatur 


Wie bereits in Abschnitt 3.3.3 erwähnt, sind die Aggregatelager eines Kraft- 
fahrzeuges je nach Betriebs- und Umgebungsbedingungen des Fahrzeugs 
variierenden Temperaturen ausgesetzt. Nachdem im letzten Abschnitt der 
Einfluss des Tragkörperwerkstoffes analysiert wurde, soll im folgenden der 
Einfluss der Temperatur auf das Bauteilverhalten analysiert werden. 
Abbildung 3.11 zeigt die dynamische Transfersteifigkeit für hydraulisches 
Topflager mit einer Tragfeder aus einer Elastomermischung von 53 Sh für 
fünf verschiedene Temperaturen zwischen 288 K und 333 K. Es ist zu erken- 
nen, dass sich das gesamte Steifigkeitsniveau mit zunehmender Temperatur 
reduziert und die Bauteilresonanzen zu tieferen Frequenzen verschieben. 
Die Bauteilresonanz, welche sich bei 288K bei etwa 1600 Hz ausbildet, 
verschiebt sich mit zunehmender Temperatur bis auf etwa 1400 Hz. Die- 
ses Verhalten entspricht demjenigen, was mit der Substitution des Tragkör- 
perwerkstoffes durch eine weichere Mischung einhergeht (vgl. Abschnitt 
3.4.1). Dies deckt sich mit dem in Abschnitt 2.2.1 bereits erläuterten Zusam- 
menhang, dass eine Temperaturerhöhung bei Elastomeren eine Erweichung 
des Werkstoffes mit sich bringt. 

Ebenso bemerkenswert ist die Ausprägung einer weiteren Resonanz, wel- 
che bei tiefen Temperaturen nur geringfügig in Erscheinung tritt. Bei er- 
höhten Temperaturen bildet sich jedoch bei dem untersuchten Bauteil ein 
lokales Maximum in Betrag und Phase der Transfersteifigkeit bei etwa 
800 Hz aus. Dies bildet sich genau dort aus, wo die Tragfeder ohne hydrauli- 
sche Komponenten eine Mode aufweist (vgl. Untersuchungen in Abschnitt 
4.3.1). Dies deutet darauf hin, dass sich in der Arbeitskammer durch die 
Temperaturerhöhung eine Nachgiebigkeit in Form von Ausgasungen aus 
dem Fluid bildet. Ursächlich hierfür sind vermutlich im Hydraulikmedium 
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Abbildung 3.11: Transfersteifigkeiten eines hydraulischen Topflagers bei verschie- 
denen Bauteiltemperaturen 
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3.4 Ermittlung des Übertragungsverhaltens von hydraulischen Topflagern 


gelöste Gase, deren Löslichkeit durch die Temperaturerhöhung abnimmt. 
Da dieses lokale Maximum jedoch im Vergleich zu den beiden Maxima bei 
circa 200 Hz und 1500 Hz eine deutlich geringer Transfersteifigkeit aufweist 
und somit für die Körperschallübertragung von geringerer Relevanz ist, soll 
dieses im Folgenden nicht weiter betrachtet werden. 

Die Temperaturabhängigkeit der Transfersteifigkeit der Bauteile kann dem- 
zufolge auf die Elastomerkomponenten zurückgeführt werden und kann 
durch die Verwendung der WLF-Parameter (vgl. Abschnitt 2.2.1), welche in 
Abschnitt 3.3.3 für die verwendeten Elastomermischungen bereits ermittelt 
wurden, beschrieben werden. Hierzu ist es jedoch erforderlich, die Wech- 
selwirkungen zwischen den Elastomerkomponenten und den hydraulischen 
Komponenten zu analysieren, da letztere nicht durch die WLF-Gleichung 
beschreibbar sind. Im folgenden Abschnitt soll daher der Einfluss der hy- 
draulischen Parameter genauer identifiziert werden. 


3.4.3 Einfluss der hydraulischen Parameter 


Um den Einfluss der hydraulischen Abstimmung des Lagers auf die Trans- 
fersteifigkeiten der Bauteile zu analysieren, sollen im Folgenden verschiede- 
ne Parameter der hydraulischen Abstimmung untersucht werden. Die Aus- 
wirkung der geometrischen Eigenschaften des Kanals im Düsen-Membran- 
System ist bereits vielfach betrachtet worden (vgl. Abschnitt 2.3.3). Da das 
Schwingsystem üblicherweise bei circa 5-18 Hz abgestimmt ist, hat dieses 
für das hochfrequente Bauteilverhalten einen geringen Einfluss. Die Ein- 
flüsse der Abstimmung des Düsen-Membran-Systems sind hingegen noch 
nicht weiter untersucht. Um diese Einflüsse zu analysieren werden die Höhe 
des Gitters, zwischen dem die Membran eingelegt ist, sowie das Lochbild 
dieses Gitters variiert. Abbildung 3.12 stellt diese Parameter grafisch dar. 


Für die Gitterhöhe werden zwei Werte betrachtet (6,5 mm und 3 mm), das 
Lochbild wird in vier verschiedenen Stufen (12, 21, 30 und 36 Löcher) vari- 
iert. Die sich daraus ergebenden Varianten 1-8 sind in Tabelle 3.3 aufgeführt. 
Die sich bei der experimentellen Untersuchung der Varianten 1-8 ergeben- 
den komplexen Transfersteifigkeiten der hydraulischen Aggregatelager sind 
in Abbildung 3.13 und 3.14 dargestellt. 


Es ist zu erkennen, dass sich im dargestellten Frequenzbereich bis 1000 Hz 
bei allen untersuchten Varianten ein Resonanzeffekt einstellt, welcher zu ei- 
nem Anstieg des Betrages der dynamischen Transfersteifigkeit führt. Hier- 
bei kann festgestellt werden, dass sich mit zunehmender Bohrungsanzahl im 
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Abbildung 3.12: Darstellung der variierten hydraulischen Parameter des Düsen- 
Membran-Systems 


Anzahl Bohrungen Dicke Querschnitts- 


tin [mm] fläche in [mm?] 


Variante 1 12 6,5 398 
Variante 2 21 6,5 697 
Variante 3 30 6,5 995 
Variante 4 36 6,5 1195 
Variante 5 12 3.0 398 
Variante 6 21 3.0 697 
Variante 7 30 3.0 995 
Variante 8 36 3.0 1195 


Tabelle 3.3: Geometrische Eigenschaften der untersuchten Hydraulikvarianten 
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Abbildung 3.13: Dynamische Transfersteifigkeit der untersuchten Hydraulikvarian- 
ten 1-4 
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Abbildung 3.14: Dynamische Transfersteifigkeit der untersuchten Hydraulikvarian- 
ten 5-8 
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Düsen-Membran-System die Resonanz zu höheren Frequenzen verschiebt. 
Dies bedingt eine höhere Frequenzlage des Maximums in der Phase und ei- 
ne Verschiebung des Anstiegs sowie des Maximums des Betrages der dyna- 
mischen Transfersteifigkeit. Bei den Varianten mit zwölf Bohrungen (Va- 
riante 1 und 5) liegt der Betrag der dynamischen Steifigkeit am unteren 
Ende des betrachteten Frequenzbereiches deutlich höher, als bei den Va- 
rianten mit höherer Bohrungsanzahl. Dies kann auf die Änderung des hy- 
draulischen Übersetzungsverhältnisses zwischen der Offenfläche des Ent- 
koppelgitters und der Tragfederfläche und die damit verbundene Änderung 
des über die Membran entkoppelten Schwingweges zurückgeführt werden. 
Auch die Abnahme des Betrages unmittelbar vor der Resonanzfrequenz bei 
diesen beiden Varianten kann über diese Reduktion des wirkenden Entkop- 
pelspieles begründet werden. Aufgrund des gewählten Prüfablaufes mit ei- 
ner konstanten Schwingbeschleunigungsamplitude von 30 m/s? kommt es 
zu einer quadratischen Abnahme des Schwingweges über der Frequenz. Bei 
den Varianten mit zwölf Bohrungen kann dadurch bei tiefen Frequenzen die 
Volumenänderung in der Arbeitskammer nicht mehr vollständig durch eine 
Bewegung der Membran aufgenommen werden und es resultiert eine höhe- 
re Steifigkeit. Nimmt der Schwingweg nun über der Frequenz ab, sinkt die 
dynamische Steifigkeit ab, bis das eingestellte Entkoppelspiel unterschritten 
wird. 

Beim Vergleich der Varianten mit identischem Lochbild, aber unterschiedli- 
cher Gitterhöhe (jeweils die Paarungen Variante | und 5; 2 und 6; 3 und 7; 4 
und 8) kann festgestellt werden, dass die Resonanzfrequenz bei einer gerin- 
geren Gitterhöhe höher liegt als bei einer größeren Dicke des Gitters. Eben- 
so zeigt sich, dass bei einer höheren Gitterhöhe die Maximalwerte der Phase 
höher liegen als bei der entsprechenden Variante mit geringerer Gitterhöhe. 
Zusammenfassend kann geschlussfolgert werden, dass mit einer zunehmen- 
den Gitterhöhe die Frequenzlage des Resonanzeffektes absinkt und dass mit 
einer Zunehmenden Anzahl von Bohrungen die Frequenzlage ansteigt. 

Für eine Auslegung von hydraulischen Aggregatelagern ist es sinnvoll, die 
Frequenzlage dieses Resonanzeffektes so einzustellen, dass die für den je- 
weiligen Antrieb dominante Motorordnung bei Maximaldrehzahl noch un- 
terhalb dieses Anstieges im Betrag der dynamischen Steifigkeit liegt. 
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4 Modellbildung des 


Ubertragungsverhaltens von 
Aggregatelagern 


Die experimentellen Analysen des vorherigen Kapitels zeigen, dass die 
dynamischen Transfersteifigkeiten von Buchsenlagern und hydraulischen 
Topflagern vor allem von den Materialeigenschaften des verwendeten Ela- 
stomerwerkstoffes dominiert werden und sich im betrachteten Frequenzbe- 
reich Resonanzeffekte ausbilden. Im nachfolgenden Kapitel soll auf Basis 
der herausgearbeiteten Einflussparameter (Werkstoff, Temperatur, hydrauli- 
sche Abstimmung bei hydraulischen Lagern) eine Modellbildung des Bau- 
teilverhaltens betrieben werden. 


4.1 Grundlagen der Modellierungsstrategien 


Um die Modellierung des Bauteilverhaltens in diesem Kapitel durchfüh- 
ren zu können, sollen im folgenden Abschnitt zunächst die Grundlagen der 
Finiten-Elemente-Methode dargelegt werden und anschließend noch auf die 
Bedatung des in Abschnitt 2.2.2 vorgestellten Materialmodells eingegangen 
werden. 


4.1.1 FE-Modellierung 


Um eine Modellierung des Systemverhaltens von Aggregatelagerelementen 
durchführen zu können, welche Kontinuumsschwingungen der Einzelkom- 
ponenten abbilden kann, ist eine Beschreibung mit Hilfe von finiten Ele- 
menten (FE) besonders geeignet. Die theoretischen Grundlagen, welche die 
Basis für kommerziellen FE-Programmen bilden, sollen im folgenden Ab- 
schnitt kurz erläutert werden. 

Grundlage für die Herleitung ist die Gleichgewichtsbedingung der Elasto- 
dynamik [FB09, Kle12]: 


Vio+b=-pü auf Q (4.1) 
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mit: Vs dem symmetrischen Gradientenoperator; © dem Spannungsvektor; 
b dem Vektor der Volumenlasten; ti der Beschleunigungsvektor; Q das 
gesamte betrachtete Volumen 


Oxx 

oO 

o ð o yy 

UL ER E 

> a a ð — zz 
Vi=lo 2 6 2.0.2 =| (4.2) 

002 90 2 2 “ 

Oz Ox ð 
x y Or: 
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Zusatzlich gilt der Zusammenhang zwischen der Verschiebung und der Deh- 
nung [HW10]: 


Ey 2 0 0 
E 0 2 0 
en i 
E 0 0 FE 
e= Vu=|, 3 7 Uy (4.3) 
Vey ty x 0 a 
d d z 
ez 0 & d 
ð d 
e a l & 


mit: € dem Dehnungsvektor 


Fiir die Materialeigenschaften gilt das lineare Hooke’sche Gesetz im dreidi- 
mensionalen Fall [HW10, FB09]: 


o=Ee (4.4) 


mit: E der materialabhängigen, symmetrischen Elastizizätsmatrix 
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4.1 Grundlagen der Modellierungsstrategien 


Für den Fall eines isotropen Materialverhaltens hat die Matrix E die Form: 


1=V. v v 0 0 0 

v 1-v v 0 0 0 

E v v l-v 0 0 0 

~ (+v)(1=2¥) | o 0 o LY o oO 
0 0 0 o + o0 

0 0 0 0 Oo > 


(4.5) 
Zusätzlich zur Gleichgewichtsbedingung müssen die Randbedingungen auf 
dem Rand des betrachteten Bereiches erfüllt werden: 


Spannungsrandbedingung: tn=t auf T, (4.6) 
Verschiebungsrandbedingung: u=ü auf T, (4.7) 


mit: T dem dreidimensionalen Spannungstensor; n dem Normalenvektor; t 

dem Spannungsvektor am Neumann-Rand; u dem Verschiebungsvektor; ü 

der Verschiebungsvektor auf dem Dirichlet-Rand; I, dem Neumann-Rand; 
T„ dem Dirichlet-Rand 


Nach dem Prinzip der virtuellen Verschiebungen [Kle12, HW10] gilt für 
jeden Spannungszustand, dass die virtuelle Arbeit der äußeren Kräfte 6W, 
gleich der virtuellen Arbeit der inneren Kräfte 6W, ist: 


OW, — 6Wı =0 (4.8) 


Dieses Prinzip kann nach [Kle12] und mit Gleichung 4.1 geschrieben wer- 
den zu: 


| be! odQ— | bu! bdQ— f ŝul tdr, =0 (4.9) 
Q Q I 


Aus dem Prinzip der Virtuellen Verschiebungen lässt sich mit Hilfe von 
Gleichung 4.3 und 4.4 das elastische Potential x des Systems ableiten zu: 


T= J (Vsu)"E(Vsu)d2- | uodo- f u tdr; (4.10) 
2Jo Q I; 
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Betrachtet man zusätzlich wirkende Einzelkräfte auf das System, so ergibt 
sich das elastische Potential zu [Kle12]: 


1 = nF = 
T= 5 (Vsu)"E(Vsu)a2-— | u'ban - f utd; — Yul Fr 

2 JQ Q Ñ fo 

(4.11) 


mit: Fg dem Kraftvektor am Ort k; ug der Verschiebung am Ort k;Fx der 
Anzahl der Punkte, an denen äußere Kräfte angreifen 


Die finite Elemente Methode beruht darauf, das Gesamtsystem in Teilele- 
mente zu zerlegen und innerhalb eines Elementes für das Verschiebungsfeld 
eine Ansatzfunktion zu verwenden. Je nach verwendetem Elementtyp kann 
diese Ansatzfunktion variieren. In den meisten Fällen werden jedoch Poly- 
nome ersten oder zweiten Grades verwendet. Das Verschiebungsfeld ergibt 
sich mit diesen Ansatzfunktionen zu: 


uxu? = N°d‘ (4.12) 


mit: uf dem Verschiebungsfeld im Element e; N° der Ansatzfunktion im 
Element e; d° dem Knotenverschiebungsvektor im Element e 


Die Dehnungen im Element ergeben sich demnach zu: 
exe’ = Vsu’ = VsN°d’ = Bed’ (4.13) 
mit: B° der Dehnungs-Verschiebungs-Matrix im Element e 


Das elastische Potential für ein Element kann durch einsetzen von Glei- 
chung 4.13 und 4.12 in Gleichung 4.11 ermittelt werden zu: 


TE = ; | dT BT EC Be d dQ! u f dT NeTpedo® 
Qe Qe 
né (4.14) 
= $ A7 NTUN- Y d" N° (4) Fe 
t k=1 


Anhand von Gleichung 4.14 werden folgende, die Notation vereinfachende 
Matrizen eingeführt [FB09]: 
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e Elementsteifigkeitsmatrix: 


gem | BT E B°dQ (4.15) 
Qe 


e Kraftvektor der Volumenlasten: 


p= f NTa (4.16) 
Qe 


e Kraftvektor der Oberflächenlasten: 


= f ; N* far? (4.17) 


e Kraftvektor der Einzelkräfte: 


Np _ 
fe = } N° (xi) Fx (4.18) 
k=1 
e Elementkraftvektor 
F= f+ fé+ fi (4.19) 


Durch ausklammern der Knotenverschiebungsvektoren aus den einzelnen 
Termen und Einsetzen von Gleichungen 4.15-4.19 in Gleichung 4.14 ergibt 
sich für das elastische Potential: 


1 
n = 5d! Kd" ay (4.20) 
Zur Lösung von dynamischen Problemen ist zusätzlich zum elastischen Po- 


tential eines Elements die kinetische Energie T° des jeweiligen Elements zu 
bestimmen. Diese ergibt sich nach Schwarz [Sch84] zu: 


1 
Te =5 f pT ù do (4.21) 
Qe 


Hierbei ist [FB09]: p 
ù= N‘d(t) (4.22) 


Führt man in Gleichung 4.21 die Abkürzung 
M= [ PNT Neda’ (4.23) 
Qe 
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mit der Element-Massenmatrix M ein, so ergibt sich die Lagrange-Gleichung 
mit Gleichung 4.21 und Gleichung 4.20 zu: 


i's ae | 
L=T-n= z4” Md" = s Kd -df (4.24) 


Die Lagrangesche Bewegungsgleichung 


d f 6L ôL 
(=) =d (4.25) 


ergibt letztendlich das Differentialgleichungssystem des Elements 
M’d+K‘d-f=0 (4.26) 


Die Differentialgleichungssysteme der einzelnen Elemente können über 
Assemblierungsmatrizen dann zum globalen Differentialgleichungssystem 
zusammengeführt werden. Somit kann das betrachtete Feldproblem in ei- 
ne Differentialgleichungssystem überführt werden, welches anschließend 
durch einen entsprechenden Algorithmus gelöst werden kann. 


4.1.2 Bedatung der Materialmodelle 


Das Körperschallübertragungsverhalten von Aggregatelager hängt, wie im 
vorherigen Kapitel gezeigt, vor allem von der Elastomerkomponente ab. 
Ein wichtiger Schritt bei der Modellbildung des Bauteilverhaltens ist da- 
her die Abbildung des Werkstoffverhaltens des jeweiligen Elastomers. Die 
Modellierung des Werkstoffes soll im Folgenden erläutert werden. Zur Mo- 
dellierung des Werkstoffverhaltens der für Aggregatelager verwendeten 
Elastomere werden zunächst die Transfersteifigkeiten von Bauteilen un- 
terschiedlicher Naturkautschukvulkanisate (NR-Vulkanisate) experimentell 
ermittelt. Anschließend soll das Verhalten dieser Bauteile mit Hilfe einer 
FE-Simulation nachgestellt werden. 

Über diese numerische Analyse dieser Bauteile und ein Least-Square-Fitting- 
Verfahren wird ein isotropes, viskoelastisches Materialmodell bedatet, so- 
dass die Simulation das Verhalten der experimentell untersuchten Bauteile 
einschließlich der entsprechenden Eigenmoden abbildet. Hierbei werden 
sowohl das Speicher- als auch das Verlustmodul sowie die Querkontrakti- 
onszahl für vier Elastomervulkanisate angepasst. Diese vier verschiedenen 
Vulkanisate stellen eine üblicherweise für Aggregatelagerelemente verwen- 
deten Bereich von Naturkautschukmischungen dar. 
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Abbildung 4.1: durch Fitting ermitteltes komplexes Elastizitätsmodul verschiedener 
NR-Mischungen über der Frequenz 


Auf Basis der in Abschnitt 2.2 dargelegten Werkstoffeigenschaften von Ela- 
stomeren, werden für das Least-Square-Verfahren folgende Annahmen ge- 
troffen: 


e der Betrag des komplexen Elastizitätsmoduls nimmt über der Fre- 
quenz linear zu 


e die Querkontrakzionszahl ist im betrachteten Frequenzbereich mono- 
ton fallend 


Die sich aus dem Fitting-Verfahren ergebenden komplexen E-Moduli sind 
in Abbildung 4.1 dargestellt. Abbildung 4.2 zeigt vergleichend die expe- 
rimentell ermittelten, sowie die numerisch mithilfe des Fitting-Verfahrens 
approximierten komplexe Transfersteifigkeit verschiedener Mischungsvari- 
anten des für das Fitting verwendeten Getriebelagers. Es ist zu erkennen, 
dass der Betrag der Transfersteifigkeit mithilfe der numerischen Analyse 
nachgebildet werden kann und die Übereinstimmungen zwischen experi- 
menteller und numerische Analyse insbesondere im Frequenzbereich bis 
zwei Kilohertz hoch sind. Das Verhalten, welches nicht durch Eigenmoden 
der Bauteile bestimmt ist, wird dabei am genauesten abgebildet. Im Bereich 
der Eigenfrequenzen wird die Frequenzlage bis auf wenige Prozent Abwei- 
chung sowohl in der Frequenzlage als auch in der Überhöhung durch die 
Resonanzen hinreichend genau abgebildet. Die Abweichungen zwischen 
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Abbildung 4.2: Vergleich der experimentell und numerisch ermittelten Transferstei- 
figkeiten eines Getriebelagers mit verschiedenen Naturkautschuk- 
mischungen 
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numerischer und experimenteller Analyse sind hierbei geringer als die Va- 
rianz zwischen verschiedenen Chargen der Elastomerkomponenten. 

Der Phasenverlauf der numerisch ermittelten Werte entspricht bis zur ersten 
Eigenfrequenz sehr genau dem jeweiligen experimentell ermittelten. Ober- 
halb der ersten Eigenfrequenz der untersuchten Varianten wird die Differenz 
der numerischen und experimentell ermittelten Transfersteifigkeit für jede 
der untersuchten Naturkautschukmischungen mit zunehmender Frequenz 
immer größer. Für Frequenzen oberhalb von zwei Kilohertz ergeben sich 
sowohl für den Betrag als auch für die Phase größere Abweichungen, ins- 
besondere für die Varianten mit Härten von 55 Sh und 60 Sh. Die Abwei- 
chungen zwischen Simulation und Messung nehmen jedoch bei allen unter- 
suchten Varianten oberhalb von 2000 Hz deutlich zu. Diese Abweichungen 
lassen sich unter Berücksichtigung der genannten Randbedingungen für das 
komplexe Modul und die Querkontraktionszahl auch nicht weiter reduzie- 
ren. Die Ursache für diese Abweichungen liegen daher vermutlich in einer 
Abweichung von der Annahme eines isotropen, homogenen Materials. 

Zur Untersuchung der Homogenität der Materialien soll die Streuung der 
Shore-Härte innerhalb eines Bauteils ermittelt werden. Hierzu werden Schnit- 
te angefertigt und das Bauteil entlang dieser Schnittflächen in einem Raster 
von 3mm mit einem Durometer abgetastet. Hieraus ergeben sich 19 Mes- 
spunkte über der Querschnittfläche des untersuchten Bauteils. 

Abbildung 4.3 zeigt die statistische Verteilung der Shorehärte über diese 
Messpunkte für die vier untersuchten NR-Mischungen in einem Boxplot. 
Es ist zu erkennen, dass die Messwerte im Allgemeinen oberhalb des Nenn- 
wertes liegen. Des weiteren kann Abbildung 4.3 entnommen werden, dass 
die Streuung zu höheren Nenn-Shore-Härten hin zunimmt. Dies drückt sich 
weniger stark in den minimalen bzw. maximalen Werten der jeweiligen Rei- 
hen aus, vielmehr jedoch in den 25 % und 75 % Perzentilen. 

Da bei steigender Frequenz die Wellenlänge der Schwingphänomene immer 
kleiner wird, gewinnen lokale Steifigkeitsverteilungen im Bauteil zuneh- 
mend an Bedeutung für das Gesamtsystemverhalten. Diese lokalen Effekte 
verändern das Schwingverhalten des Systems gegenüber einer numerischen 
homogenen Betrachtung. Durch diese Abweichung der Steifigkeitsvertei- 
lung vom Mittelwert in Gummibauteilen kommt es mit zunehmender Fre- 
quenz zu einer Abweichung zwischen experimentell und numerisch ermit- 
telter Transfersteifigkeit. Diese Inhomogenitäten und deren Ursache soll im 
Zuge einer Validierung der Bauteilsimulation in Abschnitt 4.3.1 noch wei- 
tergehend untersucht werden. 
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Abbildung 4.3: Darstellung der Streuung der Shore(A)-Härte innerhalb eines Bau- 
teils 


4.2 Modellierung des Bauteilverhaltens von 
Buchsenlagern 


Die im vorherigen Abschnitt parametrierten Materialmodelle sollen anhand 
der in Abschnitt 3.3.1 untersuchten Varianten von 4-Steg-Buchsenlagern 
validiert werden und dabei die Eignung zur Vorhersage der dynamischen 
Transfersteifigkeit von Elastomeren mit Hilfe von FE-Modellen überprüft 
werden. 


4.2.1 Modellbildung und Simulationsrandbedingungen 


Zur Ermittlung der dynamischen Transfersteifigkeit mit Hilfe eines FE- 
Modells wird analog zum Vorgehen am Prüfstand eine Verschiebung auf 
den Kern in die Z-Richtung des Buchsenlagers aufgebracht und am einge- 
spannten Außenring die Reaktionskraft ermittelt. Aus dem Verhältnis der 
komplexen Reaktionskraft und der aufgebrachten Verschiebung lässt sich 
anschließend die dynamische Transfersteifigkeit ermitteln. Zur Modellie- 
rung der Elastomerkomponente wird das im vorherigen Abschnitt erläuter- 
te viskoelastische Materialmodell verwendet. Abbildung 4.4 zeigt das FE- 
Modell des untersuchten Buchsenlagers sowie das zugehörige Koordinaten- 
system und stellt die Anbindungsstellen des Außenrings und des Kerns dar. 
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Die betrachteten Bauteile weisen eine gänzliche andere Geometrie auf, als 
diejenigen, welche für die Bedatung der Materialmodelle verwendet wur- 
den. Durch dieses Vorgehen soll gezeigt werden, dass eine Vorhersage von 
Bauteileigenschaften rein auf Basis der Geometrie und unter Verwendung 
der im vorherigen Abschnitt erarbeiteten Materialmodell erfolgen kann. 


Abbildung 4.4: Darstellung des FE-Modells des untersuchten 4-Steg- 
Buchsenlagers 


4.2.2 Simulationsergebnisse und -validierung 


Die mit dem im vorherigen Abschnitt erläuterten Modell erzeugten kom- 
plexen Transfersteifigkeiten sind in Abbildung 4.5 vergleichend mit den 
experimentell erzeugten Daten dargestellt. Es lässt sich erkennen, dass die 
numerisch erzeugten Daten eine hohe Übereinstimmung mit den experi- 
mentellen Daten haben. Sowohl der Verlauf des Betrages und der Phase 
als auch die Absolutwerte der dynamischen Transfersteifigkeit zeigen sich 
sehr ähnlich. Das numerische Modell liefert im Minimum der Phase noch 
geringere Werte als die experimentellen Untersuchungen. Ursächlich hier- 
für ist vermutlich die im Vergleich zu den zur Parametrierung verwendeten 
Materialien geringfügig andere Werkstoffdämpfung. Die Materialdämpfung 
kann ohne eine signifikante Shore-Härten-Änderung des Werkstoffes vari- 
ieren, sodass in den experimentell untersuchten Buchsenlagern eine etwas 
höhere Dämpfung vorhanden sein könnte. Diese Materialdämpfung hätte im 
Speziellen in den Resonanzen eine signifikante Auswirkung auf das Bau- 
teilverhalten, vor allem auf die Phase der dynamischen Transfersteifigkeit. 
Hierdurch weicht der Phasenwinkel zwischen Anregung und Kraftantwort 
in den Resonanzeffekten zwischen Versuch und Simulation leicht voneinan- 
der ab. Die untersuchte Variante mit 55 Sh zeigt einen leicht höheren Betrag 
der dynamischen Transfersteifigkeit in den experimentellen Daten und der 
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Abbildung 4.5: Vergleich der experimentell und numerisch ermittelten Transferstei- 
figkeit von Buchsenlagern 
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Abfall der Phase erfolgt bei höheren Frequenzen als bei der numerischen 
Betrachtung dieser Variante. Es lässt sich daher schlussfolgern, dass das 
Vulkanisat bei dieser Variante ein etwas höheres Elastizitätsmodul als das 
zur Parametrierung des Materialmodells verwendete Pendant. Diese Ab- 
weichungen sind jedoch im Rahmen der Fertigungstoleranzen akzeptabel, 
da aufgrund der Herstellung und der Verwendung eines Naturkproduktes als 
Rohstoff für die Vulkanisate leichte Schwankungen unvermeidbar sind. 


Die numerische Analyse der Bauteile erlaubt es die Schwingform derjenigen 
Betriebsmode der Buchsenlager zu analysieren, welche die Transfersteifig- 
keit beeinflusst. Die Schwingform ist in ihren beiden maximalen Auslen- 
kungen (bei Phasenwinkel +7/2) in Abbildung 4.6 dargestellt Der Abbil- 


“Lx 


Abbildung 4.6: Darstellung der die Transfersteifigkeit von Buchsenlagern beein- 
flussenden Mode 


dung kann entnommen werden, dass es sich bei der Betriebsmode um eine 
transversale Biegeschwingung der einzelnen Stege handelt. Die Schwin- 
gungen aller vier Stege sind dabei in die Belastungsrichtung gleichphasig. 
Diese Biegeschwingung tritt für alle vier Stege bei der gleichen Frequenz 
auf, wodurch es zu einer starken Überhöhung der Reaktionskräfte kommt. 
Oberhalb dieser Resonanzfrequenz kommt es zu einer dynamischen Ent- 
kopplung der Gummistege, wodurch die Transfersteifigkeit deutlich absinkt 
und geringer ist als unterhalb der Biegemode. Es kommt hierbei zu einer 
Isolationswirkung der Gummistege durch eine überkritische Frequenzlage. 
Dieser Effekt scheint für die Optimierung der Übertragungseigenschaften 
günstig und sollte dahingehend ausgelegt werden, dass Frequenzbereiche, 
in denen die Karosserieeigenschaften Empfindlichkeiten aufweisen, ober- 
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halb dieser Resonanzfrequenz liegen. 


4.3 Modellierung des Bauteilverhaltens von hydraulisch 
dämpfenden Lagerelementen 


Nachdem in vorherigen Abschnitt die Modellierung von Elastomer-Buch- 
senlagern erfolgt ist, soll im Folgenden die Modellierung von hydraulisch 
gedämpften Topflagern vorgestellt werden. Hierbei soll zunächst analog zur 
Modellierung der Buchsenlager gezeigt werden, dass sich die Elastomer- 
komponenten mit dem verwendeten Materialmodell im FE-Modell abbilden 
lassen. Anschließend sollen die Ergebnisse der Parameterstudie zum Ver- 
halten der Hydraulikkomponenten genutzt werden, um eine Modellbildung 
der Hydraulikeigenschaften durchzuführen. Die Ergebnisse dieser beiden 
Teilmodelle sollen daraufhin kombiniert werden, um das Gesamtsystemver- 
halten abbilden zu können. 


4.3.1 Modellierung der Elastomerbauteile 


Die Tragfeder eines hydraulischen Aggregatlagers besteht, wie in Abschnitt 
2.3.3 vorgestellt, im wesentlichen aus drei Bestandteilen. Diese drei Be- 
standteile sind: ein metallischer Kern zur Befestigung des Bauteils an ei- 
ner Motorkonsole, ein daran anvulkanisierter, lasttragender Elastomeranteil 
und einem diesen Elastomeranteil umschließenden Außenring, welcher zur 
Fixierung der Tragfeder im Gehäuse des Bauteils dient. Im folgenden Ab- 
schnitt soll erläutert werden, wie diese Tragfeder durch eine numerische Si- 
mulation abgebildet werden kann. Die Ergebnisse sollen anschließend durch 
einen Abgleich mit den experimentellen Daten validiert werden. 


Modellbildung und Simulationsrandbedingungen 


Um die Transfersteifigkeit der Tragfeder zu bestimmen, wird nur der Elasto- 
meranteil in der Simulation abgebildet, da sowohl Kern als auch Außenring 
nur zur Anbindung an die sich jeweils anschließende Fahrzeugstrukur be- 
nötigt werden und das Übertragungsverhalten nur in sehr geringem Umfang 
beeinflussen. Der Kern befindet sich sowohl bei den Untersuchungen auf 
dem Prüfstand als auch in den Simulationen auf der angeregten Seite. Die- 
ser stellt somit eine Masse dar, welche beschleunigt werden muss und somit 
eine Kraft hervorruft. Am Prüfstand wird diese zusätzliche Masse durch 
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eine entsprechend größere Kraft durch den elektrodynamischen Shaker be- 
schleunigt. In der Simulation würde dadurch die Kraft steigen, welche zur 
Anregung des Systems notwendig ist. Die Kraft, die auf der passiven Sei- 
te übertragen wird, wird dadurch jedoch nicht nicht beeinflusst. Durch die 
Nichtbetrachtung des Kerns in der Simulation lässt sich die Anzahl der zu 
berücksichtigenden Elemente und Freiheitsgrade reduzieren, was zu einer 
Reduktion von Rechenzeiten führt. Abbildung 4.7 zeigt einen Schnitt durch 
eine solche Tragfeder ohne Kern und Außenring und stellt die Anbindungs- 
stellen dar, an denen diese Komponenten anvulkanisiert sind. 

Zur Ermittlung der Transfersteifigkeit wird der Anbindungsbereich des 


Anbindung Kern 


Anbindung Außenring 
“Lx 
Abbildung 4.7: Schnittdarstellung durch eine Tragfeder eines hydraulisch dämpfen- 
den Aggregatlagers 


Kerns mit einer vorgegebenen Verschiebung in z-Richtung beaufschlagt. 
Der Anbindungsbereich des Außenrings wird fest eingespannt und dort 
die Reaktionskraft ermittelt. Durch dieses Vorgehen wird in einem ersten 
Schritt die Motorvorlast auf die Tragfeder aufgebracht und in einem zwei- 
ten Schritt eine sinusförmige Verschiebungen mit Frequenzen zwischen 100 
und 2500Hz aufgebracht. Aus der Verschiebung und der komplexen Re- 
aktionskraft kann anschließend die Transfersteifigkeit des Elastomeranteils 
ermittelt werden. 

Für die Modellierung des Elastomers soll dasselbe Materialmodell verwen- 
det werden, welches auch für die numerische Simulation der Buchsenlager 
verwendet wird. Die Elastomermischungen der untersuchten Tragfedern 
stammen dabei nicht vom selben Rohstofflieferanten, wie die der in Ab- 
schnitt 4.2 betrachteten Buchsenlagern. Es lässt sich daher ebenfalls die 
Übertragbarkeit der verwendeten Materialmodelle über verschiedene Naur- 
kautschukmischungen validieren. 
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Simulationsergebnisse und -validierung 


Die Ergebnisse der Ermittlung der Transfersteifigkeiten nach der im vor- 
herigen Abschnitt beschriebenen numerischen Methode ist in Abbildung 
4.8 dargestellt. Im Vergleich zu den numerisch ermittelten Ergebnissen sind 
die experimentell ermittelten Transfersteifigkeiten von Tragfedern eines hy- 
draulischen Motorlagers abgebildet. Diese experimentellen Daten sind nach 
der in Abschnitt 3.2 beschriebenen Methodik an Tragfedern ohne Gehäuse 
und Hydraulik, jedoch mit Kern und Außenring, ermittelt. 

Aus Abbildung 4.8 kann entnommen werden, dass auch die Elastomer- 

komponenten hydraulischer Topflager durch das in Abschnitt 4.1.2 bedatete 
Materialmodell beschrieben werden können. Sowohl der Verlauf des Be- 
trages als auch der Verlauf der Phase der dynamischen Transfersteifigkeit 
werden bis etwa 2000 Hz in Bezug auf die Lage von Maxima und auch der 
Absolutwerte sehr genau nachgebildet. Im betrachteten Frequenzbereich 
zeigen sich bei allen untersuchten Varianten zwei charakteristische Reso- 
nanzeffekte, welche mit zunehmender Shore-Härte zu höheren Frequenzen 
verschoben werden. Das untersuchte Bauteil mit 45 Sh zeigt eine zusätzliche 
dritte Resonanz, welche bei den übrigen aufgrund der erhöhten Frequenz- 
lagen der Bauteilmoden oberhalb des Messbereichs liegt. Ab etwa 2000 Hz 
zeigt sich jedoch, wie auch schon bei der Materialmodellbedatung und der 
numerischen Analyse der Buchsenlager, dass die experimentellen und die 
numerischen Werte stärker voneinander abweichen. 
In der numerischen Analyse zeigen sich die Bauteilmoden, welche die 
Transfersteifigkeit der Elastomeranteile stark beeinflussen, als Transversal- 
schwingung entlang der Mantelhöhe auf der rotationssymmetrischen Gum- 
misehne der Tragfeder. Die Schwingformen, die zu diesen Moden gehören, 
sind in Abbildung 4.9 dargestellt. Hierbei ist jeweils die maximale bzw. 
minimale Auslenkung dargestellt, welche jeweils einen Phasenabstand von 
180° zueinander bzw. +90° zur aufgeprägten Anregung aufweisen. 


Bei der Betrachtung der Schwingformen stellt sich auf der Länge der Man- 
telhöhe (&-Koordinate des in Abbildung 4.7 dargestellten Koordinatensy- 
stems) für die erste Mode eine Halbwelle in transversale Richtung, bei der 
zweiten Mode eine Vollwelle mit einem Schwingungsknoten auf der Mitte 
der Mantelhöhe der Tragfeder ein. Die dritte Mode weist nochmals einen 
Schwingungsknoten mehr auf und zeigt somit eine Schwingform, welche 
drei Halbwellen in die transversale Richtung über der Mantelhöhe der Trag- 
feder entspricht. 
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Abbildung 4.8: Vergleich der experimentell und numerisch ermittelten Transferstei- 
figkeit der Tragkörper des untersuchten Hydrolagers 
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(a): Mode 1 


(b): Mode 2 


(c): Mode 3 


Abbildung 4.9: Darstellung der Knotenverschiebungen bei den die Transfersteifig- 
keit von Tragfedern beeinflussenden Moden 


Validierung der Ergebnisse durch Laser-Doppler-Vibrometrie 
Durch eine Untersuchung der Bauteile mit Hilfe eines Laser-Doppler- 
Vibrometers (LDV) in dem bereits fiir die experimentellen Analyse ver- 
wendeten Versuchsaufbaus können diese Betriebsschwingformen ebenfalls 
experimentell ermittelt werden!. Hierzu wird mittels eines eindimensiona- 
len LDVs an einem Ausschnitt der Mantelfläche der Tragfeder während der 
experimentellen Analyse die Oberflächenschnelle ermittelt. Abbildung 4.10 
stellt diesen Messbereich und das dabei verwendete Koordinatensystem dar. 
Der Messbereich erstreckt sich über den gesamten Bereich 0 <hy< hm max 
und —25° <@y< 25°. Eine Begrenzung des Messbereichs ist aufgrund der 
rotationssymmetrischen Bauweise und der Betrachtung mit Hilfe eines ein- 
dimensionalen LDVs notwendig. Der Messbereich wird dabei so gewählt, 


! Teile dieser Untersuchungen sind bereits in einer vom Autor der vorliegenden Arbeit betreu- 
ten Abschlussarbeit behandelt worden, vgl. [Tre17] 
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dass der durch den angularen Versatz entstehende Messfehler auf maximal 
10 % begrenzt wird. 


Mit der LDV-Untersuchung lässt sich die Oberflächenschnelle des Tragkör- 
pers bei den einzelnen Bauteilmoden auf der Außenseite des Tragkörpers 
ermitteln. Hierdurch können die numerisch ermittelten Schwingformen qua- 
litativ mit den experimentell ermittelten verglichen werden. Abbildung 4.11 
zeigt diese Oberflächenschnellen für die untersuchte Variante mit 45 Sh 
für den in Abbildung 4.10 dargestellten Ausschnitt in einer Seitenansicht, 
sodass jeweils der Ausschnitt 0° < @y < 25° abgebildet ist. Die Schwing- 
formen unterscheiden sich im Wesentlichen nicht zwischen den einzelnen 
Elastomermischungen, wobei die Frequenzlage der Moden, wie bereits zu- 
vor diskutiert, mit zunehmender Shore-Härte ansteigt. 

Der Darstellung der Oberflächenschnelle in Abbildung 4.11 kann entnom- 
men werden, dass die erste Schwingform einer transversalen Schwingung 
der Tragfedersehne mit einer halben Wellenlänge entspricht. Die Oberflä- 
chenschnelle weist ihr Maximum bei @y = 0 auf, was durch die Analyse 
mit einem eindimensionalen LDV begründet ist. Mit zunehmendem Ab- 
stand nimmt aufgrund des sich einstellenden Winkelfehlers die Schnelle ab. 
Die zweite und dritte Eigenschwingform weisen ein zur numerischen Ana- 
lyse analoges Verhalten auf, da die Anzahl der Schwingungsknoten über der 
Sehnehnlänge hy mit jeder Mode um einen Knoten ansteigt. Hierdurch zeigt 


Abbildung 4.10: Messbereich des Tragfedermantels der LDV-Untersuchung 
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Abbildung 4.11: Darstellung der experimentell ermittelten Schwingformen der 
Tragfedern bei den Eigenfrequenzen bei 90° (links) und —90° 
Phasendifferenz zur Anregung [Tre17] 


sich das bereits zuvor diskutierte und in Abbildung 4.9 auf Basis der nume- 
rischen Analyse dargestellte Verhalten. Die zweite Eigenschwingform weist 
eine ganze Transversalwelle über der Sehnenhöhe auf, die dritte Mode drei 
Halbwellen. Bei der Skalierung der Oberflächenschnelle ist zu beachten, 
dass durch die verwendete Anregung einer über der Frequenz konstanten 
Beschleunigungsamplitude die Schnelle linear über der Frequenz abnimmt. 
Dies hat zur Folge, dass höhere Moden eine geringere Oberflächenschnelle 
aufweisen als die niederen und die Skalierung der Graphiken entsprechend 
unterschiedlich ausfallen. 

Bei Betrachtung der Oberflächenschnellenverteilung bei der zweiten und 
dritten Mode kann festgestellt werden, dass das Maximum bzw. Minimum 
der Schnelle nicht wie erwartet bei @y = 0° auftritt, sondern bei der zweiten 
Mode im oberen Schwingungsbauch bei @y ~ 10° und bei der dritten Ei- 
genmode bei Pu ~ 18°. Diese Abweichung der Schnelleverteilung von der 
zu erwartenden deutet widerum darauf hin, dass es aufgrund von Inhomo- 
genitäten im Elastomer zu einer nicht rotationssymmetrischen Schwingung 
der Tragfeder kommt. 
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Untersuchung von Materialinhomogenitäten 

Um die Einflüsse der Materialinhomogenität weiter zu untersuchen, soll die 
in Abbildung 4.11 beobachtete Asymmetrie der transversalen Schwingung 
anhand von Durometeranalysen nach DIN 7619-1 [DINO2] zur Härtevertei- 
lung in den einzelnen Bauteilen erfolgen. Hierzu soll zunächst der Verlauf 
der Shore-Härte der für die Ermittlung der experimentellen Schwingformen 
verwendeten Tragfeder erfolgen. Eine Auswertung des Härteverlaufs über 
der @y Koordinate auf der Mantelinnen- und Außenseite ist in Abbildung 
4.11 dargestellt. Hierin ist zu erkennen, dass die Außenseite der Tragfeder 
eine höhere Shore-Härte aufweist als die Innenseite. Bei Qy = 6.25° zeigt 
sich ein Minimum in der Shore-Härte, was sich mit der Lage des Maximums 
der Auslenkung bei der zweiten Bauteilmode in Abbildung 4.11 deckt. 
Durch die Ergebnisse der Durometeruntersuchungen können die asymme- 


42 1 L 1 L 1 1 
-25 -18.25-12.5 -6.25 0 625 12.5 18.25 25 
Umfangskoordinate @ in [°] 


e Messpunkte Außen * Messpunkte Innen 
Fit Umfang Außen Fit Umfang Innen 


Abbildung 4.12: Darstellung der Shore-Härte-Streuung über den Umfang der Trag- 
feder [Tre17] 


trischen Schwingungsauslenkungen aus den LDV-Untersuchungen plausi- 
bilisiert werden und die Abweichungen zwischen Bauteilsimulation und ex- 
perimenteller Analyse nachvollzogen werden. Zur weiteren Analyse der Ur- 
sache der Streuung der Werkstoffhärte sollen die zur experimentellen Ana- 
lyse verwendeten Tragfedern bei @y = 90° und bei Pu = —30° entlang des 
Tragkörperquerschnitts in radialer Richtung auf eine Streuung der Shore- 
Härte untersucht werden. Hierzu wird der Querschnitt in 35 äquidistante 
Messpositionen eingeteilt, welche sich auf einem Rechteck mit einer Höhe 
von 10mm und einer Breite von 15mm befinden. Hierbei wird zur Erfüllung 
von DIN ISO 7619-1 (vgl.[DINO2]) ein Abstand zu den freien Rändern und 
zum anvulkanisierten Kern und Außenring eingehalten. Eine Darstellung 
der Messbereiche und die Lage des verwendeten Koordinatensystems ist 
Abbildung A.2 in Anhang A.2 zu entnehmen. Zur Darstellung der Ergebnis- 
se werden die diskreten Messpunkte durch zweidimensionale Splineelemen- 
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te interpoliert. Die Ergebnisse dieser Analyse sind in Abbildung 4.13 für die 
untersuchten Shore-Härten zwischen 45 und 60 Sh dargestellt. Da es sich bei 
den im Rahmen dieser Arbeit untersuchten Bauteilen im Prototypen-Teile 
handelt, wird bezüglich der Inhomogenitäten noch eine weitere Variante mit 
einer nominellen Härte von 43Sh dargestellt (Abbildung 4.13 a), welche ei- 
nem Serienproduktionsprozess entstammt. 

Der Abbildung kann entnommen werden, dass die Härte der untersuch- 


a 49 


47 
45 


10 


= 
a 
nN 
tatari 


y in [mm] 


56 
5z 
542 
53 


y in [mm] 


Abbildung 4.13: Verteilung der Shore-Härte über den Querschnitt der Tragfedern 
bei @y = 90° für Tragfedern der nominellen Shore-Härte: a): 
43Sh (Serienproduktion); b): 45Sh; c): 50Sh; d): 55 Sh; e): 60Sh 
[Tre17] 


ten Bauteile über dem Querschnitt einer starken Streuung unterworfen sind. 
Beim Vergleich der untersuchten Varianten zeigt sich, dass sich diese Streu- 
ung bei den Prototypenteilen in ähnlicher Größenordnung (ca. +2Sh Abwei- 
chung vom nominellen Wert) einstellt, wie bei dem untersuchten Bauteil aus 
dem Serienprozess. Die an den untersuchten Bauteilen abgeleiteten Ergeb- 
nisse können daher auch auf Bauteile aus einer Serienproduktion übertragen 
werden. 

Zunächst kann festgestellt werden, dass die Werte am inneren, freien Rand 
des Tragfederquerschnitts (y = 0mm) bei allen untersuchten Bauteilen ge- 


110 


4.3 Modellierung des Bauteilverhaltens von hydr. Lagerelementen 


ringer sind als auf dem übrigen Tragfederquerschnitt. Diese Abweichung 
ist mit der Prozessführung bei der Vulkanisation zu begründen. Die Vulka- 
nisationswerkzeuge sind mehrteilig aufgebaut, um eine Entformbarkeit des 
fertigen Erzeugnisses sicherstellen zu können. Die Innenseite der Tragfeder 
ist dabei ein separater Teil der Vulkanisationsform, welches für den Vul- 
kanisationsprozess beheizt wird. Bei der Produktion wird nach der Vernet- 
zung zunächst der äußere Teil des Vulkanisationswerkzeuges entfernt und 
das Bauteil anschließend aus der Form entnommen. Dieses Vorgehen hat 
zur Folge, dass die Innenseite des Bauteils einer längeren Temperaturein- 
wirkung ausgesetzt ist und es zu einer Reversion der Vernetzung kommen 
kann (vgl. Abschnitt 2.3.4). Das Bauteil kühlt an der Oberseite (y = 10mm) 
schneller aus und neigt an dieser Stelle weniger zur Reversion. Untersu- 
chungen zur Abweichung der Härte zwischen Innen- und Außenseite sind 
in Abbildung A.3 in Anhang A.2 ausgeführt. Der Kern der Tragfeder ist 
eine massive Aluminiumstruktur und kühlt aufgrund seiner hohen thermi- 
schen Trägheit deutlich langsamer aus als das Elastomer und sorgt somit an 
der Stelle x = Omm für einen weiteren Wärmeeintrag, was an diesem Rand 
ebenfalls für eine geringere Härte sorgt. 

Diesen Reversionseffekten überlagert sind die durch Mischungsaufberei- 
tung und Mastikation erzeugten Materialinhomogenitäten, welche zusätz- 
lich von den Parametern des Einspritzvorganges der Elastomermischung in 
die vorgeheizte Vulkanisationsform beeinflusst werden. Hierdurch kommt 
es zu einer weiteren Streuung der Härteverteilung, welche auch innerhalb 
eines Bauteils für eine nicht rotationssymmetrische Ausbildung der Inho- 
mogenitäten sorgt. Bei Betrachtung von Abbildung 4.13 d) kann dieser Ein- 
fluss beobachtet werden, da es bei x = 0mm nicht zu der Abnahme der Härte 
kommt, welche durch die Reversionseffekte anzunehmen wären. 

Zur weiteren Untersuchung der durch die Mischungsverarbeitung entste- 
henden Inhomogenitäten soll die Korrelation der Härteverteilung Hoo- bei 
Qu = 90° und H_30 bei Pu = —30° herangezogen werden. Nach [Sac04] 
ergibt sich die Korrelation der beiden Vektoren der Dimension 35 x 1 Hoge 
und H,,——30° zu 


cov( Hoo, H_30°) 
/ var (Hoe )var( H_30° ) 


P (Hoe , H_30°) = (4.27) 


Die Ergebnisse dieser Korrelationsanalyse sind fiir den Vergleich der bei- 
den Schnittebenene jeder Shore-Härte in Tabelle 4.1 dargestellt. Unter der 
Annahme der Homogenität der Bauteile wäre eine Korrelationswert von 
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43Sh 45Sh 50Sh 55Sh 60Sh 
P (Hoo, H_30) -0,1567 -0,0569 -0,4555 -0,1464 -0,3374 


Tabelle 4.1: Korrelation der Härteverteilungen in den Elastomermischungen bei 
Qu = 90° und Pu = —30° 


1 zu erwarten. Da die errechneten Korrelationswerte zwischen den beiden 
Schnittebenen negativ sind und zusätzlich zwischen den untersuchten Prüf- 
lingen stark streuen, kann darauf geschlossen werden, dass für die beiden 
untersuchten Querschnitte kein eindeutiger Zusammenhang besteht. Dies 
lässt darauf schließen, dass die durch die Mischungsaufbereitung und den 
Injektionsprozess in die Vulkanisationswerkzeuge erzeugten Materialinho- 
mogenitäten einen größeren Einfluss auf die Streuung der Härte der Bauteile 
haben als die Reversionseffekte. Eine Vorhersage der Verteilung der Shore- 
Härte im Bauteil a priori ist daher nicht möglich. Die numerische Modellie- 
rung der Bauteile kann daher nur von einem homogenen Material ausgehen. 
Dies hat zur Folge, dass lokale Schwingungseffekte, welche bei Frequen- 
zen oberhalb der zweiten Bauteilmode zunehmend an Bedeutung gewinnen, 
nicht vollständig in einer numerischen Simulation abgebildet werden kön- 
nen. Eine numerische Vorhersage des Bauteilverhaltens der Elastomertrag- 
federn von hydraulischen Topflagern kann jedoch bei den betrachteten Bau- 
teilen bis 2 kHz durchgeführt werden. Die Übereinstimmung der Ergebnisse 
zwischen einer experimentellen Analyse und der numerischen Vorhersage 
nimmt jedoch mit zunehmender Lokalisation der Schwingungseffekte und 
somit mit steigender Frequenz ab. 


Validierung der Simulationsergebnisse durch analytische Vergleichs- 
rechnung 

Zur zusätzlichen Validierung der Ergebnisse der Betriebsschwingformana- 
lyse soll im Folgenden eine analytische Betrachtung des Systems durchge- 
führt werden und anschließend die ermittelten Schwingformen der analy- 
tischen, numerischen und experimentellen Ergebnisse verglichen werden. 
Zur analytischen Betrachtung der Tragfedern sollen diese mit Hilfe ei- 
nes quaderförmigen Kirchhoff’schen Plattenelementes beschrieben werden. 
Dieses Vorgehen erweist sich zum Beispiel bei der Modellierung von Ma- 
schinengehäusen und Karosseriebauteilen als dienlich. Als grundlegende 
Annahme wird hier zu Grunde gelegt, dass das betrachtete Element eine 
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Schubstarrheit und eine zu vernachlässigende Drehträgheit aufweist (vgl. 
[GHW07, Irr13]). Zur Beschreibung der orts- und zeitabhängigen Verschie- 
bung w in ¢-Richtung eines sich in &- und n-Richtung (vgl. Koordinatensy- 
stems aus Abbildung 4.7) erstreckenden Plattenelementes ergibt sich nach 
Gross [GHW07] und Irretier [Irr13] die Verschiebungsdifferentialgleichung 
vierter Ordnung im Raum und zweiter Ordnung in der Zeit: 


. atw atw atw 
phw 4 K( Sa Hagen tm) = p(x,y,t), (4.28) 
tK= Eh? 
mS = a2) 


Hierbei ist K die Biegesteifigkeit des Plattenelementes, h die Höhe des Plat- 
tenelementes, p die Dichte des Plattenwerkstoffes, E und v das Elastizitäts- 
modul und die Querkontraktionszahl des Plattenwerkstoffes und p die äuße- 
re Belastung des differentiellen Plattenelementes. Zur Berücksichtigung der 
Einspannungsrandbedingungen der Tragfeder wird die in Abbildung 4.14 
einseitig mit fester Einspannung am Rand AQ(E = 0,7) gelagerte Platte 
der Länge / und Hohe h, welche am gegenüberliegenden Rand 0Q(& =1, n) 
durch ein gelenkiges Loslager gelagert ist. 

Durch die beiden freien Ränder 9Q.(&,n = 0) und 0Q(€,n = b), entfällt 


Abbildung 4.14: zweiseitig gelagerte Rechteckplatte[Tre17] 


die Abhängigkeit der Differentialgleichung 4.28 von der n-Komponente 
und kann durch einen Bernoulli-Produktansatz und die Berücksichtigung 
der Randbedingungen bei 0Q(§ =0,n) und AQ(E =1,n) gelöst werden. 
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Als Gesamtbewegung der Rechteckplatte ergibt sich nach [Irr13] die Super- 
position aller Eigenschwingformen des Plattenelements zu: 


oo 


w(é,n,t) =), Wi(&,n)cos(@;t) (4.29) 


i=1 


Hierbei enstprechen W; den Schwingformen der einzelnen Moden, welche 
sich zur zeitabhängigen Gesamtschwingung des Systems überlagern. Für 
die Schwingformen W; gilt dabei: 


W;(&,n) = a; [(cos(ki§ ) — cosh(k;é ) + (sinh(k;&)—sin(k;&))|, (4.30) 


__ cosh(k;l) — cos(k;l) 
~ sinh(kjl) — sin(kjl) ` 


Hierbei sind K; = k;l die Eigenwerte des zur Ermittlung der Eigenschwing- 
formen notwendigen Eigenwertproblems mit den zugehörigen Eigenkreis- 
frequenzen @;, mit 
K? |K 

Oi = Pon (4.31) 
Bei der Betrachtung der analytischen Lösung der Differentialgleichung kann 
festgestellt werden, dass die Eigenfrequenzen der Lösung nicht den Eigen- 
frequenzen der experimentellen und numerischen Ergebnissen entsprechen. 
Ursächlich hierfür ist die Vereinfachung des Tragfederquerschnitts zu ei- 
nem Rechteck bei der Betrachtung als Kirchhoff-Platte. Eine analytische 
Lösung des Systems kann daher nicht zur Ermittlung des Modalverhaltens 
des Gesamtsystems verwendet werden. Zur Validierung der numerischen 
Ergebnisse kann jedoch trotz dieser Unschärfe ein Vergleich der Schwing- 
formen der ersten Eigenmoden der Tragfeder herangezogen werden. Abbil- 
dung 4.15 stellt diese zur besseren Vergleichbarkeit in normierter Form dar. 
Diese Normierung bezieht sich sowohl auf die &-Richtung als auch auf die 
¢-Richtung des KOS aus Abbildung 4.7. 


Der Darstellung kann entnommen werden, dass sowohl die Schwingformen 
der numerischen Modellierung als auch diejenigen der analytischen Be- 
schreibung durch die Kirchhoff-Platte eine hohe Ubereinstimmung mit den 
experimentell ermittelten Schwingformen haben. Dies gilt insbesondere fiir 
die Schwingformen der ersten beiden Moden. Bei der dritten Mode zeigt 
sich eine deutliche Abweichung der analytischen Lösung von den experi- 
mentell und numerisch ermittelten Schwingformen. Ursächlich hierfür ist 
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die Vernachlässigung der Schubsteifigkeit des Tragfedermantels durch die 
Modellierung als Platte. Ebenso zeigt sich, dass die Diskrepanz zwischen 
der experimentell und der numerisch ermittelten Schwingform bei den hö- 
heren Eigenschwingformen insbesondere bezüglich der Lage der Schwin- 
gungsknoten zunimmt. Dies ist, wie schon in der Ergebnisdiskussion der 
Bauteilsimulationen erörtert, auf die Inhomogenitäten der Elastomerkör- 
per zurückzuführen. Für die dritte Eigenschwingform zeigt sich, dass die 
Schwingungsbäuche bei & = 0,83 nicht den Maximalwert € = Gmax an- 
nehmen. Ursächlich hierfür ist der angulare Versatz des Laser-Doppler- 
Vibrometers zur Oberflächennormalen des Tragfedermantels, welcher mit 
steigenden &-Werten zunimmt. Die Auswertung der numerischen Schwing- 
formen erfolgt zur besseren Vergleichbarkeit in die gleiche Richtung und 
zeigt daher die gleiche Charakteristik. 


-Messdaten — num. — anal. 


Abbildung 4.15: Darstellung der normierten Betriebsschwingformen[Tre17] 
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4.3.2 Modellierung der Hydraulikeigenschaften 


Wie im vorherigen Abschnitt gezeigt, lassen sich die Eigenschaften des 
Tragkörpers mit Hilfe einer FE-Simulation in einem Frequenzbereich bis 
etwa 2000 Hz abbilden. Um das Verhalten der hydraulischen Lager jedoch 
vollständig abbilden zu können, ist eine Modellbildung der Hydraulikeigen- 
schaften erforderlich. Im Folgenden soll daher zunächst der Einfluss der 
hydraulischen Parameter, der in Abschnitt 3.4.3 experimentell untersucht 
wurde, modelliert werden und diese Modellbildung anschließend mit der 
FE-Simulation gekoppelt werden. 


Modellbildung 


Die in Abschnitt 3.4.3 experimentell erarbeiteten Ergebnisse haben gezeigt, 
dass bei den untersuchten Bauteilen eine Resonanz zwischen ca. 200 und 
500 Hz vorliegt, welche sich durch die Abstimmung der Gittergeometrie 
im Membran-Düsensystem verschieben lässt. Hierbei zeigt sich, dass eine 
zunehmende Gitterhöhe diese Resonanzfrequenz nach unten, eine erhöhte 
Anzahl von Gitteröffnungen die Resonanzfrequenz nach oben verschiebt. 
Dieses Verhalten ist ähnlich zu dem in Abschnitt 2.3.3 diskutierten Re- 
sonanzeffekten im Dämpfungskanal. Es deutet darauf hin, dass mit einer 
zunehmenden Gitterhöhe eine schwingende Masse erhöht wird und durch 
die Vergrößerung der Offenfläche ebenso eine Masse vergrößert wird je- 
doch gleichzeitig auch das hydraulische Übersetzungsverhältnis verändert 
wird. Eine Variation der Tragfedersteifigkeit hat ebenso einen Einfluss auf 
das Schwingsystem (vgl. Abschnitt 3.4.1), was darauf hindeutet, dass das 
für die Hydraulikresonanz notwendige Federelement die Elastomertragfe- 
der darstellt. Auf Basis dieser Ergebnisse kann ein einfaches Ersatzmodell, 
welches diesen Resonanzeffekt abbildet, erstellt werden. 

Beteiligt an diesem Schwingsystem sind dabei eine schwingende Fluidmas- 
se, welche sich in dem Gitter oberhalb und unterhalb der Membran befindet 
und die Blähsteifigkeit der Tragfeder, gegen die diese Masse schwingen 
kann. Zusätzlich hat auch die Membran eine gewisse Federwirkung und die 
dissipativen Effekte durch die hydraulischen Verluste beim Ein- und Aus- 
strömen aus dem Gitter beeinflussen das Schwingsystem. Parallel zu diesem 
System wirkt noch die Grundsteifigkeit der Tragfeder als zusätzlicher Kraft- 
pfad des Systems. Aus diesem Zusammenhang ergibt sich das in Abbildung 
4.16 dargestellte mechanische Ersatzsystem. 


116 


4.3 Modellierung des Bauteilverhaltens von hydr. Lagerelementen 


mit: cr: Komplexe Grundsteifigkeit der Tragfeder; cz: komplexe 
Blähsteifigkeit der Tragfeder; cy: Membransteifigkeit; ky: 
Dämpfungskonstante des Düsen-Membran-Systems; r: hydraulisches 
Übersetzungsverhältnis 


Abbildung 4.16: Darstellung des mechanischen Ersatzsystems zur Abbildung des 
hydraulischen Verhaltens von Aggregatelagern [WGP17] 


Das Ubersetzungsverhiltnis r ergibt sich dabei als Quotient aus der Offen- 
fläche des Gitters und der Pumpfläche der Tragfeder. Die Masse m entspricht 
der Fluidmasse, welche sich in den Gitteröffnungen befindet und errechnet 
sich für dem Fall zylindrischer Gitteröffnungen zu [WGP17]: 


MF uid = AGitter (t + h)PFinid 
T (4.32) 
= "Bohrung 74 Bohrung (t + h)PrFiuid 


r= ATragfeder 


(4.33) 
AGitter 


mit: Bohrung: Anzahl der Gitterbohrungen; dgonrung: Durchmesser der 
Gitterbohrungen; Priud: Dichte des Hydraulikmediums; t: Höhe des 
oberen Gitters; h: Höhe des unteren Gitters; A Gitter: Offenfläche des 
Gitters; ATyagfeder: Pumpfläche der Tragfeder 
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Aus der Betrachtung des Kräftegleichgewichts an der Masse m ergibt sich 
die Differentialgleichung, welche das Schwingsystem beschreibt: 


1 1 
mł» +kyky — -Cp(&ı — —X) = 0 (4.34) 
r r 
1 1 
mx + km&o + xXo(cm + 2) = 72841 (4.35) 


Die sich daraus ergebende Eigenfrequenz des ungedämpften Systems ap 
errechnet sich zu: 


Cp/r tem _ a 
oy = [sel Few _ ja 


Fa (4.36) 


Atrag ATragfeder 
C 
cB + (an AGitter 


AGitter 


ATragfeder 
(p) AGitter (t +h)PFluid 


Und somit die Eigenfrequenz @, des gedämpften Systems zu: 


2 
+ € 
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Gleichungen 4.36 und 4.33 zeigen, dass mit zunehmender Masse die Ei- 
genfrequenz des Systems absinkt. Dies passt zu dem in der Parameterstudie 
gezeigten Verhalten, dass die Resonanzfrequenz mit zunehmender Masse 
absinkt. Durch eine zunehmende Anzahl von Bohrungen wird diese Masse 
ebenfalls erhöht, es ändert sich jedoch auch das hydraulische Übersetzungs- 
verhältnis. Da dieses quadratisch in die Eigenfrequenz eingeht, steigt die 
Eigenfrequenz mit einer zunehmenden Offenfläche an. 


Simulationsergebnisse und -validierung 


Um mit dem vorgestellten Modell das Verhalten der hydraulischen La- 
ger abbilden zu können, ist eine Parametrierung der Variablen notwen- 
dig. Die Masse kann hierbei aus den geometrischen Parametern und der 
Dichte des Hydraulikmediums nach Gleichung 4.32 einfach ermittelt wer- 
den. Die Grundsteifigkeit der Tragfeder und die zugehörige Blähsteifigkeit 
kann durch einfache Versuche ermittelt werden. Für die Grundsteifigkeit 
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reicht es aus, die Tragfeder ohne ein Hydraulikmedium zu charakterisieren. 
Die Blähsteifigkeit lässt sich ermitteln, wenn die Transfersteifigkeit des hy- 
draulischen Lagers mit einem vollständig verschlossenen Düsen-Membran- 
System ermittelt wird. Die Membransteifigkeit hingegen lässt sich expe- 
rimentell nur schwer ermitteln, da hierbei eine Wechselwirkung mit der 
Offenfläche des Gitters besteht. Abhängig von der Offenfläche muss sich 
die Membran aufgrund der Änderung des hydraulischen Übersetzungsver- 
hältnisses mehr oder minder stark in die Öffnungen blähen, was zu einer 
Nichtlinearität der Federsteifigkeit cy im Modell führt. Da sich ebenso der 
Volumenstrom durch die Gitteröffnungen in hohem Maße verändert, wirkt 
sich diese Nichtlinearität auch auf das im Modell enthaltenen Dämpfungs- 
glied aus. Da die experimentelle Ermittlung der Parameter mit sehr hohem 
Aufwand verbunden ist, wird zur Ermittlung der beiden verbleibenden Pa- 
rameter ein Least-Square-Fitting verwendet, sodass das Maximum und der 
Anstieg des Betrages der dynamischen Transfersteifigkeit möglichst gut ab- 
gebildet wird. 

Mit Hilfe der so bestimmten Parameter läst sich aus dem vorgestellten einfa- 
chen Ersatzmodell die Übertragungsfunktion in Form der komplexen Trans- 
fersteifigkeit bestimmen. Der Vergleich zwischen den experimentell ermit- 
telten Werten und den Ergebnissen des Ersatzmodells ist in Abbildung 4.17 
für die experimentell untersuchten Varianten 2 und 4 dargestellt. 

Beim Vergleich der experimentell und aus dem vorgestellten Modell ermit- 
telten Steifigkeiten zeigt sich, dass sich das Verhalten der Hydraulikeigen- 
schaften bis etwa 700Hz sehr genau beschreiben lässt. Die Frequenzlage 
des Maximums sowie die sich ergebenden Werte in Phase und Betrag der 
dynamischen Steifigkeit stimmen sehr gut überein. Die Werte der Phase 
werden im Vergleich zu den experimentell ermittelten Werte durch die Mo- 
dellergebnisse leicht unterschätzt. Insgesamt kann festgestellt werden, dass 
die Modellbildung bis etwa 700 Hz eine hohe Korrelation zu den experimen- 
tellen Ergebnissen aufweist. Durch die Verwendung eines Einfreiheitsgrad- 
Systems und die Anpassung mittels eines Fitting-Verfahrens ist dieses Ver- 
halten zu erwarten, da hierdurch der Resonanzeffekt für eine einzelne Reso- 
nanz sehr exakt nachgebildet werden kann. Für höhere Frequenzen weichen 
die Modellergebnisse jedoch immer weiter von den experimentellen Daten 
ab. Auch der erneute Anstieg der Phase ab etwa 800 Hz, welcher auf einen 
weiteren Resonanzeffekt hindeutet, kann mit dem verwendeten Modell nicht 
abgebildet werden. 
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Abbildung 4.17: Vergleich der experimentell und numerisch ermittelten Trans- 


fersteifigkeit des hydraulischen Verhaltens von Aggregatelagern 
[WGP17] 
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Diskussion 


Der Resonanzeffekt, welcher sich ab etwa 700-800Hz abzuzeichnen be- 
ginnt, kann mit den Ergebnissen aus Abschnitt 4.3.1 auf eine Kontinuums- 
schwingung der Tragfeder zurückgeführt werden. Dieses Verhalten kann oh- 
ne weiteres nicht durch ein einfaches mechanisches Ersatzmodell beschrie- 
ben werden. Das Resonanzverhalten, welches sich bei allen untersuchten 
hydraulischen Lagern gezeigt hat, kann jedoch mit dem aufgezeigten Er- 
satzmodell beschrieben werden. Im folgenden soll daher das zur Beschrei- 
bung der Transfereigenschaften der Elastomerkomponenten verwendete FE- 
Modell erweitert werden, sodass das Verhalten der Hydraulik mit abgebildet 
wird. Hierdurch kann ein Modell erzeugt werden, welches sowohl die Kon- 
tinuumsschwingungen der Elastomeranteile als auch des Schwingverhaltens 
der Hydraulik im hochfrequenten Bereich berücksichtigt. 


4.3.3 Kopplung der Modellierungsverfahren 


Um den Einfluss der hydraulischen Komponenten abzubilden, soll das FE- 
Modell aus Abschnitt 4.3.1 erweitert werden. Damit dies erfolgen kann, 
wird das Düsen-Membran-System der in Abschnitt 3.4 experimentell ana- 
lysierten Bauteile in das FE-Modell mit aufgenommen. Da der lange hy- 
draulische Dämpfungskanal oberhalb seiner Abstimmfrequenz, welche sich 
bei etwa 5-18 Hz befindet, dynamisch abgekoppelt ist (vgl. Abschnitt 2.3), 
wird er bei der Modellierung der hochfrequenzen Transfersteifigkeit nicht 
berücksichtigt. Die nachgiebigen Elemente, welche für das Verhalten der 
Hydraulikeigenschaften im vorherigen Abschnitt identifiziert wurden, sol- 
len jedoch in das Modell integriert werden. Hierzu zählen die nachgiebige 
Membran und die bereits im Modell enthaltene Tragfeder. Auch derjenige 
Anteil des Düsen-Membran-Systems, welcher die Arbeitskammer begrenzt, 
muss mit abgebildet werden, da über diesen die durch den Fluiddruck 
in der Arbeitskammer induzierten Kräfte abgestützt werden. Die Düsen- 
Membran-Systeme bestehen in der Regel aus zwei Einzelteilen, welche 
jeweils die Arbeits- und Ausgleichskammer begrenzen. Für die Modellbil- 
dung genügt derjenige Teil, welcher auf der Seite der Arbeitskammer liegt, 
da von der Ausgleichskammer aus keine Druckkräfte ausgehen, solange der 
Rollbalg als ideal weich angesehen wird. Aufgrund seiner geringen Steifig- 
keit und der Tatsache, dass sich oberhalb der Hydraulikresonanzen nahezu 
kein Volumenstrom mehr durch das Düsen-Membran-System einstellt, wird 
im Folgenden von dieser Annahme ausgegangen. Daher existieren in der 
Ausgleichskammer auch keine Fluidbewegungen mehr. Dies hat zur Folge, 
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dass ebenfalls keine Trägheitskräfte mehr wirken. Diejenigen Trägheitskräf- 
te, welche aus der schwingenden Fluidsäule in den Bohrungen des Hydrau- 
likgitters in der Arbeits- und Ausgleichskammer entstehen (vgl. Abschnitt 
4.3.2), werden durch eine Punktmasse berücksichtigt (dieses Vorgehen wird 
im Folgenden noch detaillierter erläutert). Durch die genannten Vereinfa- 
chungen kann die Ausgleichskammer bei der folgenden Modellbildung ver- 
nachlässigt werden, was die Modellkomplexität deutlich reduziert. 


Modellbildung 


Zur Modellbildung wird das bereits in Abschnitt 4.3.1 verwendete FE- 
Modell der Elastomertragkörper um die Düsenplattenhälfte, welche die Ar- 
beitskammer abschließt, und die zwischen den Gittern befindliche Membran 
erweitert. Der untere Teil des Gitters wird nicht mit in das Modell aufge- 
nommen, da bei den sehr kleinen Belastungsamplituden, welche im aku- 
stisch relevanten Frequenzbereich auftreten, die Auslenkungen der Mem- 
bran so gering sind, dass diese nicht mit dem Gitter in Kontakt kommt. Eine 
Modellierung von Kontaktbedingungen zwischen dem Gitter und der Mem- 
bran ist daher nur an den Randbereichen der Membran notwendig, wo diese 
über eine Dichtkante gegen das Gitter die Arbeitskammer abdichtet. Das 
sich durch diese Modellierung ergebende FE-Modell ist in Abbildung 4.18 
dargestellt. 

Durch die beschriebene Erweiterung des Modells ergibt sich die Arbeits- 


(a): Fluid gefüllte Kavität; (b): Einspannposition der Tragfeder; (c): 
Einspannposition der Düsenplatte; (d): Einspannposition der Membran 


Abbildung 4.18: Darstellung des FE-Modells zur Abbildung des Gesamtsystemver- 
haltens 


kammer des untersuchten hydraulischen Lagers. Durch die Vernachlässi- 
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gung des langen Tilgungskanals stellt die Arbeitskammer ein abgeschlosse- 
nes System dar. Diese ermöglicht die Modellierung einer mit einem Fluid 
gefüllten Kavität, welche die Bewegung der Tragfeder in einen Innendruck 
im hydraulischen System wandelt. Hierdurch kommt es zur Kopplung der 
Tragfederschwingung und der daraus resultierenden Druckkräfte im Sy- 
stem, welche sich auf der Düsenplatte abstützen. Das Fluid wird als kom- 
pressibles, massebehaftetes Medium dargestellt, welches die resultierenden 
Druckkräfte homogen in der gesamten Kavität verteilt. Von dieser homo- 
genen Druckverteilung wird ausgegangen, da die Schallgeschwindigkeit in 
Glykol-Wasser-Gemischen so hoch ist, dass die Wellenlänge einer Druck- 
schwingung deutlich höher ist als die Ausdehnung der Arbeitskammer bzw. 
die erste Eigenfrequenz einer stehenden Welle innerhalb der Arbeitskammer 
eine Größenordnung oberhalb des betrachteten Frequenzbereichs liegt. 

Da das Düsen-Membran-System in einem real existierenden Bauteil kraft- 
schlüssig verbunden ist, erfolgt für die Abbildung im FE-Modell eine fe- 
ste Einspannung der oberen Düsenplatten-Hälfte (Fläche (c) in Abbildung 
4.18). Darüber hinaus wird die untere Dichtkante der Membran, welche 
sich in der Realität auf der unteren Düsenplatten-Hälfte abstützt und über 
diesen Kraftpfad Reaktionskräfte auf dem Gehäuse abstützt, an der Fläche 
(d) in Abbildung 4.18 ebenfalls fest eingespannt. Durch diese Modellierung 
ergeben sich drei Kraftpfade (Frfragfeder, FDiisenplatte und FMembran), welche 
aufsummiert die gesamte Reaktionskraft des Bauteils auf die aufgeprägte 
Elongation am Tragfederkern ergeben. Aus dieser Summenkraft und der be- 
kannten Elongation kann somit die Transfersteifigkeit des Bauteils ermittelt 
werden. 

Das im Lager enthaltene Hydraulikmedium wird durch das kompressible, 
massebehaftete Medium im Modell berücksichtigt, die Strömungsverhält- 
nisse des Mediums in der Arbeitskammer jedoch nicht. Durch diesen Um- 
stand kommt es dazu, dass die hydraulische Übersetzung, welche sich zwi- 
schen dem Gitter und der Pumpfläche der Tragfeder ergibt (vgl. Abschnitt 
4.3.2), mit der Modellierung des Fluids in dieser Form nicht abgebildet wer- 
den kann. 

Um das in Abschnitt 4.3.2 dargestellte Verhalten der hydraulischen Ent- 
kopplung und der damit verbundenen Hydraulikresonanz des Entkoppelgit- 
ters in das FE-Modell integrieren zu können, wird zusätzlich eine Punkt- 
masse in das System eingebracht, welche über eine kinematische Kopplung 
mit der Membran verbunden wird und die mittlere Bewegung der Mem- 
bran mit ausführt. Dieser Massepunkt wird entsprechend des hydraulischen 
Übersetzungsverhältnisses nach den in Abschnitt 4.3.2 dargestellten Zu- 
sammenhängen mit einer Masse beaufschlagt. Durch diese Modellierung 


123 


4 Modellbildung des Übertragungsverhaltens von Aggregatelagern 


kann auch das Resonanzverhalten der Hydraulikentkopplung mit abgebildet 
werden. 


Simulationsergebnisse und -validierung 


Mit Hilfe des im vorherigen Abschnitts beschriebenen Modells kann das 
Gesamtsystemverhalten des untersuchten hydraulischen Topflagers nachge- 
bildet werden. Die Ergebnisse der numerischen Analyse für drei der bereits 
in Abschnitt 4.3.1 und 3.4 untersuchten Tragfedervarianten mit 45 Sh, 53 

Sh und 62 Sh sind in Abbildung 4.19 vergleichend zu den experimentell er- 
mittelten Daten dargestellt. 

Beim Vergleich der experimentell und numerisch ermittelten Daten ist er- 
kennbar, dass sich die bereits in Abschnitt 4.3.2 diskutierte Hydraulikreso- 
nanz durch die Implementierung der hydraulisch wirkenden Masse in das 
FE-Modell sowohl in Betrag als auch in der Phase der dynamischen Trans- 
fersteifigkeit nachbilden lässt. Die Verschiebung dieser Resonanz und der 
Resonanz zwischen 1500 Hz und 2000 Hz mit zunehmender Shore-Härte zu 
höheren Frequenzen ist auch in den numerische ermittelten Daten erkenn- 
bar. Für die Variante mit 62 Sh liegt die zweite Eigenmode des Systems in 
der numerisch ermittelten Kennlinie bei einer etwas höheren Frequenz als 
in den experimentellen Daten. Es ziegt sich, dass die Resonanzüberhöhung 
des Betrages der dynamischen Transfersteifigkeit der zweiten Bauteilmode 
für alle drei dargestellten Varianten höher ausfällt als die experimentell er- 
mittelte. Dies lässt sich auf die dissipativen Effekte des Fluids zurückführen, 
welche sich aufgrund der Modellierung mit nur einem Freiheitsgrad für den 
Innendruck in der Fluidkavität nicht hinreichend abbilden lassen. Trotz der 
Vernachlässigung dieser Dämpfungseffekte zeigen die numerischen Ergeb- 
nisse eine hohe Übereinstimmung mit den experimentell ermittelten Trans- 
fereigenschaften. 

Mit Hilfe des FE-Modells kann durch die hohe Übereinstimmung zwischen 
den Modellergebnissen und den experimentellen Daten im Folgenden eine 
weitere Analyse des Systems und der Wechselwirkungen in diesem statt- 
finden. Hierzu sollen zunächst die Schwingformen des Gesamtsystems be- 
trachtet werden. Diese sind für die beiden Eigenmoden bei etwa 250 Hz und 

1500 Hz in Abbildung 4.20 dargestellt. 

Bei Betrachtung dieser Eigenschschwingformen fällt auf, dass für die er- 
ste Mode eine gekoppelte Schwingung aus einem Aufblähen der Tragfeder 
und einer Blähbewegung der Membran stattfindet. Diese Kopplung deckt 
sich mit den aus der Modellbildung in Abschnitt 4.3.2 gewonnenen Er- 
kenntnissen. Die Fluidsäule in den Entkoppelbohrungen schwingt demzu- 
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—— 45Sh Experimentell - - - 45Sh Numerisch 
—— 53Sh Experimentell - - - 53Sh Numerisch 
—— 62Sh Experimentell - -- 62Sh Numerisch 


Abbildung 4.19: Vergleich der experimentell und numerisch ermittelten Transfer- 
steifigkeit des Gesamtsystemverhaltens von hydraulisch dämpfen- 
den Aggregatelagern 
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(a) Mode 1 


(b) Mode 2 


Abbildung 4.20: Darstellung der die Transfersteifigkeit von hydraulischen Topfla- 
gern beeinflussenden Moden 


folge zwischen der Blähfeder der Membran und der Tragfeder. Nach einer 
dynamischen Abkopplung dieses Schwingers kommt es zu einer quasi ver- 
schlossenen Arbeitskammer, da durch die Entkoppelbohrungen nahezu kein 
Volumenausgleich mehr möglich ist. Demzufolge stellt sich bei der zweiten 
Mode im betrachteten Frequenzbereich eine Eigenschwingform ein, welche 
nur noch durch die Eigenschwingform der Tragfeder beeinflusst wird. Es 
kommt zu der bereits in Abschnitt 4.3.1 diskutierten Biegeschwingung der 
Tragfeder. Diese liegt jedoch bedingt durch die versteifende Wirkung des 
Fluids in der Arbeitskammer gegenüber dem ungefüllten Zustand bei einer 
höheren Frequenz. Die bereits in Abschnitt 4.3.1 festgestellte erste Eigen- 
mode der Tragfeder tritt bei einer Kopplung mit dem hydraulischen System 
daher doppelt auf. Durch die angekoppelte Fluidmasse im Entkoppelgitter 
liegt sie zum einen gegenüber dem ungefüllten System bei einer niedrigeren 
Frequenz. Durch die Abkopplung dieser Masse und die versteifende Wir- 
kung des Fluids bei einer erhöhten Frequenz tritt nochmals die Schwing- 
form mit einer halben Wellenlänge auf der Länge der Gummisehne auf. 
Abbildung 4.21 stellt die experimentell mit Hilfe eines Laser-Doppler- 
Vibrometers ermittelten Schwingformen des hydraulischen Lagers mit ei- 
ner Tragfeder aus einer 45 Sh Elastomermischung dar. Hierbei sind die 
Schwingformen der drei sich im Messbereich befindlichen Bauteilmoden 
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Abbildung 4.21: Darstellung der experimentell ermittelten Schwingformen der 
Tragfedern bei den Eigenfrequenzen bei 90° (links) und —90° 
Phasendifferenz zur Anregung[Tre17] 


bei ca. 225 Hz, 1500 Hz und 2100 Hz dargestellt. Es zeigen sich die gleichen 
Schwingformen ab, welche auch numerisch ermittelt wurden. Für die ersten 
beiden Moden ergeben sich jeweils wieder Schwingungen mit einer halben 
Wellenlänge entlang der Sehnenlänge der Tragfeder. Dies bestätigt die nu- 
merisch ermittelte Aufdopplung dieser Mode durch die Massenankopplung 
der Fluidsäule auf der einen Seite und die versteifende Wirkung des Fluides 
auf der anderen. Die höheren Moden mit Vielfachen von halben Wellen- 
längen auf der Sehnenlänge treten in der numerischen Simulation genauso 
auf. Die Modellqualität ist jedoch auch bedingt durch die in Abschnitt 4.3.1 
analysierten Materialinhomogenitäten für den sehr hohen Frequenzbereich 
eingeschränkt. 

Aus den Ergebnissen des FE-Modells kann jedoch bis etwa 2 kHz eine de- 
tailliertere Analyse des Übertragungsverhaltens erfolgen. Im Gegensatz zur 
experimentellen Analyes kann mit Hilfe der Simulation eine Analyse der 
Einflüsse der einzelnen Kraftpfade durchgeführt werden. Die Reaktionskraft 
des Gesamtsystems setzt sich dabei aus dem Anteil der durch die Tragfeder 
übertragenen Kräfte, der sich auf der Düsenplatte abstützenden Druckkräfte 
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und der sich auf der Membran abstützenden Druckkräfte zusammen. 
Bei Betrachtung der in Abbildung 4.22 dargestellten Einzelkomponenten 
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Abbildung 4.22: Vergleich der Anteile der Reaktionskraft des modellierten hydrau- 
lischen Topflagers auf eine Einheitsauslenkung von Imm 


der Gesamtkraft zeigt sich, dass der Hauptübertragungspfad der Kraft die 
Düsenplatte darstellt. Die sich darüber abstützenden Druckkräfte sind deut- 
lich größer als die übrigen Kräfte. Ebenso zeigt sich, dass über die Mem- 
bran nur eine sehr geringe Kraft übertragen wird, was aufgrund ihrer hohen 
Nachgiebigkeit und der Tatsache, dass dieser Kraftpfad oberhalb der ersten 
Mode dynamisch entkoppelt ist, sinnhaft erscheint. Eine Möglichkeit, die 
Kraftübertragung über die Düsenscheibe zu reduzieren, besteht darin, die 
Düsenscheibe elastisch im Gehäuse des Lagerelementes zu befestigen, um 
eine Entkopplung dieses Kraftpfades zu realisieren (vgl. [Dr.16]). Hierbei 
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besteht jedoch ein Konflikt mit der Forderung nach einer Dichtheit des Sy- 
stems auch bei hohen übertragenen Lasten. Eine weitere Möglichkeit der 
Optimierung bietet sich durch das Einbringen eines definierten Luftvolu- 
mens in die Arbeitskammer. Hierdurch wird die Fluidkompressibilität er- 
höht und die Druckkräfte, welche auf die Düsenplatte wirken, nehmen ab. 
Betrachtungen hierzu sind in Anhang A.3 auf Basis von Simulationen mit 
verschiedenen Fluidkompressibilitäten ausgeführt. 


4.3.4 Fazit 


Die Ergebnisse der Modellbildung zeigen, dass sich das Körperschallüber- 
tragungsverhalten von Buchsenlagerelementen durch eine FE-Modellierung 
mit einem viskoelastischen Materialmodell darstellen lassen. Durch Erwei- 
terung dieses Modells mit einer Fluid-Struktur-Interaktion kann ebenso das 
Verhalten hydraulisch dämpfender Lager mit berücksichtigt werden. Hier- 
bei muss besonderes Augenmerk auf die Modellierung der Überströmquer- 
schnitte von Arbeits- zu Ausgleichskammer gelegt werden (vgl. Abschnitt 
4.3.2). Die Modellgüte kann hierbei durch die hohe Übereinstimmung von 
Simulationsergebnissen mit experimentell erzeugten Messdaten als ausrei- 
chend gut für eine valide Bauteilbeschreibung angenommen werden. Dar- 
über hinaus kann durch die Analyse der Bauteile und die Beschreibung 
durch das FE-Modell festgestellt werden, dass die Eigenschaften der Gum- 
mitragfeder einen sehr großen Einfluss auf das Verhalten eines hydrauli- 
schen Aggregatlagers hat. Hierbei zeigt sich die erste Biegeeigenfrequenz 
der Tragfeder als wichtigster Parameter für das Gesamtsystem. Zur geziel- 
ten Auslegung dieser Komponente, die je nach Fahrzeugrandbedingung ver- 
schiedenen Zielsetzungen unterworfen sein kann, ist es erforderlich, die Pa- 
rameter zu identifizieren, welche das Verhalten beeinflussen. Hierzu soll im 
folgenden Kapitel eine Design-Space-Exploration durchgeführt werden, um 
die Stellgrößen zur Beeinflussung der ersten Biegeeigenfrequenz der Trag- 
körper zu finden. Darüber hinaus soll zusätzlich aufgezeigt werden, in wel- 
chem Umfang die erste Biegeeigenfrequenz bei Einhaltung der aus dem je- 
weiligen Einsatzzweck im Fahrzeug resultierenden Randbedingungen vari- 
iert werden kann. 
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Wie in Abschnitt 4.3.1 und 4.3.3 gezeigt, wird das Körperschallübertra- 
gungsverhalten von hydraulischen Topflagern vor allem durch die Eigen- 
schwingformen der Tragfeder beeinflusst. Auf Basis dieser Ergebnisse soll 
nun im Folgenden untersucht werden, wie sich diese Bauteilmoden gezielt 
beeinflussen lassen. Hierzu sollen die geometrischen Eigenschaften einer 
Tragfeder soweit reduziert werden, dass verschiedene Geometrien automa- 
tisiert erzeugt und auf ihre Eigenfrequenz hin untersucht werden können. 
Anschließend besteht die Möglichkeit, mit Hilfe eines entsprechenden Opti- 
mierungsalgorithmus mögliche Designs auf bestimmte Kriterien hin zu op- 
timieren und auf Basis dieser Ergebnisse Zusammenhänge zwischen den 
Design-Parametern und den Zielwerten und Randbedingungen zu ermitteln. 
Im Folgenden werden diese Schritte detailliert erläutert. 


5.1 Abstraktion der Tragfedergeometrie und 
Modellbildung 


Zur automatisierten Analyse von Tragfedergeometrien wird die in der Praxis 
am häufigsten verwendete, rotationssymmetrische Tragfedergeometrie auf 
sieben wesentliche Design-Parameter heruntergebrochen und als reduzier- 
tes Modell erzeugt. Die Tragfedergeometrie wird dabei als rotationssym- 
metrisch angenommen. Zur Abstimmung der Quer- und Längssteifigkeiten 
kommen auch nicht-rotationssymmetrische Bauteile zum Einsatz, was je- 
doch im Folgenden nicht berücksichtigt wird. Die Kernaussagen, welche 
mit einer rotationssymmetrischen Tragfeder erzeugt werden gelten auch 
für diejenigen Bauteile, welche über verschiedene Grundsteifigkeiten in 
den einzelnen Raumrichtungen verfügen. Der in Abbildung 5.1 dargestell- 
ten, für die automatisierte Erzeugung von Tragfeder-Varianten verwendeten 
Skizze können die sieben wesentlichen geometrischen Design-Parameter 
entnommen werden. Diese sind der Außen- und Innenradius der Tragfeder, 
die Dicke der Elastomersehne an den Schnittstellen zur Periphärie an der 
Außen- und Innenseite (Kern und Außenring), der Aufstellhöhe der Elasto- 
mersehne und die Wölbung auf der Ober- und Unterseite der Elastomerseh- 
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ne. 
Für die Dimensionierung dieser Parameter müssen einige logische und bau- 


r;: Innenradius; rz: Außenradius; h;: Dicke der Elastomersehne innen; ha: 
Dicke der Elsatomersehne außen; hs: Aufstellhöhe der Elastomersehne; Ao: 
Wölbung auf der Oberseite; hu: Wölbung auf der Unterseite 


Abbildung 5.1: Skizze der für die Design-Raum-Erkundung reduzierten Tragfeder- 
geometrie 


raumbedingte Randbedingungen beachtet werden. Die Aufstellhöhe muss 
mindestens so groß gewählt werden, dass bei einer Einfederung der Trag- 
feder bis in den Anschlag die Tragfeder nicht in die Horizontale ausgelenkt 
wird. Hierdurch wird eine Anstellung der Elastomersehne entgegen der 
Lasteinwirkung sichergestellt, was diese vor Zugkräften schützt. Diese sind 
für die Betriebsfestigkeit des Bauteils nachteilig und sorgen für eine ver- 
stärkte Rissausbreitung [Tre15]. Des weiteren sind der Innenradius und der 
Außenradius nicht beliebig wählbar, da zum einen noch ausreichend Bau- 
raum für die Anbindung im Zentrum der Tragfeder, dem Kern, verfügbar 
sein muss, zum anderen auch der Bauraum des Motorlagers im Motorraum 
nicht beliebig vergrößert werden darf. Für die Sehnendicke an der Außen- 
und Innenseite gilt die Randbedingung, dass die Betriebskräfte über die an- 
vulkanisierten Flächen noch sicher übertragen können werden müssen und 
der Bauraum nicht beliebig groß werden darf. Die Wölbung der Tragfeder 
auf der Ober- und Unterseite kann dazu dienen, die Steifigkeitsverältnisse 
der Hoch- zu den Querrichtungen einzustellen und ist auf der Außenseite 
wiederum nur durch den zur Verfügung stehenden Bauraum begrenzt. Auf 
der Innenseite wird die Wölbung durch eine ausreichende Dimensionierung 
der Arbeitskammer begrenzt, da diese auch bei einer vollständigen Einfe- 
derung noch ein gewisses Restvolumen aufweisen muss. Auf Basis dieser 
Randbedingungen ergeben sich für die einzelnen Design-Parameter die in 
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Tabelle 5.1 dargestellten Grenzen. 
Neben den geometrischen Randbedingungen existieren noch weitere aus 


Parameter untere Grenze obere Grenze 


[mm] [mm] 
ri 5 100 
Fa 10 100 
hi 1 100 
ha 1 100 
hs 10 60 
ho 
hy 0 5 


Tabelle 5.1: Grenzen für die Design-Parameter des Tragfedermodells 


der Grundauslegung von Aggregatelagerungen hervorgehende Randbedin- 
gungen. Das Abstützen quasistatischer Gewichts-, Massen- und Betriebs- 
kräfte führt zu Forderungen bezüglich der Grundsteifigkeit in die einzelnen 
Raumrichtungen, welche ebenso in eine Optimierung der Bauteile eingehen 
müssen. 

Zusätzlich bestehen bezüglich der physischen Umsetzbarkeit eines Desi- 
gns noch Nebenbedingungen für die möglichen Tragfedergeometrien. Die 
logische Voraussetzung, dass der Innenradius stets geringer sein muss als 
der Außenradius, ist eine solche Nebenbedingung. Darüber hinaus erfordert 
eine tieffrequente Abstimmbarkeit des Kanals im Düsen-Membran-System 
ein entsprechendes hydraulisches Übersetzungsverhältnis (vgl. Abschnitt 
2.3.4) und bedingt somit die Forderung nach einer Mindestpumpfläche der 
Tragfeder. Diese Nebenbedingungen müssen ebenfalls in das Vorgehen zur 
Optimierung der Tragfedergeometrie mit eingehen. 


5.2 Multikriterielle Optimierung 


Die geometrischen Eigenschaften der Tragfeder eines hydraulischen Aggre- 
gatlagers bestimmen nicht nur dessen Körperschallübertragungsverhalten 
in Form von bauteilspezifischen Eigenmoden sondern auch die statischen 
Grundsteifigkeiten. Eine Optimierung alleinig auf Basis der Transfereigen- 
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schaften hätte zur Folge, dass die Grundsteifigkeiten unter Umständen der- 
art verändert werden, dass das daraus entstehende Aggregatlager nicht mehr 
für den Fahrzeugeinsatz geeignet ist. Um auch diese Anforderungen mit 
in die Optimierungsaufgabe zu integrieren, bedarf es einer multikriteriellen 
Optimierungsstrategie. Ziel dieser Optimierungsaufgabe ist es daher eine 
Kompromisslösung zu finden, bei der die erste Eigenfrequenz der Tragfeder 
möglichst hoch liegt und gleichzeitig die statischen Steifigkeiten des Bau- 
teils in dem für den Fahrzeugeinsatz geeigneten Bereich liegen. Da sowohl 
die statischen Steifigkeiten in die radialen Richtungen und die Hochrich- 
tung sowie die Eigenfrequenz des Bauteils als Zielfunktion zu definieren 
sind, kann ein Vergleich der sich ergebenden Lösungen nicht ohne weiteres 
erfolgen. Zu diesem Ziel eignet sich das Konzept der Pareto-Optimalität. 
Zitzler et al. [ZTL*03] definieren die Pareto-Optimalität eines Lösungs- 
punktes so, dass dieser Lösungspunkt in keinem seiner Zielfunktionswerte 
verbessert werden kann, ohne gleichzeitig einen anderen zu verschlechtern. 
Durch Kenntnis der Lage dieser pareto-optimalen Lösungen kann anschlie- 
Bend der für den jeweiligen Einsatzzweck dienlichste Kompromiss ausge- 
wählt werden. Das Ziel ist es daher, das notwendige Wissen über die Pareto- 
Menge P bereitzustellen. Die Anzahl der möglichen Punkte auf der Pareto- 
Menge ist dabei nicht begrenzt und es können theoretisch unendlich viele 
pareto-optimale Lösungen gefunden werden. Im Falle einer zweidimensio- 
nalen Optimierung stellt die Pareto-Menge eine Pareto-Front dar, welche 
bei mehrdimensionalen Optimierungen einer Hyperfläche entspricht. 

Bei dem im Rahmen dieser Arbeit betrachteten Optimierungsproblem be- 
steht der Design-Raum aus den in Tabelle 5.1 dargestellten Parametern, wel- 
che über die im vorherigen Kapitel erarbeitete FE-Simulation in den Bild- 
bzw. Zielraum überführt werden. Die FE-Simulation entspricht dabei einer 
Abbildung vom Design-Raum in den Bildraum. Der Bildraum bezeichnet 
die Menge der Abbildungen des gesamten Design-Raums. Der Zielraum 
ist eine Teilmenge dieses Bildraums und enthält nur diejenigen Lösungen, 
welche die Rand- und Nebenbedingungen erfüllen. Die Pareto-Menge ist im 
Bildraum stets ein Rand des Zielraums [MGMO06]. Dieser Zusammenhang 
ist in Abbildung 5.2 graphisch dargestellt. 


Die Suche nach dieser Pareto-Front ist das Ziel einer multikriteriellen Op- 
timierung (MOP). Durch die Wahl eines entsprechenden Algorithmus kann 
auf eine vollfaktorielle Analyse des Design-Raumes verzichtet werden (vgl. 
[MGM06]). Die Auswahl einer geeigneten Lösung aus der Pareto-Menge 
kann dabei erst nach der Optimierung der Zielfunktionen erfolgen. Aus die- 
sem Grund spricht Miettinen [Mie08] in diesem Zusammenhang auch von 
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Designraum Bildraum 
x2 a 
zulässige Asian 
Menge 


Pareto-Menge Pareto-Menge 


X1 fi 
(a) Designraum (b) Bild- und Zielraum 


Abbildung 5.2: Darstellung von Design- und Bildraum [Vail6] 


einem a posteriori Optimierungsverfahren. Ein Verfahren, welches vor der 
Abbildung in den Bildraum eine Gewichtung der einzelnen Zielfunktionen 
zueinander vornimmt hingegen entspräche nach Miettinen einem a priori 
Verfahren. Hierbei würde ein Optimierungsalorithmus direkt eine zu präfe- 
rierende Lösung liefern. 

Die a priori Methoden weisen jedoch den entscheidenden Nachteil auf, dass 
durch die Wahl der Gewichte der einzelnen Zielfunktionen die Erkundung 
des Design-Raumes beeinflusst wird und dieser somit unter Umständen 
nicht vollständig ausgeleuchtet wird [KJS11]. Da im Rahmen dieser Arbeit 
das Ziel verfolgt wird, die Einflüsse einzelner geometrischer Parameter der 
Tragfedern auf die Eigenfrequenzen dieser Bauteile zu analysieren, kom- 
men a priori Verfahren nicht in Betracht. Hierzu müsste eine Gewichtung 
oder Priorisierung der einzelnen Ziele anhand eines Beispielauslegungsfal- 
les erfolgen. 


5.2.1 Definition des Optimierungsproblems 


Zunächst soll das mit einem geeigneten Optimierungsalgorithmus zu lösen- 
de Optimierungsproblem definiert werden. Hierzu ist es notwendig, die De- 
signvariablen, die Zielfunktionen und die Nebenbedingungen zu spezifizie- 
ren und anschließend einen geeigneten Algorithmus auszuwählen. 
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Designvariablen 


Diejenigen Designvariablen, welche den Designraum des Optimierungspro- 
blems aufspannen, werden in der Regel zu einem Vektor zusammengefasst. 
Auf Basis der in Abschnitt 5.1 verallgemeinerten Tragfedergeometrie erge- 
ben sich die bereits in Tabelle 5.1 mit ihren Grenzen aufgeführten Design- 
Parameter. Diese lassen sich somit im Designvariablenvektor xpy darstel- 
len: 


TDV = ha (5.1) 


hu 


Da es sich bei den Designvariablen um geometrische Längen handelt, sind 
diese kontinuierlich variierbar. Ein vollfaktorielles Versuchsdesign, bei dem 
der Design-Raum vollständig in den Bildraum überführt wird scheidet damit 
per se für die Analyse aus. 


Zielfunktion 


Die Zielfunktion(en) f eines multikriteriellen Optimierungsproblems wer- 
den so gewählt, dass die Optimerung aus einer Minimierung der Zielfunkti- 
on(en) besteht. Eine Maximierungsaufgabe kann nach Ehrgott [Ehr05] über 
eine entsprechende Invertierung des Vorzeichens in eine Minimierungsauf- 
gabe überführt werden. Die Zielfunktion(en) sind dabei Funktionen des 
Vektors der Design-Variablen. Liegen mehrere Zielfunktionen vor, so wer- 
den diese im Zielfunktionsvektor f zusammengefasst werden. 


F(zov) = (5.2) 


fa (eov) 
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Die bei dem in dieser Arbeit betrachteten Optimierungsproblem verwen- 
deten Zielfunktionen sind zum einen die erste Eigenfrequenz der Tragfeder, 
welche für die Körperschallübertragungseigenschaften des gesamten Lagers 
maßgeblich ist (vgl. Abschnitt 4.3), die statische Bauteilsteifigkeit in die ra- 
diale Richtung der Tragfeder Cstat aq und die statische Bauteilsteifigkeit in 
die axiale Richtung der Tragfeder Cstat,ax. Die erste Eigenfrequenz fEigen 
soll dabei maximiert werden, wobei die statischen Steifigkeiten möglichst 
gering gehalten werden sollen. Eine Erhöhung der statischen Bauteilsteifig- 
keiten kann jederzeit über eine Anpassung der Gummimischung erfolgen, 
was wiederum zu einer Erhöhung der ersten Eigenfrequenz beiträgt. Eine 
möglichst hohe Eigenfrequenz der Tragfeder ist sinnvoll, da mit zunehmen- 
der Frequenz die Karosserieeigenschaften in der Regel bessere Randbedin- 
gungen bereitstellen und somit das durch das Lagerelement übertragene An- 
triebsgeräusch reduziert wird (vgl. Abbildung 2.5Abschnitt in 2.1.2). 

Die Zielfunktionen lauten daher in dem für diese Arbeit betrachteten Opti- 
mierungsproblem: 


fi (pv) E —fEigen (5.3) 
h (pv) = Cstat rad (5.4) 
P (xpv) = Cstat ax (5.5) 


Die Zielfunktion fı ist hierbei das, wie bereits beschrieben, durch Vorzei- 
cheninvertierung erzeugte Minimierungsäquvalent der Maximierungsziel- 
funktion der ersten Eigenfrequenz. 


Nebenbedingungen 


Wie bereits in Abschnitt 5.1 erläutert, existieren für die einzelnen Design- 
Variablen Grenzen, die diese nicht unter- und überschreiten dürfen. Die- 
se werden durch Ungleichungen als Nebenbedingungen in die multikri- 
terielle Optimierung integriert. Es ergeben sich somit für jede beidseitig 
beschränkte Design-Variable zwei Ungleichungen. Sämtliche Ungleichun- 
gen können wiederum in einem Vektor zusammengefasst werden [MGM06, 
Ehr05]. Darüber hinaus existieren gekoppelte Nebenbedingunen der einzel- 
nen Design-Variablen, wie zum Beispiel die Forderung, dass der Außenra- 
dius stets größer als der Innenradius sein muss. Diese können nach Ehrgott 
[Ehr05] in einem Ungleichungssystem der Form 


Azpy <b (5.6) 
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dargestellt werden. Die Matrix A wird dabei so gewählt, dass die entspre- 
chend gekoppelten Design-Variablen berücksichtigt werden und die jewei- 
ligen Grenzen im Vektor b enthalten sind. Diejenigen Punkte des Design- 
Raumes, welche die Nebenbedingungen erfüllen, werden als zulässige Men- 
ge x bezeichnet [Vail6]. 


5.2.2 Beschreibung des zum Lösen des 
Optimierungsproblems verwendeten Algorithmus 


Zur Lösung des im vorherigen Abschnitt beschriebenen Optimierungspro- 
blems soll ein Partikel-Schwarm-Optimierungsalgorithmus (PSO) einge- 
setzt werden. Dieser Optimierungsalgorithmus wurde erstmals von Kenne- 
dy und Eberhart [EK95] beschrieben und basiert auf einem Schwarm von 
Partikeln, welche den Suchraum erkunden, wobei jeder Partikel von den 
Verbesserungen profitiert, welche die anderen dabei entdecken. Eine PSO 
ahmt das Verhalten von Vogelschwärmen in der Natur nach, welche auf 
der Suche nach einem optimalen Landeplatz sind. Hierzu werden zunächst 
einzelne Partikel generiert und an zufälligen Stellen im Design-Raum ver- 
teilt, welche sich in zufällige Richtungen des Design-Raumes bewegen. 
Die Anzahl der Dimensionen j des Positions- und Richtungsvektors jedes 
Partikels entspricht dabei der Anzahl an Design-Variablen. Die Bewegungs- 
richtung jedes Partikels wird dann schrittweise so verändert, dass dieses 
sich in Richtung der besten bekannten Position von sich selbst und allen 
anderen Partikeln bewegt. Auf diese Weise erkundet jeder Partikel seine 
unmittelbare Umgebung und findet eventuell eine neue beste Position und 
beeinflusst damit die anderen Partikel im Schwarm [PC10]. Die Position 
und Geschwindigkeit jedes Partikels x € R/ und v € Rİ werden dabei in je- 
dem Optimierungsschritt nach zwei Formeln aktualisiert. Diese beiden sind 
nach Shi und Eberhart für die Berechnung des k-ten Berechnungsschrittes 
[ESOO]: 


VE OVk-1 + Opp (MP — Le-1) + Vere (Gg — LK-1) (3.7) 
TR Le-1 + VUE (5.8) 


Hierbei bestimmt œ, wie groß der Anteil der bisherigen Geschwindigkeit an 
der neuen Geschwindigkeit ist, und wird daher auch als Trägheitsgewicht 
bezeichnet. Der Einfluss der bisher bekannten besten Position p des betrach- 
teten Partikels und des gesamten Schwarmes g wird mit den sogenannten 
Verhaltensparametern @, und g gesteuert, welche bei der Definition des 
Optimierungsalgorithmus gewählt werden. Hierdurch kann die Konvergenz 
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des Algorithmus gesteuert werden. Eine zu schnelle Konvergenz kann dabei 
unter Umständen dafür sorgen, dass sich die Partikel nur einem lokalen Ma- 
ximum nähern und das globale Maximum nicht gefunden wird. Die Parame- 
ter rp und rg sind zur Vermeidung von zu schneller Konvergenz daher sto- 
chastische Parameter (r,,r, € [0,1]) und werden bei jedem Optimierungs- 
schritt zufällig bestimmt. In jedem Optimierungsschritt wird anschließend 
die neue Position x; bestimmt, indem zur letzten Position die Geschwin- 
digkeit des Partikels addiert wird (vgl. Formel 5.8) [PC10]. Sollte der Par- 
tikel durch seine neue Position die Grenzen des Design-Raumes verlassen, 
so wird dieser wieder in den Design-Raum bewegt. Darüber hinaus kann es 
sinnvoll sein, den Betrag der Geschwindigkeit zu begrenzen, sodass ein Par- 
tikel sich nicht zu schnell im Design-Raum bewegt. Die Partikel-Schwarm- 
Optimierung ergibt sich demnach zu dem in Abbildung 5.3 in Pseudo-Code 
dargestellten Algorithmus. 

Der Vorteil dieses Verfahrens ist, dass es nicht gradientenbasiert ist und so- 
mit die Differenziation der Abbildungsfunktion, welche den Design- in den 
Bildraum überführt, nicht berechnet werden muss [PC10]. Dies stellt vor al- 
lem dann einen entscheidenden Vorteil dar, wenn der Gradient dieser Funk- 
tion nicht ermittelt werden kann. Ein weiterer Vorteil ist die Parallelisierbar- 
keit der Abbildungsfunktion, welche in einem Optimierungsschritt für jeden 
Partikel ausgewertet werden muss. 


5.2.3 Diskussion der Ergebnisse 


Zur Durchführung des Design-Space-Exploration wird der im vorherigen 
Abschnitt vorgestellte Optimierungsalgorithmus mit dem FE-Modell der 
Tragfeder verknüpft und eine automatisierte Auswertung der Abbildungs- 
funktion durch die FE-Berechnung implementiert. Die Rechenzeit für die 
Auswertung der Abbildungsfunktion beträgt nur wenige Sekunden und ein 
manuelles Eingreifen ist nicht erforderlich. Somit ist eine große Zahl von 
Optimierungsschritten in einer relativ kurzen Zeit darstellbar. Hierdurch 
kann eine große Zahl von Punkten im Design- und Bildraum entstehen. Zur 
Ergebnisdarstellung werden nur diejenigen Datenpunkte betrachtet, welche 
über eine statische, axiale Steifigkeit zwischen 215 N/mm und 255 N/mm, 
was für eine Auslegung auf 235 N/mm einer Toleranz von circa 10 % ent- 
spricht. Für die Simulationen wurde das Materialmodell mit 45 Shore ver- 
wendet, sodass die sich ergebenden Steifigkeiten noch über die Shore-Härte 
variiert werden könnten. Hierdurch würden jedoch ebenfalls die Eigenmo- 
den der Tragfedern beeinflusst werden. 

Abbildung 5.4 a) und b) zeigen die Lage der erzeugten Datenpunkte als 
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Eigenfrequenz in Hz 


Eigenfrequenz in Hz 


Eigenfrequenz in Hz 
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Abbildung 5.4: Ergebnisse der Optimierung 
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Abbildung 5.3 : Darstellung des Partikel-Schwarm- 
Optimierungsalgorithmus in Pseudo-Code nach [PC10] 
i=0; 
Initialisiere eine Anzahl n von Partikeln mit zufalligen Positionen 
und Geschwindigkeiten im Design-Raum; 
while i < imax do 
if Abbruchkriterium erfüllt then 
Optimierung Abbrechen 
else 
foreach Partikel mit Position x und Geschwindigkeit v do 
Aktualisiere die Geschwindigkeit des Partikels nach 
Gleichung 5.7; 
if v > Vmax then 
| Begrenze Geschwindigkeit des Partikels; 
Bewege Partikel nach Gleichung 5.8 an neue Position; 
if Partikel außerhalb der Design-Raum-Grenzen then 
| Bewege Partikel zurück auf den Rand; 
Bestimme Zielfunktion f(a) des Partikels x; 
if f(a) < f(p) then 
Aktualisiere beste bisherige Position des Partikels 
P & T; 
if f(x) < f(g) then 
Aktualisiere beste bisherige Position des gesamten 
Schwarms p - 2; 


gefärbte Flächen im Bildraum. Hierin ist ebenfalls der Verlauf der Pareto- 
Front erkennbar. Abbildung 5.4 a) stellt die erste Eigenfrequenz über der 
statischen, axialen Steifigkeit der erzeugten Tragfedergeometrien dar. Hier- 
bei kann festgestellt werden, dass die erste Eigenfrequenz mit einer zu- 
nehmenden statischen Axialsteifigkeit zunimmt. Es fällt dabei auf, dass ein 
nichtlinearer Zusammenhang besteht und für eine Axialsteifigkeit zwischen 
210 und 240 N/mm nur ein geringer Anstieg der ersten Eigenfrequenz er- 
zielbar ist. Eine deutliche Zunahme der Eigenfrequenz ergibt sich für axiale 
Steifigkeiten oberhalb von 240 N/mm. Ein ähnlicher Zusammenhang ergibt 
sich für die radialen Steifigkeiten in Abbildung 5.4 b). Auch hier ergibt sich 
ein signifikanter Anstieg der Eigenfrequenzen erst für Steifigkeiten oberhalb 
von 250 N/mm. Die höchste erste Eigenfrequenz, welche durch die Optimie- 
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rung ermittelt werden kann beträgt 816 Hz, was jeweils bei der durch die 
Randbedingung maximal zulässigen statischen Steifigkeit einhergeht. Die- 
ser Zusammenhang ist insofern plausibel, als dass eine erhöhte Steifigkeit in 
der Regel mit einer höheren ersten Eigenfrequenz eines Bauteils einhergeht. 
Vergleicht man jedoch die Ergebnisse der PSO mit den experimentellen 
Ergebnisse für die Tragfederuntersuchungen aus Abschnitt 4.3.1, so liegt 
die erste Eigenfrequenz des Bauteils mit 43 Sh und einer statischen Grund- 
steifigkeit von 250 N/mm in axialer und radialer Richtung bei ca. 600 Hz. 
Demzufolge kann bezüglich dem Stand der Technik entsprechenden Bautei- 
len eine Steigerung der ersten Eigenfrequenz um mehr als 200 Hz erreicht 
werden, was einer Erhöhung um mehr als 30 Prozent entspricht. 

Abbildung 5.4 c)-g) stellt die erste Eigenfrequenz über den Design-Variablen 
bzw. geometrischen Größen der jeweiligen Tragfedergeometrien dar. In Ab- 
bildung 5.4 kann erkannt werden, dass ein quadratisch-reziproker Zusam- 
menhang zwischen der Sehnenlänge / (J = / (ra — ri)? +h?) und der ersten 
Eigenfrequenz besteht. Dieser Zusammenhang kann nach Gleichung (4.31) 
aus der Vergleichsrechnung mithilfe der Kirchhoff’schen Plattentheorie in 
Abschnitt 4.3.1 erwartet werden. Ebenso zeigt sich, dass eine Verringerung 
der Sehnenlänge unter circa 40 mm zu keiner weiteren Zunahme der Eigen- 
frequenz mehr führt, da die Einhaltung der Steifigkeitsgrenzen über eine 
entsprechende Dickeänderung keine Erhöhung der Eigenfrequenz mehr zu- 
lässt. Zur Interpretation und Diskussion sind in den Abbildungen 5.4 die 
Punkte A bis E eingezeichnet, um die Wechselwirkungen der einzelnen Pa- 
rameter zwischeneinander darzustellten. Bei Betrachtung der Punkte A und 
B zeigt sich eine sehr unterschiedliche Sehnenlänge jedoch die gleichen 
statischen Steifigkeiten in axialer und radialer Richtung und eine nahezu 
identische Eigenfrequenz. Hier muss zur Einhaltung der Grenzwerte für die 
geforderten Steifigkeiten eine Änderung der Dicke der Gummisehne erfol- 
gen, was durch die Lager der Punkte A und B in Abbildung 5.4 g) erkannt 
werden kann. 

Abbildung 5.4 d) kann entnommen werden, dass sich eine negative annä- 
hernd lineare Korrelation zwischen der ersten Eigenfreuenz und der Auf- 
stellöhe einstellt. Dies resultiert aus dem linearen Zusammenhang zwischen 
der Sehnenlänge und der Aufstellhöhe. 

Bei Betrachtung von Abbildung 5.4 e), in welcher die Eigenfrequenz über 
dem Innenradius aufgetragen ist, zeigen sich drei, die Menge des Bildrau- 
mes beschränkende Grenzen. Die erste dieser drei verläuft zwischen den 
Punkten E und D und wird durch die untere Grenze der radialen Steifigkeit 
bestimmt. Dies bedeutet, dass es nicht möglich ist, mit geringen Innenradien 
eine hohe Eigenfrequenz zu erreichen, ohne die geforderte minimale radia- 
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le Steifigkeit von 220 N/mm zu unterschreiten. Die zweite Grenze verläuft 
zwischen den Punkten A bzw. B und dem Punkt D, diese wird durch die 
Pareto-Front aus Abbildung 5.4 a) und b) beschränkt. Die dritte Grenze für 
Innenradien, welche größer als 50mm sind, wird durch die obere Steifig- 
keitsgrenze von 270 N/mm in radialer und 255 N/mm in axialer Richtung 
bestimmt. 

Ein ähnliches Bild ergibt sich für die Abhängigkeit der Eigenfrequenz vom 
Außenradius, welche Abbildung 5.4 f) entnommen werden kann. Für Au- 
Benradien kleiner als ca. 35 mm ist es nicht möglich die axialen Steifig- 
keitsvorgaben zu erreichen, was aus der Lage des Punktes E geschlossen 
werden kann. Die radialen Steifigkeitsvorgaben begrenzen die Fläche dann 
in etwa zwischen den Punkten E und D (vgl. Lage der Punkte in Abbil- 
dung 5.4 b)). Anschließend begrenzt die maximale Steifigkeitsvorgabe in 
axialer Richtung den Bildraum an der Grenze, auf der die Punkte A und B 
liegen. Die Begrenzung zwischen den Punkten A und C kann entweder auf 
eine Kombination aus Randbedingungsvorgaben oder auf eine unzureichen- 
de Exploration des Design-Raumes abseits der Pareto-Front zurückgeführt 
werden. Da der Punkt C eine große Entfernung zu dieser aufweist und die 
PSO das Ziel hat, möglichst schnell gegen gute Lösungen mit maximaler 
Eigenfrequenz zu konvergieren, wird der Design-Raum in großem Abstand 
zur Pareto-Front unter Umständen nicht vollständig erkundet. 

Abbildung 5.4 g) zeigt die Eigenfrequenz über der Aufdickung d (d = 
hy +ho) der untersuchten Tragfedern. Hierin ist zu erkennen, dass die Eigen- 
frequenz nahezu unabhängig von der Aufdickung der Tragfeder eingestellt 
werden kann, lediglich für sehr große Aufdickungen größer als 7 mm nimmt 
die Eigenfrequenz wieder ab. An der Lage des Punktes E kann erkannt wer- 
den, dass zur Einhaltung der Steifigkeitsgrenzen die Aufdickung bei sehr 
langen Sehnen und niedriger Aufstellhöhe erhöht werden muss, was per se 
mit einer geringeren ersten Eigenfrequenz einhergeht. 


5.2.4 Fazit 


Zusammenfassend kann geschlussfolgert werden, dass die erste Eigenfre- 
quenz des Tragkörpers über die geometrischen Design-Parameter stark be- 
einflusst werden kann und zur Maximierung der Biegeeigenfrequenz kom- 
pakte Tragkörpergeometrien zu bevorzugen sind, welche eine kurze Sehnen- 
länge haben und eine geringe Aufstellhöhe. Ebenso kann geschlussfolgert 
werden, dass erhöhte Steifigkeiten in die axiale und radiale Richtung förder- 
lich sind für eine Erhöhung der Eigenfregeunz. Dies kann - sofern bei der 
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Grundauslegung der Lagerung berücksichtigt — genutzt werden um die Mo- 
den des hydraulischen Topflagers auf die Erfordernisse seitens der Karosse- 
rie oder der Motoranregung eines Fahrzeuges anzupassen. Ebenso kann aus 
der Design-Raum-Analyse ein Optimierungspotential bezüglich dem Stand 
der Technik aufgezeigt werden, welches theoretisch bis zu 30 % beträgt. 
Dieses Potential wird sich aus Gründen der Betriebsfestigkeit der Bautei- 
le nicht voll ausschöpfen lassen, lässt aber den Schluss zu, dass bezüglich 
einer Eigenfrequenzabstimmung der Bauteile noch Verbesserung möglich 
ist. Wie in Abschnitt 2.1.2 aufgezeigt, ist jedoch eine reine Fokussierung 
auf die Maximierung der Eigenfrequenz nicht zweckdienlich, da auf diese 
Weise unter Umständen eine Empfindlichkeit der Karosseriestruktur mit ei- 
ner Eigenfrequenz des Motorlagers zusammenfällt. Eine Abstimmung der 
Systeme aufeinander und eine gezielte Auslegung des Motorlagers in einen 
Frequenzkorridor, in welchem sowohl die Motoranregung niedrig und die 
Karosserieempfindlichkeiten (NTF und Eingangssteifigkeit) günstig sind, ist 
hierzu ein geeignetes Mittel. Die verwendete PSO kann bei Kenntnis dieser 
Randbedingungen auch dazu genutzt werden die Eigenfrequenz gezielt ein- 
zustellen. Durch die Ergebnisse der in der vorliegenden Arbeit durchgeführ- 
ten Design-Raum-Exploration kann jedoch ein Abstimmfenster aufgezeigt 
werden, in welchem sich Tragkörper hydraulischer Topflager bewegen kön- 
nen. 
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6.1 Zusammenfassung der Arbeitsschwerpunkte 


In der vorliegenden Arbeit wurde zunächst eine Methode aufgezeigt, den 
Körperschallübertragungsweg von antriebserregten Schwingungen mit Hil- 
fe mechanischer Vierpole zu modellieren. Aus dieser Modellierung konnten 
die wesentlichen Eigenschaften von Aggregatelagerelementen herausgear- 
beitet werden und eine experimentelle Methode zur Bestimmung dieser vor- 
gestellt werden. Mit Hilfe dieser Methode wurden Parameterstudien an kon- 
ventionellen und hydraulischen Lagerelementen durchgeführt und die we- 
sentlichen Einflussparameter auf das Übertragungsverhalten identifiziert. In 
einem nachfolgenden Schritt wurde eine Modellierung der Bauteile mit Hil- 
fe finiter Elemente durchgeführt und auf diese Weise ein Systemverständnis 
geschaffen. Aufbauend auf dieser Modellierung wurden im Rahmen einer 
Design-Space-Exploration die geometrischen Parameter der Elastomerkom- 
ponente von hydraulischen Topflagern zum gezielten Systemdesign ermit- 
telt. 


6.2 Wissenschaftlicher Beitrag 


Aufbauend auf die Arbeit von Sell [Sel05] wurden in der vorliegenden Ar- 
beit die Körperschalltransfereigenschaften von elastischen Elementen im 
akustisch relevanten Frequenzbereich analysiert. Die vorliegende Arbeit be- 
schäftigt sich dabei im speziellen mit hydraulischen und konventionellen 
Aggregatelagerelementen und schließt dabei nicht nur eine experimentelle 
Analyse mit ein, sondern führt darauf aufbauend auch eine Modellbildung 
der Bauteile durch. Die Betrachtung dieser Bauteile im Kilohertzbereich 
und eine anschließende Modellbildung findet in keiner wissenschaftlichen 
Arbeit im Bereich der Körperschallübertragung in der Fahrzeugtechnik bis- 
her Beachtung. 

Bei der Modellbildung wird im Besonderen der Einfluss der Interaktion zwi- 
schen Hydraulikmedium und den Elastomerkomponenten auf die Bauteilre- 
sonanzen untersucht. Diese Interaktion konnte mithilfe eines FE-Modells 
und einem viskoelastischen Materialmodell für die für Aggregatelager üb- 
licherweise verwendeten Naturkautschukelastomere abgebildet werden. In 
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keiner wissenschaftlichen Arbeit wird dabei bisher diese Fluid-Struktur- 
Interaktion in einer Modellbildung berücksichtigt. Durch diese Modellbil- 
dung konnten direkt die Gestaltungsparameter der Hydraulikeinheit identifi- 
ziert werden, mit denen sich das Bauteilverhalten beeinflussen lässt. Um die- 
se Design-Paramter auch für die Elastomer-Tragfedern zu ermitteln, wurde 
im nachfolgenden Schritt eine Design-Space-Exploration durchgeführt. So 
konnten konkrete Design-Richtlinien erarbeitet werden, mit denen sich das 
Bauteilverhalten in eine entsprechende Richtung beeinflussen lässt. Durch 
die Kenntnis dieser Parameter und der Verwendung des in der vorliegen- 
den Dissertation erarbeiteten Bauteilmodells kann bei der Neuauslegung 
von Bauteilen das Körperschallübertragungsverhalten bereits in einer digi- 
talen Phase mit berücksichtigt werden. Darüber hinaus kann auf akustische 
Phänomene, welche bei der Triebstrangintegration ins Fahrzeug auftreten, 
direkt reagiert werden und über Detailänderungen an den einzelnen Kom- 
ponenten Resonanzen des Aggregatelagers verschoben werden. 


6.3 Ausblick 


Die im Rahmen dieser Arbeit vorgestellte experimentelle Methode zur Mo- 
dellierung des Transferpfades und zur Ermittlung der Körperschalltransfe- 
reigenschaften kann nicht nur für Aggregatelagerelemente, sondern für ver- 
schiedene schwingungstechnische Anforderungen verwendet werden. Am 
Beispiel eines Kraftfahrzeuges kann dies zum Beispiel auch die Körper- 
schallweiterleitung von Reifen-Fahrbahngeräuschen über die Radführung 
sein. Die finite Elemente Modellierung von Elastomerbauteilen, mit der ge- 
koppelten Fluid-Struktur-Simulation kann auch für weitere Bauteile, wie 
z.B. hydraulische oder konventionelle Lenkerlager verwendet werden. Mit 
einer gezielten Auslegung dieser Bauteile kann die Körperschallübertra- 
gung über die Radführung ebenfalls gezielt ausgelegt werden. 

Die Ergebnisse der Design-Space-Exploration können verwendet werden, 
um eine gezielte Auslegung der Eigenfrequenzen der untersuchten Aggre- 
gatelager auf die Anforderungen der Fahrzeugumgebung, insbesondere die 
Körperschalltransfereigenschaften der Karosseriestruktur, anzupassen. Auf 
diese Weise kann ein gezieltes Design des per Körperschall übertragenen 
Antriebsgeräusches vorgenommen werden. Hierdurch kann zum einen ein 
Sound-Cleaning betrieben werden, zum anderen aber auch ein gezieltes 
Sound-Design vorgenommen werden, um nur einen gewissen Frequenzin- 
halt des Antriebsgeräusches am Insassenohr hörbar zu machen. 
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A Anhang 


A.1 Erläuterung zum direkten Prüfverfahren zur 
Bestimmung von Transfersteifigkeiten 


1:Schwingerreger; 2:höhenverstellbare Traverse; 3:Prüfrahmen; 
4:Kraftmesssystem (eingangsseitig); 5:Prüfkörper; 6:Kraftmesssystem 
(ausgangsseitig); 7:starre Grundplatte 


Abbildung A.1: schematischer Aufbau eines Prüfstandes nach dem direkten Ver- 
fahren [DIN11b] 


Beim Messprinzip nach dem direkten Verfahren kommt ein möglichst stei- 
fer Prüfrahmen zum Einsatz. Dieser Prüfrahmen hat zur Aufgabe, die aus- 
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gangsseitige Bewegung des Prüfkörpers möglichst gering zu halten. Der 
schematische Aufbau eines Prüfstandes nach diesem Verfahren ist in Ab- 
bildung A.1 dargestellt. Der Schwingerreger (1) ist in dem in Abbildung 
A.1 dargestellte Aufbau in einer Quertraverse (2) angebracht. In der Regel 
kommen bei Prüfständen nach dem direkten Verfahren servohydraulische 
Schwingerreger zum Einsatz. Durch diesen Schwingerreger wird eine De- 
formation auf den Prüfkörper (5) aufgebracht. Durch die Anbindung der 
Ausgangsseite des Prüfkörpers an die Grundplatte des Prüfstandes (7) wird 
dessen Schwingweg x2 sehr klein und eine Berechnung der Eingangs- und 
Transfersteifigkeit nach Gleichung 3.11 und 3.12 wird ermöglicht. 

Da dieses Verfahren nur in einem Frequenzbereich valide Ergebnisse liefert, 
in dem der Prüfrahmen und die Traverse als starres Element wirken, kann ei- 
ne Bauteilprüfung nur unterhalb der Eigenfrequenz der ersten Betriebsmode 
des Prüfrahmens erfolgen. Diese Eigenfrequenzen liegen in der Regel in in 
einem Frequenzbereich von einigen hundert Hertz. Um eine experimentelle 
Analyse der Bauteile im akustisch relevanten Frequenzbereich durchführen 
zu können, wird daher für die experimentellen Analyse in dieser Arbeit auf 
das indirekte Messverfahren nach [DINO6] zurückgegriffen. 
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A.2 Analyse der Härte-Streuung in den Elastomerbauteilen 


Abbildung A.2: Darstellung der Koordinatensysteme und Messbereiche der zur 
Analyse der Härte-Streuung innerhalb der Bauteile verwendeten 
Schnitte 
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A.2 Analyse der Härte-Streuung in den 
Elastomerbauteilen 


O 6.25 125 18.75 25 31.25 37.5 43.75 50 
Umfangskoordinate @ in [°] 


O 6.25 125 18.75 25 31.25 37.5 43.75 50 
Umfangskoordinate @ in [°] 


Härte [Sh] 
MANA 
CORMNMWAMN 


0 6.25 125 18.75 25 31.25 37.5 43.75 50 
Umfangskoordinate @ in [°] 


O 6.25 125 18.75 25 31.25 37.5 43.75 50 
Umfangskoordinate @ in [°] 
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e Messpunkte Außen * Messpunkte Innen 
____ Fit Umfang Außen Fit Umfang Innen 
Tragfedern mit der nominellen Shore-Härte: a): 43 Sh; b): 45 Sh; c): 50 Sh; 
d): 55 Sh; e): 60 Sh 


Abbildung A.3: Shore-Härten-Verteilung auf der Mantelhorizontalen der Tragfe- 
dern auf der Innen- und Außenseite 
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A.3 Numerische Untersuchungen zum Einfluss der 
Fluidkompressibilität in hydraulischen 
Lagerelementen 


Bei der Betrachtung der erhöhten sich auf der Düsenplatte abstützenden Re- 
aktionskräfte der ersten Biegemode der Tragfeder in einem hydraulischen 
Topflager in Abschnitt 4.3.3 ist als ein Lösungsansatz zur Reduktion die- 
ser Überhöhung eine Reduktion der Fluidkompressibilität aufgezeigt. Zur 
näheren Untersuchung des Einflusses des Kompressionsmoduls des Hy- 
draulikmediums sollen auf Basis des im Rahmen dieser Arbeit erarbeiteten 
FE-Modells Prinziprechnungen angestellt werden. Hierzu wird der Kom- 
pressionsmodul des Fluids auf 70 % und 50 % des ursprünglichen Wertes 
reduziert. Die Simulationsergebnisse sind in Abbildung A.4 dargestellt. 
Beim Vergleich der drei dargestellten Varianten wird ersichtlich, dass sich 
die Überhöhung der dynamischen Transfersteifigkeit, welche durch die 
Kontimuumsschwingung der Tragfeder hervorgerufen wird, mit einer Ab- 
senkung des Kompressionsmoduls reduziert. Die Nachgiebigkeit des Flui- 
des wird hierdurch erhöht und die Reaktionskräfte auf die Düsenplatte neh- 
men ab. Ebenso kann festgestellt werden, dass die Eigenfrequenz der ersten 
Kontinuumsmode der Tragfeder mit einer Reduktion des Kompressions- 
moduls abnimmt. Dies kann durch die geringere versteifende Wirkung des 
Fluids bei einem geringeren Kompressionsmodul begründet werden. 

Als Ergebnis kann daher festgehalten werden, dass eine Reduktion des 
Fluidkompressionsmoduls sich günstig auf die Transfersteifigkeiten von 
hydraulischen Topflagern auswirkt. Für eine technische Umsetzung die- 
ses Prinzips muss jedoch auf die Wechselwirkung mit der tieffrequenten 
Abstimmung zur Bedämpfung von Aggregatemoden geachtet werden. Eine 
Absenkung des Komressionsmoduls zieht zwangsläufig eine Reduktion der 
Dämpfungswirkung im tieffrequenten Bereich nach sich. Die Umsetzung 
zur Reduktion kann zum einen über eine Substitution des Hydraulikmedi- 
ums erreicht werden oder über das Einbringen einer definierten Kompres- 
sibilität in die Arbeitskammer in Form von Lufteinschlüssen. Bei letzterer 
muss darauf geachtet werden, dass die Nachgiebigkeit in der Arbeitskam- 
mer verbleibt und nicht im Laufe des Betriebes in die Ausgleichskammer 
wandert. Eine Lösung hierzu ist durch den Autors dieser Arbeit in Form 
eines Patentes geschützt. (vgl. [Dr.17]) 
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Abbildung A.4: numerisch ermittelte dynamische Transfersteifigkeit mit verschie- 
denen Fluidkompressionsmoduli 
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A.4 Ausgewählte Datenpunkte aus den Ergebnisse der Design-Space-Exploration 
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Im Zuge der Emissionsreduktion eines Kraftfahrzeuges werden sowohl innermo- 
torische, wirkungsgradsteigernde Maßnahmen ergriffen als auch die Triebsträn- 
ge zunehmend elektrifiziert. Dies führt zu einem verstärkten Schwingungseintrag 
bei hohen Frequenzen und bedingt eine Verbesserung des Isolationsverhaltens 
des Körperschallpfades vom Antrieb bis zum Fahrerohr. Hierzu ist es notwendig, 
die Übertragungseigenschaften der Lagerelemente auf die im Fahrzeug vorherr- 
schenden Randbedingungen anzupassen. Dies betrifft insbesondere die Körper- 
schallempfindlichkeiten der Karosserie und das Anregungsspektrum des Antriebs. 


In dieser Arbeit wird daher das Körperschallübertragungsverhalten von Aggre- 
gatelagern im akustisch relevanten Frequenzbereich experimentell und numerisch 
untersucht. Hierzu wird zunächst ein Prüfstand zur Ermittlung der dynamischen 
Körperschalltransfereigenschaften vorgestellt und anschließend die wichtigsten 
Einflussparameter auf das Bauteilverhalten untersucht. In einer anschließenden 
Modellbildung wird das Systemverhalten durch ein FE-Modell abgebildet. 
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